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Summary 



The progenitors of core-collapse supernovae are stars with an initial mass greater than about 8M Q . Un- 
derstanding the evolution of these stars is necessary to comprehend the evolution and differences between 
supernovae. 

We have constructed new and unique opacity tables to increase model accuracy during the latest stages 
of stellar evolution. We have investigated how initial mass, initial composition and mass loss affects the 
progenitors and their populations. There are many prescriptions for mass loss. Different research groups 
use their preferred rates. We have compared 12 different prescriptions and determined which provides the 
best fit to observations. We use our preferred mass-loss scheme to make suggestions as to the source of the 
differences between supernova types from our progenitor models. 

Binary evolution is considered in order to search for low luminosity SN progenitors and progenitor types not 
possible from single stars. Removal of the hydrogen envelope is more common and we find quite different 
hydrogen deficient SN progenitors. We discuss the implications of our binary models for ultra-luminous 
X-ray sources and gamma-ray bursts. We present an estimation of the mass distribution for black holes at 
various metallcities showing that massive black holes are not formed until very low metallicities. Finally we 
combine the single star and binary results to determine their relative populations and compare to observa- 
tions. However it is not possible to draw many firm conclusions because of the uncertainty in observations 
to date. 
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Chapter 1 



Introduction 



1.1 Looking at stars 



Pumbaa: "Timon, ever wonder what those sparkly dots are up there?" 

Timon: "Pumbaa, I don't wonder; I know." 

Pumbaa: "Oh. What are they?" 

Timon: "They're fireflies that got stuck up in that bluish-black thing". 

Pumbaa: "Oh, gee. I always thought they were balls of gas burning billions of miles away." 

Timon: "Pumbaa, with you, everything's gas." 

The Lion King. 



Everything starts with gas. Slowly over time diffuse interstellar gas collapses under its own gravity, forming 
a rough sphere that gradually contracts and heats up. Eventually the centre becomes hot enough to ignite 
nuclear reactions and forms a gravitationally confined nuclear fusion reactor, a new born star, the giver of 
life. The elements we need to exist were created in stars during their long lifetimes and sometimes their 
fiery ends. It is difficult to accept the concept that once the nuclei in our bodies were in such an extreme 
environment. 

Stars are the most important objects in the Universe. Even though they make up less than 3% of the matter 
they make up 100% of the visible matter. They are a major apparatus for studying the universe and its 
evolution. This includes the remnants left from their deaths and the discarded material from stellar winds 
during their life. The main reason why we began studying the stars must have been their phenomenal beauty 
in the night sky. On a clear and dark night it is easy to become captivated if you do something very simple, 
look up. Mankind has, with little doubt, always been drawn to the sky since we started thinking and shaping 
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Figure 1.1: Discovery rate of SNe up to 2003. Before 1990 supernovae are binned into decade totals, after 
1990 binned into annual totals. 

our environment. We initially thought that the sky was fixed and any changes we observed would directly 
affect our fate and destiny. At least this allowed enterprising astronomers to secure funding by always 
providing a good spin on astronomical events for the rulers of their countries. 

Once permanent written records of the sky were made we had to abandon the idea that the sky was unchang- 
ing forever. It was realised that the sky did evolve: stars move relative to each other, they are not as fixed 
as we thought. Planets would be the most obvious stars wandering across the sky. However things changed 
forever when the once in a lifetime spectacular event occurred, a supernova. Luckily this happened when two 
of the most famous astronomers in history, Tycho Brahe and Johannes Kepler, were alive. The supernovae 
of 1572 and 1604 captivated them both and they began to provide theories on what they were. Some of the 
ideas were wilder than others. The most interesting idea was that the supernova of 1604 was seen at the same 
time as a conjunction between a number of planets and it was thought that this cosmic ballet gave birth to 
the new star. 

We now know that supernovae are in fact the death of stars. The events have been observed in history as 
far back as the beginning of the second millennium according to current records. If we assume that Homo 
sapiens have been around for 300,000 years then there have been about 6,000 supernovae in our galaxy 
during our presence on the planet. Assuming we can only see nearby supernovae then maybe 2,000 of these 
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would have been noticeable in the sky . We have now observed over 2000 supernovae in nearby galaxies by 
diligent observing schemes. Figure 1.1 shows the number of supernovae seen over the time we have been 
searching for them. The future of supernova studies looks bright. The growth in the area has primarily been 
driven by cosmologists using type la supernovae to determine the geometry of our Universe but the study of 
supernovae themselves is now starting to be seen as the new important output of the observations. The study 
of supernovae will filter down to the study of stellar evolution, limiting the possible outcomes forcing us to 
refine our models and approach a comprehension of stars and their fiery deaths. 

1.2 Stars and their Evolution 

"Space is big. You just won't believe how vastly hugely, mind-bogglingly big it is. I mean, you may think 
it's a long way down the road to the chemists, but that's just peanuts to space." 
The Hitch-Hikers Guide to the Galaxy, by Douglas Adams. 

The scale of stars is beyond everyday experience. The temperature, radius, luminosity and density are all far 
beyond the environmental parameters humans can survive. What is most confounding is the large range of 
variation between stars themselves. Four factors determine how stars evolve, initial mass, initial composition, 
mass loss and duplicity. 

Initial mass determines the pressures and temperatures the core of a star achieves and therefore how far 
the nucleosynthesis and evolution progress. Initial composition affects the evolution only slightly moving 
the mass ranges for different behaviour to occur. We take the composition of stars to be the amounts of 
hydrogen, helium and metals, metals being everything that isn't hydrogen or helium. When the Universe 
first formed the only constituents were hydrogen and helium and a small fraction of other light elements and 
isotopes. The first generation of stars produced carbon, oxygen and other elements. This pollution of the 
Universe with fusion waste caused the stars to change. Increasing contamination increases the opacity of 
stellar plasmas and makes stars less compact. The extra carbon, nitrogen and oxygen enable a star to burn 
hydrogen more effectively by catalysing the reaction. 

Mass loss reduces the mass of a star during its evolution. In the most extreme cases this can completely 
remove the hydrogen (from a star) leaving it as a naked helium star. Mass loss can be achieved by stellar 
winds or by a companion if the star is in a binary. The latter can also provide the opportunity for mass gain 
with mass transferred from one to another. We describe duplicity in more detail below. Figures 1.3, 1.4 
and 1.5 contain three Hertzsprung-Russel (HR) diagrams showing examples of the evolution of stars with 
different initial masses. We describe the main different classes of star below. 
'Of these 588 would have been type la, 235 type Ibc and 1,177 would have been type II. 
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Figure 1.2: Flow diagram for single star evolution. Green - paths available for low mass stars, blue - paths 
available for intermediate mass stars and red - paths gyjailable for high mass stars. 



1.2 STARS AND THEIR EVOLUTION 



1.2.1 Low-Mass Stars, 0.08 < M/M Q < 0.8 



These are stars that burn hydrogen very slowly on the main sequence. Once all hydrogen is exhausted and 
the helium core is formed the burning moves into a shell around the core. This causes the star to become a 
red giant. The exact reason why a star becomes a red giant is unknown but it is thought to be due to nuclear 
burning in a shell. As the star becomes a red giant a deep convective envelope forms that mixes material 
from the core, mainly helium, to the surface. This is the event of first dredge-up. The core is supported 
by electron degeneracy pressure. Shell burning continues until the envelope is ejected by a wind during the 
giant phase and a helium white dwarf is formed. The evolution timescale of these stars is longer than that 
of the age of the universe so very few of these stars will have evolved all the way to helium white dwarfs 
unless they were in a binary system. In figure 1.2 we show the evolution paths of these stars, low-mass star 
evolution paths are shown in green. Figure 1.3 shows the HR diagram for a O.5M star. It slowly evolves 
on the main sequence and moves up the giant branch. Mass loss grows with ascension of the giant branch 
and slowly removes the envelope. Eventually removal of the envelope causes the giant to collapse to a white 
dwarf. This is when the star moves across the HR diagram to the white dwarf cooling track along which the 
star slowly loses heat to come into thermal equilibrium with the cosmic microwave background radiation. 



1.2.2 Intermediate Mass Stars, 0.8 < M/M Q < 8 



These stars evolve as before but do ignite helium in the core. In the lowest mass stars this happens in a 
degenerate core and a helium flash occurs. In degenerate material the pressure only depends on density. 
Therefore as the temperature increases the reaction runs away. Energy output increases rapidly and the 
region continues to expand until the degeneracy is removed and the burning continues quiescently. For more 
massive stars helium ignites non-degenerately. The helium burning results in the star dropping down the red 
giant branch. The presence of a central energy source removes the cause for the star to be a giant. 

The burning continues until the production of a carbon oxygen (CO) core. The CO core grows as the helium 
burning shell burns outward. With the end of central helium burning the star re-ascends the giant branch on 
the asymptotic giant branch, AGB. Second dredge-up, similar to first dredge-up, now occurs and again mixes 
nuclear processed material to the surface. In the most massive stars the convective envelope can penetrate 
deeply into the helium core to assist the helium burning shell to catch up with the hydrogen burning shell. 
The star's luminosity increases as the hydrogen burning shell is now at a higher temperature. Eventually the 
two burning shells are in close proximity. The arrangement is unstable and thermal pulses occur. The cycle 
is that the helium burning shell ignites, quickly burning the helium in the thin region up to the hydrogen 
burning shell. This pushes the hydrogen burning shell into lower temperature regions and the luminosity of 
the star decreases. The helium burning only occurs for a few decades. Afterwards the hydrogen burning 
shell reignites and burns for around 10,000 years producing more helium until the helium burning shell 
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Figure 1.3: HR diagram for a low mass star, Mjnitiai = 0.5M G . The surface temperature changes little during 
the main sequence and giant branch however once the envelope is removed and the star becomes a white 
dwarf to move on to the cooling track. 
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Figure 1.4: HR diagram for an intermediate mass star, an AGB star, Minimal = 5M . The star moves up the 
main sequence along the Hertzspung gap before moving up the giant branch until helium ignites at which 
point it moves around the blue loop before reascending the giant branch and undergoing thermal pulses at 
the tip of the branch. 
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reignites and the cycle begins again. The detail of these thermal pulses is extremely complex and uncertain 
and leads to third dredge-up that mixes nuclear burnt material to the surface at each thermal pulse. Because 
of this AGB stars are important sources of nucleosynthesis. Eventually they loose their envelopes forming 
CO white dwarfs. 

Towards the maximum mass of this range carbon ignition is possible in the core. This occurs in a shell 
around the centre in the less degenerate region. Slowly the degeneracy of the inner regions is lifted and the 
burning continues inwards. Once this brief phase of burning is complete the star becomes a Super-AGB star. 
It is identical to an AGB star but has an oxygen neon (ONe) core. Then there is an overlap with high mass 
stars. In this work if the stars lose their envelope and form ONe white dwarfs they are intermediate-mass 
stars. 

In figure 1.4 the 5M star is an AGB star. We see it evolves on the main sequence. At the end of core 
hydrogen burning it moves across the Hertzspung gap as its surface cools and then ascends the giant branch. 
Once helium ignites the star drops down and goes round the blue loop before reascending the giant branch. 
Eventually at the top of the giant branch the star undergoes thermal pulses and its final position varies slightly. 
The core grows very slowly during thermal pulses and mass loss by the star's wind will remove the hydrogen 
envelope to leave a white dwarf as the star's remnant. 



1.2.3 High-Mass Stars, M > 8M ( 



These stars are the progenitors of supernovae and do not form white dwarfs. Their remnants are neutron stars 
or, in the most extreme cases, black holes. The bottom end of this range merges with intermediate mass stars 
and they are Super-AGB stars that go supernova. ONe cores can undergo core-collapse by electron capture 
on to 24 Mg and other isotopes present. With increasing mass second dredge-up does not occur before the end 
of nuclear burning reactions and burning continues all the way to the formation of iron group elements at the 
centre of the star. Since iron is the most stable element no further reactions are exothermic and thus the core 
collapses forming a neutron star or black hole. The release of the energy from the core collapse initiates the 
supernovae. 

The evolution to this end point has many different paths. The essential factors are initial mass and mass loss. 
These are linked since to a first approximation more massive stars undergo more severe mass loss. The mass 
loss can be so severe that all hydrogen is removed from the star so that it becomes a naked helium star. If the 
initial mass is greater than around 20M Q the helium star is in fact a Wolf-Rayet star. These have strong mass 
loss from their surface and have extremely high surface temperatures, over an order of magnitude greater 
than the sun. There are a sequence of Wolf-Rayet stars, WN, WC and WO, where the second letter indicates 
the dominant element in the spectrum, nitrogen, carbon or oxygen. If the mass loss is only mild the final 
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Surface Tornperaturi* i log(~RK) 

Figure 1.5: HR diagram for three massive stars, Mj n i t i a i = 10, 20 and 4OM . Evolution is similar to the 
intermediate mass stars although only the 1OM goes round the blue loop and they do not undergo second 
dredge-up to experience the large change in luminosity as they reascend the giant branch. The 4OM loses 
most of its mass on the giant branch and moves back across the diagram to the blue region to become a 
Wolf-Rayet star. 
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outcome of the star is a red giant. 

In figure 1.5 we show the evolution of 10, 20 and 4OM stars. We can see the 1OM star has a similar final 
luminosity to the AGB star on the HR diagram. The 2OM star evolves in a similar manner but the final 
motion up the giant branch is smaller. The 4OM star goes up the red giant branch and experiences mass loss 
so severe that it moves back to the blue and becomes a Wolf-Rayet star. 



1.2.4 Stellar Winds and Mass Loss 

Mass loss modifies the evolution of a star by affecting the surface conditions and the total mass of the 
star. While theoretical mass-loss rates do exist (Kudritzki, 2002; Vink, de Koter & Lamers, 2001; Nugis 
& Lamers, 2002) most are empirical due to our limited understanding of the mechanism of stellar winds in 
most cases. Mass-loss rates must be divided into two main categories, those for Wolf-Rayet (WR) stars that 
have lost their hydrogen envelope, and pre-WR or hydrogen-rich evolution. The nature of mass loss changes 
once the hydrogen envelope has been removed and observations show quite different stellar winds. WR stars 
have greater mass loss than a main sequence star of the same luminosity. The mechanism for mass loss in 
OB stars is radiatively driven winds. The WR star mechanism is not known for certain. It could be driven 
by radiation, pulsations, magnetic fields or some other uncertain process (Heger & Langer, 1996; Nugis & 
Lamers, 2002). 

The results of evolution are strongly dependent on the mass-loss prescription. Most observations of the mass 
loss from stars are at solar metallicity and there are few observations at different metallicities. The commonly 
adopted scaling with metallicity is that M(Z) = M(Z ) x (Z/Z ) 5 due to the work of theoretical models 
of Kudritzki, Pauldrach & Puis (1987); Kudritzki et al. (1991). This scaling arises from the assumption that 
stellar winds are line driven and with lower surface opacity at lower metallicity there are weaker winds. 
However while there is agreement that mass loss scales in this form there is a range of suggested values 
for the exponent. It is assumed that this scaling applies at all stages of evolution even though there is some 
evidence that it varies with spectral type. For example Vink, de Koter & Lamers (2001) suggest that for hot 
stars the exponent should be 0.69. This seems likely because the exponent assumes the winds are line driven 
while mass loss from red giants is driven by some quite different unestablished mechanism. At very low 
metallicites where no observations exist we must rely on theoretical mass-loss rates of Kudritzki (2002) or 
extrapolate the empirical rates into this region. 

The most recent rates for WR are from Nugis & Lamers (2000) and these depend on luminosity and surface 
abundance. They find that the mass-loss rate depends strongly on composition as well as luminosity. This 
dependence gives lower mass-loss rates than the earlier and widely used rates of Langer (1989) which link 
the mass-loss rate to the WR star's mass and the equivalent rates that Nugis & Lamers (2000) also derive. 
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However it is common to use the rates of Langer but lower them by a factor of two or three to bring them 
into line with the more recent observations of clumping in WR winds (Hamann & Koesterke, 1998). 

1.2.5 Duplicitous Stars 

Single stars are somewhat unusual. Most stars are in binaries or higher order multiple systems and higher 
mass stars are more likely to have a companion. The presence of a companion gives the opportunity for 
new evolutionary paths. These are mass loss either by the star slowly expanding until the radius of the star 
reaches a point where the gravity of the secondary is greater than the primary and matter is transferred or 
the star could rapidly expand and a common envelope could occur, where the star's core and its companion 
orbiting in the envelope which is removed by some currently unknown process. The evidence that both occur 
is that we see systems which can only be explained if mass transfer has occurred. For example in Algol the 
less massive star is more evolved than the more massive star and this can only be explained if the mass ratio 
has changed at some point in the past. Planetary nebulae contain a large number of close binaries. We also 
see helium stars in tight orbits that can only be explained if a common envelope has occurred. 

Binary systems also provide opportunities for mass gain. If one star is losing mass it is probable that the 
other star is gaining this discarded matter. This can increase the mass of the secondary star to the point where 
it will go SNe when it would not have done so without it. There are other binary effects and the physics is 
quite uncertain and much has still to be studied in this region even though we currently have some basic 
models that fit observations reasonably well. In any study of stars an understanding of binary stars is of 
fundamental importance. 



1.3 Numerical Astronomy 

"Computers make efficient and practical servants, but I do not wish to serve under them." 
Spock, Star Trek, The Ultimate Computer. 

Whenever there is something to understand it is important to make a model. By changing the parameters of a 
model we can see how the output (of our model) changes and thus gain insight into the physical object. The 
modelling of stars started with Lord Kelvin's calculations for the lifetime of the Sun if it was powered by 
gravitational collapse. However theoretical stellar evolution started with Schwarzschild, Henyey and others 
using physical laws to create increasingly sophisticated stellar models. The process invented was similar to 
that used now, the solution of differential equations by numerical rather than analytical processes. 
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The invention of the computer during the Second World War to break codes is probably the most important 
invention that has shaped modern society. Later on the invention of the silicon transistor and affordable home 
computers pioneered by Sir Clive Sinclair accelerated the impact and about twenty years later provided a 
source of research students with great experience with computers. The first to apply computers to help our 
study of stellar evolution were Henyey et al. (1959). While then computer time was limited and at a premium, 
we now discard it, wasting it to show pretty screensavers. There is therefore vast CPU time now available. 
We are now able to run multiple models with different parameters to gain further understanding. 

This dissertation is the result of 70 years of CPU time, that is the equivalent of one computer calculating for 
70 years. This has been achieved by use of multiple computers over the shorter time of 3 years. However 
it is important that even a gargantuan amount of data is not in itself useful. Computers are incredible idiots. 
They do exactly what we tell them to do. 



1.4 Supernovae 



"I've never seen a supernova blow up but if it's anything like my old Chevey Nova 

it'll light up the night sky!" 
Fry, Futurama, Roswell that Ends Well. 



Supernovae (SNe) have most probably been observed since the dawn of astronomy. The earliest known 
surviving record of a supernova was made in the eleventh century by ancient Chinese astronomers (see 
Stephenson & Green 2002 for the details of these early observations) but the modern study of SNe began 
with Baade and Zwicky in the 1930s when they realised SNe are more luminous and rarer than the more 
common novae. The high luminosities and broad spectral lines led them to conclude that SNe are extreme 
explosions, the result of a dying star collapsing to a neutron star. This basic model still holds. Sixty-five years 
of work has vastly increased our understanding but many questions about SNe remain unanswered. 

A core-collapse SN occurs once a star has a core (usually of iron group elements) that cannot be supported 
by any further nuclear fusion reactions or electron degeneracy pressure. This leads to the core collapsing to 
a neutron star or black hole. In the process it releases a large flux of neutrinos that interact weakly with the 
surrounding stellar envelope depositing a large amount of energy (10 44 J) that can drive the star to explode. In 
massive cores supported by electron degeneracy pressure core-collapse can be initiated by electron capture 
on to nuclei. 

Important questions are, 'Which stars give rise to which SNe and leave which type of compact remnant?' 
The answers provide a test of stellar evolution theory, give information for galactic chemical evolution and 



31 



1.4 SUPERNOVAE 



predict SNe birth rates and remnant populations. To solve these questions we must use our understanding 
of stellar evolution to model stars up to the point at which they become SNe. Three main factors affect 
the answer, the initial mass and metallicity of the star and mass-loss (or mass-gain in some binary systems) 
during its lifetime. Heger et al. (2003), using the models described in Woosley, Heger & Weaver (2002), de 
Donder & Vanbeveren (2003) and Podsiadlowski, Joss & Hsu (1992) are example of studies of the effect 
of binaries. While their studies are informative each only investigates a single mass-loss scheme, does not 
investigate the lowest mass for a star to become a SN nor describes the SN progenitors and does not fully 
explore the factor of metallicity. These factors limit comparison with observations. We expand upon this 
work and add details on the progenitors of SNe. These details are important for surveys of progenitors now 
underway by various groups, e.g. Smartt et al (2003) and Van Dyk et al. (2003). 

There are two main types of SNe, those without hydrogen in their spectra, type I, and those with hydrogen 
in their spectra, type II. Type I progenitors have lost their hydrogen envelopes and are white dwarfs or WR 
stars, while type II progenitors have retained their hydrogen envelopes. This gives us a basic method to 
discriminate between them. 



1.4.1 Type II 



A star gives rise to a type II SNe if there is any hydrogen in its envelope. If enough of the original hydrogen 
is retained then it gives rise to a type IIP supernova. The light curve plateau phase is powered by a moving 
hydrogen ionisation front. Type IIL SNe only retain a small fraction of their hydrogen and the light curve 
decays linearly powered by the decay of 56 Ni. There are a few other subtypes, lib are similar to IIL but are 
thought to have less hydrogen and begin looking like a type-II SN but eventually look like a type-lb SN. Type 
Iln have narrow lines in their spectra indicating slow expansion speeds. They are thought to be the result of 
interaction with a dense circumstellar environment. Type Ilpec are where all the other SNe go if they don't 
fit into one of the above categories. 



1.4.2 Type I 



Type I SNe are divided into types la, lb and Ic. SNe la are extremely bright events thought to be the 
explosive carbon burning of a degenerate white dwarf that has reached the Chandrasekhar mass (Mch) by 
accretion; this is not a core collapse event. However types lb and Ic are core-collapse SNe in helium stars 
that have lost their hydrogen envelopes in a wind or by binary interaction (Ensman & Woosley, 1988). They 
are discriminated by the presence or absence of helium lines in the spectrum. However the differentiation 
between type lb or Ic from limited observations can sometimes be difficult so it is common to label them 
together as Ibc. The largest progenitors are thought to give no display because the core is so massive that 
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even with a large explosion energy nothing escapes the forming black hole. The exception would be if a jet 
driven SNe occurred that makes a black hole and produces an observable display (MacFadyen, Woosley & 
Heger, 2001). 



1.4.3 Hypernovae and Gamma-Ray Bursts 

Some supernovae such as SN2003dh are extremely energetic with around ten times the energy of normal 
supernovae and are called hypernovae (HNe). Gamma-ray bursts (GRBs) are objects where a tremendous 
burst of gamma-rays is observed over a few seconds outputting a similar amount of energy. It is thought 
these events are supernovae and all similar events of the same process but in slightly different progeni- 
tors. This is confirmed by observations of GRBs associated with SNe such as SN1998bw/GRB980425 and 
SN2003dh/GRB030329. It is also possible that hypernovae are the results of failed GRBs according to the 
models of MacFadyen & Woosley (1999). The relationship between SNe, HNe, GRBs and their relative rates 
is quite uncertain and much remains to be confirmed. The current situation is discussed by Podsiadlowski et 
al. (2004). 



1.5 Dead Stars 

"I always wanna feel this way, 
Just like a Phoenix from the Flames." 
Robbie Williams, Phoenix from the Flames. 

A star dies when it can no longer produce energy from nuclear fusion reactions. Low and intermediate- 
mass stars lose their envelopes leaving a cooling core while massive stars undergo supernovae. These events 
produce different stellar coipses. These dead stars are normally less luminous than stars that are producing 
energy by nuclear fusion. However the remnants are not necessarily dead objects but can be reborn to 
shine once more. If they have a binary companion that overfills its Roche Lobe it is possible for these 
remnants to accrete matter. The more compact the object the more extreme the conditions achieved in the 
accretion flow. As matter spirals towards the surface it can achieve high velocities and become heated to a 
high temperature so radiation is emitted. Neutrons stars and black holes when they accrete material from 
a companion are extreme X-ray emitting objects and many such binaries are observed. These spectacular 
objects were discovered by the first X-ray telescopes and provide direct evidence for the existence of these 
remnants. 
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We show in figure 1.6 the number of objects known in each of the remnant classes and their masses. We 
can see there is some overlap between these objects and the uncertainty is greatest in black hole masses. It 
is not surprising that most known remnants are white dwarfs since these should be more common as they 
come from low mass stars and they are more populous because the initial mass function makes low-mass 
stars more probable than higher-mass stars. 

1.5.1 White Dwarfs 

These are essentially the cores of low and intermediate mass stars once no further nuclear reactions are 
possible. They are supported by electron degeneracy pressure and calculations show that the maximum mass 
of a star supported by degeneracy pressure is the Chandrasekar mass, Mch, which is around 1.4M Q although 
the exact value varies very slightly with composition. They usually have a thin layer of hydrogen on their 
surface, the remains of the envelope, and cool slowly over time. They are roughly the size of the Earth, 
however the more massive the star the smaller the radius so that as the mass approaches Mch the size shrinks 
to around 10km. 

White dwarfs come in three flavours depending on how far the nuclear reactions in the progenitor's core 
have progressed. They can be made of helium (a He WD), a mixture of carbon and oxygen (a COWD) or a 
mixture of oxygen and neon with some sodium and magnesium (an ONeWD). HeWDs are formed by very 
low mass stars or more probably in binary stars and are not of interest in this work. COWDs are the remnant 
cores of stars that lose their envelopes on the AGB. The Super- AGB stars develop an ONe core and are the 
source of ONeWD. 

1.5.2 Neutron Stars 

"neutronium, but..." 
Scott, Star Trek, A Piece of the Action. 

Neutrons are unstable and free neutrons decay to a proton and electron with a half-life of around 8 minutes. 
They are however stable in a nucleus or neutron star which is effectively one giant nucleus held together 
and made stable by the gravitational binding energy. They are around 10km in radius and there are multiple 
models for their internal structure. It is thought that the innards might consist of a quark gluon mix rather 
than actual neutrons. 

Whatever their structure they are supported by neutron degeneracy pressure that is similar to electron de- 
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Figure 1.6: Measured masses for compact remnants. Object number is the indentifier used for each object 
ranked in order of lowest to highest mass. 



generacy pressure but since a neutron's mass is roughly 2000 times that of an electron a greater mass can 
be supported. The theoretical maximum mass for a neutron star is 3M e , above which the star collapses to 
a black hole. Neutron stars were initially believed to be hypothetical but the detection of pulsars and their 
interpretation as rapidly rotating neutron stars is the primary evidence of their existence. Their association 
with supernovae is evidenced by the presence of a pulsar in the Crab nebula. The nebula is the other type 
of supernovae remnant commonly referred to, that is the expanding ejecta from the explosion rather than the 
compact object. 

Neutron stars are formed by one of two mechanisms. In the lowest-mass stars to go supernova nuclear 
burning progresses as far as carbon burning to form an ONe core. If this core has a mass equal to Mch it 
collapses and once the central density reaches around 10 9 8 gcm~ 3 electrons are captured on to the nuclei 
present, most notably 24 Mg. The removal of electrons removes the electron degeneracy pressure and the 
collapse is hastened until a neutron star is formed. The release of neutrinos in the formation of the neutron 
star is thought to lead to the transfer of energy to the envelope causing it to be ejected. 

In more massive stars burning continues until a core of iron group elements is formed. Iron, being the 
most stable element, is the end of possible nuclear burning. Endothermic reactions occur including the 
photodisintegration of iron and others that lead to the core collapsing yet again to a neutron star and again 
it is thought the release of neutrinos leads to the ejection of the envelope and the supernovae. After the 
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supernova and the envelope has been ejected the neutron star is left alone. 



1.5.3 Black Holes 



Leela: 



Oh my god we're heading straight for a black hole! 
"Talk about a mood killer." 
Futurama, A Flight to Remember. 



Fry: 



Determination of the remnant formed at the heart of a SN is a black art. The physics is extreme and poorly 
understood. However there are many prescriptions which rely on the conjecture that a more massive core 
leads to a more massive remnant. As the remnants become more massive neutron degeneracy pressure cannot 
support the collapsing core and a black hole is formed. The physics of black holes is deceptively simple, their 
escape velocity is the speed of light therefore not even light can escape from these objects, hence the name 
black holes. However they aren't really that black. Single black holes produce (a weak) Hawking radiation. 
Black holes in a binary however may accrete material from their companion. This leads to an accretion disc 
that gets very hot and therefore very bright at X-ray wavelengths. Many such binaries are observed and there 
are low-mass and high-mass X-ray binaries named after the mass of the donor star. These objects can only 
be formed in supernovae or by accretion on to neutron stars that then subsequently collapse. 

The minimum mass for a black hole is classically 3M Q however it seems that they can form from remnants 
of only 2M G . The mass distribution of black holes is also interesting. We can only find black holes if they 
are accreting and hence luminous sources. So an important question is: what is the initial black hole mass 
function? Can we compare this to the observed population and work out whether most black holes are at the 
masses they were born with or have they accreted a substantial amount of material? 



The rest of this dissertation is split into six chapters. In chapter 2 we discuss the stellar evolution code that 
is used to produce our models. We also detail the opacity tables that were constructed to enable the study 
that has been undertaken. In chapter 3 contains a detailed look at the structure of SNe progenitors, trying 
to estimate from their properties which SNe type they might give rise to. We use observations of SNe and 



1.6 Dissertation Outline 



Data: 



"Do you think this is a wise course of action, sir?' 
"We're about to find out, Data." 
Star Trek: Nemesis. 



Picard: 
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their progenitors to help us reach our conclusions. Also in this chapter we examine the lowest mass of star 
to undergo supernova. This is a complicated area and the results greatly depend upon the physics used in our 
models. Chapter 4 contains a study into the various mass-loss prescriptions that are used and the effect they 
have on the progenitors of supernovae. This is split into three main sections: first we investigate the effect 
of including metallicity scaling for Wolf-Rayet stars; second we compare different mass-loss prescriptions 
at high metallicities (0.001 < Z < 0.05) and finally we compare different mass-loss prescriptions at low 
metallicities (10~ 8 < Z < 0.001). From this we choose a preferred scheme that best fits observations. 

Chapter 5 extends our study to include binary stars to search for low-luminosity progenitors and progenitors 
that have a different structure that we do not see from single stars. Chapter 6 combines our preferred results 
for all single star progenitors and discusses the implications of the observed progenitors to date for our 
models. We then combine the results from single stars with those from binaries so we can draw some 
conclusions about the relative fraction of single to binary stars and how closely we can match the observations 
that exist. We also discuss the distribution of progenitors on the HR diagrams indicating where single stars sit 
and discuss the relative position of binary stars. Finally we detail the remnant mass population and compare 
this to the limited observations that exist. In the conclusion, chapter 7, we summarise our results and list our 
final judgement on the progenitors of the SN types. 
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Chapter 2 



Nuclear Reactions & Opacities 



"I can't change the laws of physics; I've got to have 30 minutes." 
Scotty, Star Trek, The Naked Time. 



2.1 The Cambridge Stars Code 

The Cambridge stellar evolution code was originally written by Peter Eggleton during his time as a PhD 
student; the first paper on this code was Eggleton (1971). Since this time the code has been used, revised and 
adjusted by many people, in fact there are now at least four different stellar evolution codes that are based 
on the original. Because of this it is necessary to describe the features of the code used in this work. The 
full list of papers that detail the code are Eggleton (1971), Eggleton (1972), Eggleton, Faulkner & Flannery 
(1973), Eggleton (1973), Han, Podsiadlowski & Eggleton (1994), Pols et al. (1995), Pols et al. (1997) and 
Pols et al. (1998). 

The code is lightweight being less than 2000 lines of Fortran 77. It is simple to operate and to modify. This 
makes adding extra physics or routines to study wide ranging problems in stellar astrophysics easy and rapid. 
The code also runs remarkably fast on modern computers. Some time comparisons are detailed in appendix 
??. Stars can be evolved from main sequence to carbon burning in around 5 minutes. Therefore it is possible 
to produce large grids of stars with different masses, metallicities, mass-loss schemes, mixing schemes and 
varying other physics to study the effect on a star's evolution. 

The essential features that make it unique among stellar evolution codes have survived its long life span. They 
are the use of a self-adaptive non-Lagrangian mesh, the treatment of both convective and semiconvective 
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mixing as a diffusion process and the simultaneous and implicit solution of both the stellar structure equations 
and the diffusion equations for chemical composition. 

The fact the Eggleton code solves for the structure, composition and the adaptive mesh simultaneously, an 
extremely elegant way to calculate the evolution of stars, is responsible for making the system so simple and 
adaptable. The calculations are fully consistent and remove possible errors when rapid composition changes 
might effect the structure. It also means little or no intervention is required to evolve a star. 

In this chapter, after discussing the code, we shall detail two of the modifications we have made. First is the 
addition of more nuclear reactions and nuclei to the reaction network and second, to update and improve the 
opacity tables. Other additions, such as inclusion of various mass-loss rates and a new routine for following 
binary evolution, will be presented in later chapters. 



2.1.1 Assumptions 



As with any physical model there are always assumptions that have to be made to reduce the problem to a 
level where it might be realistically tackled. It is best to think that over time we refine our physical models 
and come closer to the truth of nature but that the truth is something we can never reach. 

Nearly all stellar evolution packages share the approximation of stars as perfect spheres. This reduces the 
problem to one dimension with the structure varying in radius only. However real stars are not spherical and 
effects such as rotation and duplicity distort the structure. It is unlikely to be important however because, in 
the majority of cases, current models do fit observations. But there are some interesting cases that require the 
need to invoke aspherity (Pinsonneault, 1997). However including extra dimensions increases the complexity 
of the problem and its solution. For example, for a ID star only 199 meshpoints are required for most simple 
stages of evolution. This increases with more dimensions and thus requires more computing power. The 
only known fully 3D code at the moment is Djehuty (Eggleton et al, 2003) although it is much slower than 
simpler ID codes. 

The other two primary assumptions are that the star is always in local thermodynamic equilibrium and in hy- 
drostatic equilibrium. The first assumption does break down in the atmosphere of Wolf-Rayet stars. However 
this has little effect on the interior hydrostatic core. The second assumption is broken when dynamical ef- 
fects occur, for example during rapid mass loss (M > O.OlM yr _1 ), burning in degenerate regions and rapid 
burning stages such as oxygen burning when nuclear reactions occur on a dynamical timescale. However 
these occur rarely in stellar evolution and we shall note when these features could effect the results. 
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2.1.2 Equations 

The basic equations of stellar structure are taught in undergraduate courses and found in most astrophysics 
text books. Solving the coupled equations is not easy. Analytically assumptions must be made such as 
uniform composition, constant opacity or a simple equation of state. Calculation of realistic stellar structures 
is a demanding process. Let us start by listing the equations of stellar structure. They are derived in terms of 
radius which is not practical for calculating the evolution of stars due to contrasting length scales for the core 
and envelope. It is therefore common practice to recast the equations with mass as the variable over which 
the equations are solved. The equations are conservation of mass, 



dm 7 dr 1 

— = \%pr L or -t- = - a j, (2.1) 

dr dm Anpr 1 



hydrostatic equilibrium, 

dp Gpm dp Gm 

-T = — or T" = -^' (2 - 2) 

dr r z dm 47lr* 

energy production, 



dh r i dL 

— L =<\npr 1 e or — = e, (2.3) 
dr dm 



radiative transport, 



dr 3 xp L r dT 3 k L r 

— = — — — or — = - — (2 4) 

dr 16a T 3 4nr 2 dm 16a T 3 16nV 



or convective transport, 

dr y-irdp dr y-ir Gm mn 

— = - -. or — = -- 7 . (2.5) 

dr 7 p dr dm y p Anr 1 

Also required is an equation of state. A simple example is 

P = ^BL + aT 4 . (2.6) 
M 
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The variable r is the radius of a sphere centred on the centre of the star, m is the mass contained within 
that sphere, p is the density at r, p is the pressure, L r is the luminosity through the sphere, e is the energy 
produced per unit mass of material, T is the temperature, a is the Stefan-Boltzman constant and y is the 
adiabatic exponent. 

The driving forces behind the evolution of stars are the nuclear reactions that alter the composition. Therefore 
we need equations that evolve the composition change due to nuclear reactions and mixing by convection or 
other processes. 



Where X, is the mass fraction of the element in question, £ is the diffusion constant or rate of mixing, X,7?</,/ 
is the rate at which Xj is destroyed by nuclear reactions and R c j is the rate at which X,- is created by nuclear 
reactions. 

These coupled differential equations have no (known) analytic solution. To solve them we turn to numerical 
methods. Originally this meant sitting down with a slide rule and working as quickly as humanly possible 
becoming an original computer. Then, thanks to the invention of the electronic brains to crack codes in 
World War 2, the process became faster and reliable. 

The process requires the differential equations to be recast as difference equations and the problem is solved 
over a finite mesh. The resolution of this mesh needs to be high enough that errors do not creep in when a 
section is under-resolved but low enough so as not to compromise calculation speed. 

The difference equations are mass conservation 




(2.7) 




(2.8) 



hydrostatic equilibrium 




(2.9) 



energy transport, 




(2.10) 
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where V = dl °|p and is calculated from radiative or convective theory in their respective regions. Then we 
have the equation for energy production, 



r^dS Gm ,„ „ ^dm\dm\ 

-^=\(^ + e v -T- + C P T^(V-V a )-)-\ , (2.11) 



we are also required to solve for the mesh, 
/ dm\ 

m ^- mk= {jk) k+ l (2 " 12) 



and the composition, 



dmi 



-{Xi,k+\ — Xi,k) -3- , +(Xi,k— ^i,fc-l 



k 



dt 



k+l 



(2.13) 



Where I = K[V r - VJ 2 M 2 / W , V r = ( jgl) = -r^^^and V„ = (fg) Values in square 



brackets are set to zero if their value is negative, t nuc is the current nuclear timescale, M the star's current 
mass and K is a suitable constant chosen as discussed in Eggleton (1972) to achieve mixing on a timescale 
similar to that from Mixing-length theory. Each of these equations is evaluated at points from the centre of 
the star to the surface. The location of these discrete points forms a mesh, the distribution of which depends 
on the mesh spacing function. 



2.1.3 The Mesh and its Spacing Function 

The mesh is dependent on the value of a variable, Q, which is a function of physical variables calculated for 
the star. Those used are the pressure, mass, temperature and radius. Mesh points are then placed with equal 
spacing in Q so that ^ is constant over the mesh at each timestep. The present version of the non-Lagrangian 
mesh-spacing is a function of local pressure P(r), temperature T(r), Lagrangian mass m(r) and radius r. A 
stellar model is well represented by 199 meshpoints although in certain cases where the evolution is complex 
such as thermally pulsing AGB stars more are required. If an increase in resolution is required meshpoints 
cannot be added by hand. It requires either raising the total number of meshpoints in the mesh or more 
commonly refining the mesh-spacing function to move meshpoints into the troublesome region. The most 
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common way to do this is to use the pressure at the point of interest because it is a monotonically decreasing 
variable from centre to the surface. 

To make the process automatic the pressure at the burning shells is found for each model and input into the 
<2 equation, 



Where a, are constants that can be altered to increase or decrease the influence of each term on the mesh- 
spacing. The different subscripts to the pressures represent the pressures at that burning shell which are 
found by looking for a large composition gradient in the element being burnt. These pressure terms produce 
great changes in the value of Q when the two pressures are near equal thus many mesh points are put at 
these points. There are other ways of forming these equations to push mesh points between burning shells. 
However in this case with multiple burning shells we need to resolve these regions more accurately than the 
space between them. The counter example would be in thermally pulsing AGB stars where it is important to 
resolve the intershell convection zones accurately. 

2.1.4 Boundary conditions 

To solve any differential, or difference, equation boundary conditions are required. The conditions are 
applied at the surface and centre of the star. At the surface, where k = 1, the boundary conditions are 
for the conservation of mass 



Q = oti InP + a 2 ln(P + Ph) + a 3 (ln(P + P He ) + ln(P + P c ))) 



/(m/M)J \ (- + A- « 6 ln ( -I— ) 

V a 4 / \c 3 J \T + a 7 J 



(2.14) 



— — Mwind + ^binary, 



(2.15) 



emission from a black body 



L r = KR 2 o B T e 4 Q 



(2.16) 



the point at which photons can escape to infinity 




(2.17) 
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We take the surface of the star to be the photosphere when the optical depth is |. The boundary conditions 
take account of radiation pressure to prevent super Eddington luminosities. In taking the outer mesh point to 
be the photosphere we assume in equation 2.17 that the gravity, g, is due to the mass of the star and therefore 
the photosphere has negligible mass. This is a valid assumption for most stars, any discrepancy will be 
greatest for red giants. 

The simplest central boundary conditions would be L = r = m = treating the centre as an infinitesimal 
point. Currently the central point is modelled as a spherical region. Therefore the condition at the centre 
where k = 199 are for mass conservation 



dm 

m = -—, (2.18) 
dk 



the radius, 



4 » 

m=-7lr i p, (2.19) 



and the central luminosity 

L = (e - TS)m. (2.20) 

With the equations and central boundary conditions it is possible with a suitable number of variables to solve 
the problem of stellar structure and evolution. 

2.1.5 Solution 

In the above equations we require 14 variables to model the stellar structure. The first eleven, already within 
the code, were lnT, lnm, g, lnr, L, X( l U), X( 4 He), X( 12 C), X( l6 0), X( 20 Ne) and In/, a quantity related to 
electron degeneracy. To increase the composition network extra variables were added X( 14 N), X( 24 Mg) and 
X( 28 Si). Nitro gen was included to resolve problems with AGB and low-metallicity stars in the CNO cycle, 
while magnesium and silicon were included to follow later stages of nuclear burning. 

The solution method is based on the Henyey technique (Henyey et al., 1959). An initial model is input, 
usually from a previous calculation. If the model does not satisfy the equations the variables are relaxed 
to a solution, hence the alternate name for this technique, the relaxation method. The relaxation method 
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iteratively calculates how to improve the values of the variables to achieve a better solution to the difference 
equations in a fully implicit manner. The iterations repeat until a certain level of accuracy in the solution is 
reached. Otherwise if a maximum number of iterations is reached the method is aborted since either it has 
insufficient resolution or the initial guess solution was too poor. 

The most important features of making sure the code converges to a solution are timestep and mesh control; 
effectively the resolution in time and space. Unfortunately owing to the nonlinearity of the partial differential 
equations this is not yet an exact science and something of a black art. To determine the next timestep the 
code compares the changes in structure from the last timestep to a constant that determines the maximum 
change required between models. If a large time step is allowed then the code can evolve a star with few 
timesteps although accuracy is reduced. Using a smaller timestep can also produce problems but the main 
drawback is that more models are required for the same period of evolution. 

By experimenting with time step control the number of breakdowns of the code can be reduced. This is 
necessary when computing large numbers of models. The main points of breakdown are ignition of nuclear 
burning in a highly degenerate region, burning shells close to the surface of a star, extreme mass-loss rates 
(greater than 10 M yr ) and short timescale nuclear burning reactions (i.e. Neon and Oxygen burning) 
requiring a low timestep (less than 10~ 4 yr). 

2.1.6 Convection and Mixing 

The treatment of convection was first described by Eggleton (1972). The mixing is modelled as a diffusion 
process so we can correctly follow mixing due to convection and equilibrium semiconvection. The composi- 
tion changes are calculated at the same time as the structure equations making this code unique. Convection 
also effects heat transport and must be considered carefully. To decide whether a meshpoint is stable to con- 
vection we must calculate two temperature gradients, the adiabatic gradient V a = ^ diogp ^) anc ^ tne ra diative 

gradient V r = ( d }^ T ) . These are related to equations 2.4 and 2.5 respectively. 

V BP / radiative 

We then use the classical Schwarzschild criterion for instability, normally expressed as 

V r >V a , (2.21) 

for convective instability to occur. Otherwise the region is stable and heat transport is by radiation and 
no mixing occurs (Z = 0). If the system is convective then mixing is treated as a diffusion process with 
I = (K(V V - V a ) 2 M 2 )/T nuc where K is a large constant calculated from mixing length theory, M the stellar 
mass and T nuc the current nuclear timescale. Where V r — V a is small but positive a semiconvective region 
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occurs. Heat transport is as calculated from mixing length theory with the length over which mixing occurs 
a free parameter in terms of the pressure scale height H p . The mixing length is set to / = aH p as in the 
formalism of Bohm-Vitense (1958). We take a = 2.0 which gives a good fit to the Sun (Pols et al., 1997). 

However there is some evidence that extra mixing beyond that given by standard ID convection is required 
in some cases as discussed by Pols et al. (1997) and Pinsonneault (1997). Extra mixing can be included in 
the form of convective overshooting. We do not parameterise the overshoot by a fixed fraction of the pressure 
scale height. Instead we alter the classical Schwarzschild criterion above to 

V r > V a - 8, (2.22) 



where 



2.5 + 20£ + 16£ 2 



(2.23) 



and £ is the ratio of radiation to gas pressure and <5 ov is the overshooting parameter that we vary to include 
convective overshooting or to remove it. Further details of the inclusion of this scheme in the Eggleton 
code can be found in Schroder, Pols & Eggleton (1997). By comparing models to observations of eclipsing 
binaries with an evolved component they found a best fit for 5 ov to be 0. 12. This leads to overshooting lengths 
l ov between 0.25 and 0.32Hp in stars of mass range 2.5 to 6.5M . For more massive stars the overshooting 
length remains around 0.3Hp. Another way of quantifying the effect of overshooting is to measure the mass 
of the convective core of a ZAMS model. For the 6.5M Q star the mass over which mixing occurs grows from 
1.6 to 2.1M . For a 10M o star the growth is from 3.2 to 4.1M e and for a 1OOM star from 83M to 87M . 

In our calculations we use 8 0V = 0. 12 to include convective overshooting and set 8 ov = to remove over- 
shooting. The amount of overshooting affects the core mass during evolution. This affects the nature of the 
remnants and also the mass of stars that undergo second-dredge up. There is evidence from asteroseismology 
observations that overshooting does occur in massive stars (Aerts et al., 2003). However the extra mixing 
owing to convective overshooting in our models can be considered to mimic any mixing process such as 
rotation or gravity wave mixing, not just convective overshooting. 



2.1.7 Equation of State 

The equation of state (EoS) is described by Eggleton, Faulkner & Flannery (1973) and Pols et al. (1995). 
The EoS is necessary to describe the relationship between pressure, density and temperature over the vast 
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range of regimes that stars encounter. There are a number of ways of including the equation of state in a 
stellar evolution code. One is to interpolate in tabulated data. Within this code the EoS is approximated by 
ingenious formulae to an accuracy of < 0.1%. They provide the details of electron density, pressure, internal 
energy, pressure ionisation and Coulomb interactions. 



2.1.8 GridCheck the Babysitter 

"I think there is a world market for maybe five computers." 
Thomas Watson (1874-1956), Chairman of IBM, 1943. 



The code is not perfect and occasionally numerical errors occur, usually when the mesh does not resolve a 
region that is evolving rapidly or an error occurs in the solution routine. This results in one of two problems, 
the code can get stuck in an infinite loop and outputs no more data or the variables go infinite and the code 
continues to run spewing out blank results indefinitely 1 . To get round this problem a babysitting program 
written in Perl was devised to check on the program while it runs, stopping it when it encounters a fault. 
Figure 2. 1 shows the program structure. 

This was necessary to investigate a large parameter space when a great number of stars must be evolved in 
parallel around the clock. It would be difficult for any human to continuously perform this task. To exploit 
the SunGrid at the IoA this program had to be developed along with the logical structure to evolve hundreds 
of stars automatically. It is now simple to change some of the parameters and see the effect on a group of 
models. While it is not close to the number of models that are possible in simplified population synthesis 
algorithms, e.g. Hurley, Pols & Tout (2000), it has the advantage of producing full stellar evolution models 
rather than approximations. This work would not have been possible without this program. 



2.2 The Composition Variables and Nuclear Reaction Network 

"Change is the essential process of all existence." 
Spock, Star Trek, Let That Be Your Last Battlefield. 



This results in a disk full of data that's as useful as a pair of big fat pants. 
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Figure 2.1: Flow diagrams for the grid check program. The blue boxes are decision statements, red arrows 
are the paths when the outcome is false, the green true. 
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2.2.1 Producing the Network 

Stars derive their energy from nuclear fusion reactions. The first and longest phase, is that of hydrogen 
burning to helium. Next helium is transmuted into carbon and oxygen, the remaining phases are then carbon 
burning, neon burning, oxygen burning and finally silicon burning. Once an iron core is formed the star is 
doomed to supernova. 

Not all these stages occur in all stars. Below certain mass limits the conditions required for reactions cannot 
be reached. Also mass loss can truncate the evolution of a star and expose the core so that it becomes a white 
dwarf before further reactions can occur. 

To study these reactions in a stellar evolution code the abundances of elements must be variables. In the 
Eggleton code we are limited to the number of elements we can study without reducing the performance of 
the code because we solve for the composition simultaneously with the structure and each requires a second 
order equation. In the original code the elements taken to be variables were 1 H, 4 He, 12 C, 16 and 20 Ne. 
Nitrogen was dealt with as the missing mass fraction required for the sum of all elements to be one. The 
reaction scheme for these elements is listed in the first column of table 2.1. Most of these reactions are 
correct, however there are approximations due to the small number of variables. Noticeably the nitrogen 
alpha addition and the later burning stages are contrived. 

To make a more accurate and realistic network we added three extra composition variables. This does slow 
down the code but not excessively thanks to modern computing power. The full list of elements followed 
is now 'H, 4 He, 12 C, 14 N, 16 0, 20 Ne, 24 Mg and 28 Si with the remaining fraction taken to be iron. The new 
reaction scheme is shown in the second column of table 2.1. We retain the nitrogen alpha addition since while 
this reaction should lead to 22 Ne only small amounts of this element are formed. Although it is an important 
neutron source in AGB stars it is likely to have a small effect on the structure. The later stages are now 
slightly improved and can be followed all the way to oxygen burning, although during these burning stages 
the code experiences numerical problems as the time step drops to less than 10~ 4 yr and the assumption of 
hydrostatic equilibrium is likely to be broken with dynamical effects altering the progress of burning. In the 
new reaction scheme there are rates stored in the code that are not currently used. This is because they occur 
after oxygen burning which we currently can only see the ignition of. 

2.2.2 Testing the results, a 25M© star 

Figures 2.3, 2.4 and 2.2 show the results of the evolution of a 25M Q star. These models were evolved, 
the mass-loss rates taken from de Jager, Nieuwenhuijzen & van der Hucht (1988). The figures display the 
internal composition structure, evolution over the HR diagram and timestep versus model number. Both 
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Table 2.1: Nuclear Reaction Network. Those schemes in blue had the rates included in the code but not 
used in the reaction network. The first five reactions, comprising the pp-chain, are assumed to be in transient 



Old network. 


New Network. 


l H (p, 6+v) 2 H (p,y) J He 


'H (p, 6+v) 2 H (p,y) J He 


3 He ( 3 He,2p) 4 He 


3 He ( 3 He,2p) 4 He 


2 He ( 4 He,7) 7 Be 


2 He ( 4 He,7) 7 Be 


7 Be(e~,v) 7 Li (p,a) 4 He 


7 Be (e~,v) 7 Li(p,a) 4 He 


7 Be (p,/) 8 B QS+v) 8 Be* (a) 4 He 


7 Be (p,7) 8 B (0+v) 8 Be* (a) 4 He 


12 C (p,j3+v) 13 C (p,7) 14 N 


12 C (p,j3+v) 13 C (p,7) 14 N 


14 N (p,j8+v) 15 N (p,/) 16 


14 N (p,/3+v) 15 N (p,7) 16 


14 N (p,/3+v) 15 N (p,a) 12 C 


14 N (p,/3+v) 15 N (p,a) 12 C 


16 (p,/3+v) 17 (p,a) 14 N 


16 (p,j3+v) 17 (p,a) 14 N 


4 He (a) 8 Be* (a,j) 12 C 


4 He (a) 8 Be* (a,y) 12 C 


12 C (a,7) 16 


12 C (a,7) 16 


14 N (a,7) 18 F {\a,y) 20 Ne 


14 N (a,7) 18 F (i«,7) 20 Ne 


16 (a,7) 20 Ne 


16 (a,7) 20 Ne 


20 Ne (a,7) 24 Mg 


20 Ne (a,7) 24 Mg 




24 Mg (a,7) 28 Si 


12 C ( 12 C,«7) 2() Ne 


12 C ( 12 C,a7) 20 Ne 


12 C ( 12 C,7) 24 Mg 


12 C ( 12 C,7) 24 Mg 


12 C ( 16 0,«7) 24 Mg 


12 C ( 16 0,«7) 24 Mg 


16 ( 16 0,a7) 28 Si (7,a) 24 Mg 


16 ( 16 0,a7) 28 Si 


20 Ne (7a) 16 


20 Ne (7,a) 16 


24 Mg (7,a) 20 Ne 


24 Mg (7,a) 20 Ne 




28 Si (7,a) 24 Mg 
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Figure 2.2: Comparison of a 25M stars.Top is the HR diagram and bottom is the evolution of the time step 
with model number. 
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Table 2.2: Comparison of models for a 25M Q star. 





INO IVldSS LOSS 


IVldSS LOSS 


IVldSS LOSS 


ivduscner ei di. 




5 OV = 0.12 


5 OV = 0.12 


S ov = 0.00 


(2002) 


Initial mass 


25.00 


25.00 


25.00 


25.136 


Final mass 


25.00 


16.27 


18.40 


13.079 


He core mass 


10.78 


10.06 


8.47 


8.317 


CO core mass 


7.32 


6.96 


5.33 


6.498 


ONe core mass 


2.14 


1.77 


1.66 


2.443 


Remnant mass 


4.00 


3.17 


1.49 


1.961 



diagrams are similar and show the hydrogen, helium, carbon and neon burning shells. In none of the stars 
do we progress to oxygen burning. Although the star with no mass loss does ignite oxygen but the model 
fails soon afterwards. Figure 2.2 shows that for this model the timestep has dropped to extremely low values 
and our assumption of hydrostatic equilibrium is likely to break down. The code becomes unstable at low 
timesteps as terms in the energy equation are large with only a small difference between them. 

In table 2.2 we compare our models to those of Rauscher et al. (2002), they use the mass-loss rates of 
Nieuwenhuijzen & de Jager (1990) which are similar to those we have used but do lead to greater mass loss 
in models. Using these rates in our models produces similar core masses however the final mass drops to 
12.2M e . They do not use convective overshooting and full details of their treatment of convection are given 
in Woosley, Heger & Weaver (2002). The closest agreement is with our third model without convective 
overshooting which leads to smaller values for the helium and CO cores. Our oxygen neon core mass is 
smaller because we are not able to evolve the model through to the last burning stages. The remnant masses 
are also of the same magnitude although our higher values are due to convective overshooting producing 
larger cores. 

The differences that are there are due to the way we define the core masses, the nuclear rates adopted and 
the details convection. The rates of hydrogen and helium burning are well known in comparison to the later 
burning reactions so it is not suiprising that the final core masses are different. The treatment of convection 
especially in the late stages of evolution is quite uncertain. During this time the convection will probably 
affect the energy production and vice versa. Our main conclusion is that while there are is some uncertainty 
in the physics of our models and we are not able to evolve through the final stages of evolution we do gain 
similar results to other codes for the earlier burning stages. 
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2.3 Opacity 

The main source of a star's energy is the nuclear fusion reactions occurring either in its core or in thin burning 
shells around the core. This energy is transported from the production site to the surface by radiative transfer 
or convection. The first transports energy in the form of photons, while the second is a cyclic macroscopic 
mass motion that carries the energy in bulk. 

In regions stable to convection, radiative transfer leads to the equation of stellar structure, 
dT 3krP L r 



dr 16ar 3 4Kr 2 



(2.24) 



where L r is the luminosity at radius r, p the density, a the Stefan-Boltzmann constant, T the temperature 
and Kr is the Rosseland mean opacity. 

Opacity K is a measure of the degree to which matter absorbs photons. There are four main sources of 
opacity in stars 



• Bound-bound transitions are the transitions of an electron in an atom, ion or molecule between energy 
levels which are accompanied by either the absoiption or emission of a photon. Only a photon with 
the correct wavelength can cause a given transition, the process is wavelength dependent. 

• Bound-free transitions, or photoionisation, occur when an incoming photon has enough energy to 
ionise an atom or ion and free an electron. The reverse process is the capture of an electron by an atom 
or ion. It will not occur until photons above a threshold energy are available and falls off as v~ 3 where 
V is the frequency. 

• Free-free transitions are scattering processes which occur when an electron and photon interact near 
an atom or ion. The process is also known as bremsstrahlung. Again it is proportional to v~ 3 . 

• Electron scattering is wavelength independent at low temperatures where it is Thomson scattering. 
The electron, with its low cross-section, is a small target and so only dominates at high temperatures 
when most atoms are ionised. At very high temperatures relativistic effects are important and Compton 
scattering dominates. 



From this list it is possible to see that calculating the opacity in a stellar model is a difficult process and 
depends on the composition, temperature and density of the material. Deep in the star local thermodynamic 
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equilibrium (LTE) is achieved and an average for radiative transfer over all wavelengths requires the Rosse- 
land mean opacity (Rosseland, 1924) expressed as 



where K v is the opacity at frequency v, T is the temperature and B v is the flux per unit area into unit solid 
angle per unit frequency in LTE. 

The most comprehensive opacities available today are from the OPAL (Iglesias & Rogers, 1996) or OP 
(Cunto et al., 1993) groups who have made detailed models of the above processes. They provide tables of 
the Rosseland mean opacity variation with temperature, density and composition with the metal abundance 
usually scaled to solar compositions. Iglesias & Rogers (1993) took a step forward for the OPAL project 
team by providing tables that include mixtures enhanced in carbon and oxygen (C and O) relative to the 
base solar composition. We have incorporated their full range of composition tables in the Eggleton (1971) 
evolution code. In the last implementation of this code (Pols et al., 1995) only 10 of the 265 tables available 
were used. Now we have taken the opportunity to refine our models so as to accurately follow the changes 
of opacity which occur in the later stages of evolution. We can expect the new tables to affect many types of 
stars. While they may be small for the evolution of main-sequence stars and white dwarfs we expect larger 
differences to be found in AGB and Wolf-Rayet stars. 

Asymptotic Giant Branch (AGB) stars undergo third dredge-up which mixes helium burning products to the 
surface and forms carbon stars. Using the enhanced mixture tables we shall be able to model the thermal 
pulses and envelope evolution more accurately. Wolf-Rayet (WR) stars are massive and have lost their 
hydrogen envelopes exposing the helium cores. As time progresses helium burning products are slowly 
exposed at the surface and in some cases the stars are eventually mostly composed of carbon and oxygen. 

We shall present the method of opacity table construction and detail its implementation in the Eggleton 
stellar evolution code. We shall discuss the effects on main-sequence stars, red giants and white dwarfs. We 
shall present three tests of these tables. The first is the effect of including extra carbon and oxygen on the 
structure and evolution of a low-mass population-Ill star, the second is a 5M thermally pulsing AGB star 
and the third a Wolf-Rayet star of 4OM with mass loss. 




(2.25) 
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2.4 Opacity Table Construction 

We start with the OPAL tables (Iglesias & Rogers, 1996) as the framework around which to construct the 
full tables in an approach similar to that employed by Pols et al. (1995). However we choose to use the 
variable M = p/T 6 3 rather than the density. This has the advantage that our tables can be a third smaller in 
this dimension. This partly compensates for increased memory requirement of the full tables. 

The OPAL tables extend in log 10 (r/K) from 3.75 to 8.70 in steps of 0.05 to 0.2 and in log 10 (^/gcm" 3 MK" 3 )| 
from —7 to +1 in steps of 0.5. We choose the range of the final tables to cover from 3 to 10 in log 10 (r/K) 
with increments of 0.05, and from —8 to +7 in log 10 (^/gcm~ 3 MK~ 3 ) with increments of 0.5. This gives 
us tables of 141 by 31 elements. 

To account for the effect of composition, the OPAL tables include typically 265 T and M grids for each 
metallicity. These tables are split into 5 groups with different hydrogen mass fractions X of 0.0, 0.03, 0.1, 
0.35 and 0.7. The tables in each group have different compositions in helium, carbon and oxygen. 

The lowest temperature in the OPAL tables is log 10 (r/K) = 3.75. Below this we use the tables of Alexander 
& Ferguson (1994), which extend from 3 to 4.1 in log 10 (r/K). Where the tables overlap they match well 
(Iglesias & Rogers, 1996). It should be noted that the tables provided by Alexander & Ferguson (1994) do 
not include enhanced carbon and oxygen mixtures. Therefore at temperatures when log 10 (r/K) < 4.0 our 
tables do not follow these mixtures. This is only important in the surface of AGB stars. We plan in future to 
deal with low temperature enhanced mixtures by including the work of Marigo (2002) before application to 
the structure of AGB star envelopes. 

To complete each table we fill the region from 8.70 to 10 in log 10 (r/K), where electron scattering dominates, 
according to Buchler & Yueh (1976). We then combine with a full table of the effective opacity owing to 
electron conduction for which we use the fits of Iben (1975) to the tables of Hubbard & Lamp (1969) 
and Canuto (1970). We combine the radiative opacity with the effective conductive opacity by reciprocal 
addition, 

1 1 1 

= + . (2.26) 

Keff Krad ^cond 

Before the complete tables are implemented in the stellar evolution code we perform one last process. We 
rescale the tables over the C and O plane to a more regular grid. In the OPAL grid the spacing used requires 
interpolation within irregular polygons, so to make the interpolation simpler and faster, we include extra 
tables in a similar but fully rectangular grid system. This means we have 61 tables for each hydrogen 
abundance or 305 tables for each metallicity. 
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To construct tables of different metallicity care must be taken since there are no corresponding OPAL tables 
for some of the low temperature tables and vice-versa. In these cases we construct two tables with the 
nearest two OPAL tables both including the same low temperature tables and interpolate. A similar method 
is applied for the OPAL tables. 



2.5 Implementation into the code 

The main drawback with implementation of the new tables is the increased memory requirement. Whereas 
before there were only 10 opacity tables there are now 305. The new tables are each a third of the size of the 
old tables so the overall increase is a factor of ten. This corresponds to a memory requirement increase from 
17 Mb to 170Mb but the typical specification of today's desktop computers comfortably accommodates this. 
One of the reasons why this memory requirement is so large is that we also store the spline coefficients for 
each of the tables. When we use the tables within the code we interpolate in the ffl and T plane with bicubic 
splines. This ensures the smoothness of our tables. We then use linear interpolation in the 3 composition 
dimensions to derive the final opacity value at each point. 

To test the effect of these tables we use three different interpolation schemes, 

• Method A uses 5 tables with X = 0, 0.03, 0.1, 0.35 and 0.7 and Y = 1 -X -Z. Where X is the 
hydrogen mass fraction, Y the helium mass fraction and Z the initial, solar mixture, metallicity. We 
add a further two tables, both with X = Y = 0, the first with Xq = 1 — Z and Xo = 0. The other with 
Xq = and Xo = I — Z. Where Xc and Xo are the mass fractions for the enhanced amount of carbon 
and oxygen above that included in the metallicity mass fraction so the total carbon mass fraction is 
Xc + Z(C) and for oxygen Xo + Z(0). This gives a total of 7 tables allowing simple interpolation in 
composition similar to the old method. 

• Method B is similar to method A but we include two extra tables for each hydrogen mass fraction. 
Again both have Y = 0. One has X c = 1 - X - Z and X = and the other X c = and X = 1 - X - Z. 
This leads to a total of 15 tables in all. 

• Method C incorporates all 305 tables in the opacity calculations. We have 61 tables for each hydrogen 
abundance. When interpolating in composition we use three variables X, Xc/( 1 — X — Z) and Xo/( 1 — 

X — Z). The hydrogen abundances are as for method A. We take values of 0, 0.01 , 0.03, 0. 1 , 0.2, 0.4, 0.6| 
and 1 for the carbon and oxygen grid planes. On this grid we perform a three-dimensional linear in- 
terpolation. 
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2.6 Testing & Results 

The new tables give almost identical results to the old for main-sequence stars and red giants. The timescale 
of core helium burning is altered. From method A to method C it increases by a factor of 0.1% for a 5M Q 
star. This is reassuring. Current stellar evolution models describe well the evolution of the main sequence 
and the red giant branch. However we do find that the opacity tables enhance numerical stability in the later 
stages of evolution because they resolve changes with composition in greater detail. This leads to a smoother 
variation in opacity. 

A similar result is found for white dwarfs. Between methods A and C for a carbon and oxygen white dwarf 
the radius increases by 0.2% and cooling times are affected by about 1% with the initial cooling faster with 
method C than for method A until later when the cooling times are longer with method C. While conditions 
in the atmosphere do pass through regions of the greatest difference between the two methods (the dotted 
line in figure 2.7) there is not enough mass at these temperatures and densities to make a significant global 
difference. 



2.6.1 Polluting a 0. 5M Q Population-Ill Star 

"Like Granny said, if you want a box hurled into the sun, you got to do it yourself.'" 
Hermes, Futurama, The Farnsworth Paradox. 

For our first detailed test we take a zero metallicity 0.5 M star and pollute it uniformly with carbon and 
oxygen. Such a star could have been formed just after the death of the first stars if the products of these 
stars were mainly carbon and oxygen. Or such stars could be made by the hypothetical industry of stellar 

Table 2.3: Variation of Evolutionary Time Scale with increasing carbon and oxygen mass fraction. Here 
?a,?b and tc are the evolution times to the helium flash for methods A, B and C and tz is the evolution time 
for models with an equivalent mass of scaled solar metallicity. 
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Table 2.4: Variation in physical parameters relative to method A with increasing carbon mass fraction. 
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engineering by an advanced civilizations. We choose such a low-mass star to limit the effect on the results of 
burning via the CNO cycle. Our initial model has the following properties, X = 0.7, Y = 0.3, R = 0.45 R , 
?eff = 4780 K and L = 0.96 L e . The tables provided by OPAL do not include one with X = 0.75, we could 
extrapolate to obtain the opacity but for this test we prefer to remain within the tables. From the initial 
model we perform three sub-tests. In the first we pollute the star at a constant C/O ratio of 1 , then evolve 
each model and compare the time from zero-age main sequence to the helium flash for methods A, B and C. 
We also compare a sequence of models for stars of the same mass but with a scaled solar metallicity of the 
same mass fraction as the carbon and oxygen. These results are presented in table 2.3. We see that with the 
pollution fraction increase, the timescale of hydrogen burning increases because the energy produced at the 
centre takes longer to reach the surface. However when comparing the difference between methods B & C 
we see only a comparatively small increase in timescale. When we compare with the scaled solar models we 
see a greater difference. This is due to the much larger opacity of heavier elements such as iron. 

With the second set of models we look at the radius, surface temperature and luminosity of the star on the 
main sequence as we pollute it with either only carbon or only oxygen (tables 2.4 and 2.5). We use method 
A as the base to which we compare methods B and C. A similar trend is seen in the results in tables 2.4 and 
2.5. We therefore see it is important to include the affect of enhanced C and O mixtures before hydrogen 
exhaustion. Similar mixtures are encountered in the envelopes of carbon stars of very low metallicity. 

Finally we consider the effect of varying the C/O ratio of the pollution (table 2.6) at a constant pollutant 
mass fraction of 0.001. Table 2.6 demonstrates the need to use method C to follow accurately the variation 
in opacity as the C/O ratio changes. This is important in AGB evolution during the thermal pulses when the 
C/O ratio changes gradually as carbon becomes dominant. The difference between methods B and C in this 
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Table 2.5: Variation in physical parameters relative to method A with increasing oxygen mass fraction. 
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Table 2.6: Variation in physical parameters relative to method A with varying the C/O ratio. 
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Table 2.7: Variation of timescales for a 5M star. 
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case comes about because of non-linear structure in the C and O opacities. It can only be resolved with the 
full set of tables. 

2.6.2 Thermal Pulses of an AGB Star 

Our thermally pulsing AGB model is of a 5M star. To follow a number of pulses we include convective 
overshooting during the pre- AGB evolution. This leaves the star with a larger core which reduces the rate at 
which the thermal pulses grow in strength. This means we do not find dredge-up in our models but we can 
evolve through the pulses. We obtain the average pulse period as the arithmetic mean of the times between 
the first ten helium shell flashes. Results are recorded in table 2.7. 

First we find that the different methods affect the age of the star when it starts to undergo thermal pulses. 
Second we find that methods A and B have similar interpulse periods while method C's is 3.6% less. We 
attribute this to the change in opacity due to varying He/C/O ratios in the inter-shell region which alters the 
time evolution of the pulse. 

2.6.3 Evolution of a Wolf-Rayet Star 

We have used the mass-loss rates of de Jager, Nieuwenhuijzen & van der Hucht (1988) and Nugis & Lamers 
(2000) to manufacture a Wolf-Rayet star of solar metallicity and an initial mass of 40M e . Some details are 
described by Dray & Tout (2003). Our tables are expected to make a difference because Wolf-Rayet stars 
become naked helium stars and helium burning products are mixed throughout and exposed at the surface of 
the star. We only use methods A and C and not method B because the carbon and oxygen mass fractions do 
not rise above those included in the base metallicity until after the hydrogen envelope has been removed and 
so our tables only begin to make a difference after this. 

In figure 2.5 we present the variation of opacity through the star, just after core carbon ignition, with mesh 
point. The solid line is with method A. We then take this structure and use method C to calculate the 
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opacities to give the dashed line. This model is not in hydrostatic equilibrium and the dot-dashed line shows 
the relaxed new structure in hydrostatic equilibrium from method C. Comparing the physical parameters of 
these calculations we find the main change to be an increase in radius of 3.0%. The surface temperature has 
decreased by 1.5%. 

We then made two complete runs with similar initial conditions by methods A and C. In figure 2.6 we display 
the difference between Wolf-Rayet evolution as the hydrogen envelope is removed. The greatest changes are 
in the radius and surface temperature. The star is in the WN phase between 14 to 12M and then enters the 
WC phase. In the WN phase there are only small differences but upon entering the WC phase the differences 
build up so that the temperature is lower by around 2% and the radius about 1% larger. The luminosity also 
tends to be larger by around 0.5 %. 

Figure 2.7 is a plot of the difference in opacity between methods A and C over the p and T plane when there 
is no hydrogen present. The plot is of the mean difference in the opacity calculated by methods A and C at 
all possible carbon and oxygen abundances divided by the opacity when Y = 1 — Z. From the plot we see 
that the regions of greatest difference are small. The lines show where the Wolf-Rayet model calculated by 
method C lies at certain points of its evolution. The solid line is for the star at 12.3M and the dashes is the 
same star at a later time when its mass has fallen to 6.1M . From figure 2.6 we see that the largest changes 
begin around the 12M point. This is when the model has moved into the region where the differences 
between A and C are greatest. 

We find the opacity tables alter the timescales for the late burning stages. Method C increases the helium 
burning timescale by 0.1% while the carbon-burning time scale increases by 2.5% relative to method A. 
These changes come about because the opacity varies with composition and effects the temperature structure 
and burning rates. While there is no large difference in a single variable the overall effect on structure and 
evolution is more significant. 



2.7 Conclusions 



Our major conclusion is that the changes induced by properly including opacities for varying C and O mix- 
tures are small. However they are a good thing to include in stellar evolution codes because they help 
numerical stability by removing sharp changes in the variation of opacity with composition that are encoun- 
tered when we use fewer tables and add another level of detail to models without much loss of computational 
speed. 

The main computational cost is the extra memory to store the tables. Method C requires ten times the 
memory (about 200Mb) compared to methods A and B. However there is little computational speed cost: 
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Figure 2.5: The variation of opacity throughout the Wolf-Rayet star just after core carbon ignition. The solid 
line is calculation from model A, the dashed line from model C and the dash-dotted lines from model C with 
hydrostatic equilibrium. The dotted line is the mass interior to the mesh point. Mesh point is at the surface 
and 1000 is at the centre of the star. 
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Figure 2.6: The percentage difference in model parameters between methods A and C versus mass of the 
Wolf-Rayet star. ^ 
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relative to method A, methods B and C require about 4% more time for the evolution of a Wolf-Rayet star 
even though we must evaluate four times as many opacities for each M and T value than with method A. 
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Figure 2.7: The value plotted by the contours is J s (kq — Ka) dXc dXo / J s Kq dXc dX . Here Ka and Kq are the 
opacities calculated by methods A and C. The opacities are all calculated when X = and over the surface 
(S) defined by (Xq +Xo) = to 1. The vertical artifacts are from the wide spacing of the tables in p. We 
show two Wolf-Rayet models on this plot of the same star at different times, solid line 12M and dotted 
6.1M Q . This second model is the same as shown in figure 2.5. We see that the largest differences in the 
models indeed occur as the star drifts into the region where the opacity difference is greatest. Also plotted 
is the dotted line for a 0.3M e pure C/O white dwarf model. While the line passes through regions of the 
diagram with greater values the white dwarf is only slightly affected because very little mass is in the regions 
with greatest difference. 
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Chapter 3 



The Evolution of Single-Star Supernovae 
Progenitors 

"Life is pleasant. Death is peaceful. It's the transition that's troublesome." 

Isaac Asimov. 



3.1 Introduction 

The three main factors that determine whether a single star will go supernova, and the type of supernova it 
will become, are initial mass, composition and mass loss. These factors are linked, generally a more massive 
star experiences greater mass loss and higher metallicity also enhances mass loss. These factors affect the 
composition of the ejecta, the main detail to consider in deciding the resultant SN type. The outcome of 
evolution depends on a competition between mass loss and the nuclear reactions at the centre of the star. If 
the core collapse occurs before much mass is lost a type II SN occurs. However if mass loss is severe all 
hydrogen is removed and a type I SN occurs. In the most extreme cases mass loss can affect the nuclear 
evolution of the core and lead to smaller core masses. 

Single stars lose mass in stellar winds. For red supergiant stars and Wolf-Rayet stars the mass loss mechanism 
is uncertain. There are a number of proposed mechanisms including magnetic fields, opacity effects and 
pulsations (Heger & Langer, 1996). For Wolf-Rayet stars the most promising explanation currently seems to 
be optically thick winds (Nugis & Lamers, 2002). 

Mass loss in stellar models is achieved by altering the mass boundary condition. The rate of mass loss at 
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a particular point in a stars evolution is taken from analytical expressions derived from theory e.g. Vink, 
de Koter & Lamers (2001) and Kudritzki (2002) or empirically determined from observation e.g. de Jager, 
Nieuwenhuijzen & van der Hucht (1988), Nieuwenhuijzen & de Jager (1990) and Nugis & Lamers (2000). 
The current problem is that there are a number of different prescriptions with different groups using their 
favoured rates. Also there is evidence that mass-loss rates may be lower than observed because clumping in 
the winds leads to overestimates of the mass ejected (Hamann & Koesterke, 1998). The classic case of this 
is that the WR mass-loss rates have been lowered from initially determined values by about one third due to 
this effect. There is now growing evidence that clumping is also present in the mass loss of main sequence 
OB stars. This means that there is a tendency for mass-loss prescriptions to overestimate the mass loss from 
single stars. Binary stars are likely to be a source of extra or enhanced mass-loss in stars, rapid rotation could 
be another possibility. 

In this chapter we describe the main types of single star progenitor predicted by theory, with close consider- 
ation of the minimum mass for a SN to occur. We then look in detail at progenitor models and discuss their 
evolution at three metallicities, solar, LMC and SMC. We discuss which details of these models we use to 
determine the type of SN a progenitor gives rise to. 

3.2 Low mass progenitors 

"Balderdash, I'll be the judge of who's cool using the coolometer! Good Lord! 
I'm getting a reading of over 40 MegaFonzies!" 
Professor Farnsworth, Bender should not be allowed on TV, Futurama. 

Stars with zero-age mass less than 2OM give rise to red supergiant SN progenitors. These progenitor stars 
have cool surface temperatures and radii of many hundreds of solar radii. They should be the most common 
single-star progenitors observed. We can split them into two groups distinguished by the core-collapse 
mechanism, either iron core-collapse by photodisintegration or ONeMg core-collapse by electron capture. 

3.2.1 Normal SN Progenitors 

As the normal case we refer to the most common type of SN progenitor a red supergiant where the nuclear 
fusion reactions have progressed from hydrogen burning all the way to silicon burning ending in the for- 
mation of an iron group elements core. While the internal nuclear evolution has occurred the envelope has 
expanded to a star with a radius hundreds or thousands of times that of the sun. With this expansion the 
surface temperature of the envelope drops to relatively cool temperatures of less than 4000K. 
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The pre-SN structure of the core has multiple shells with an iron core at the centre surrounded by the earlier 
nuclear fusion reactions burning in shells with the hydrogen burning shell outermost. After the formation of 
the iron core further nuclear reactions are endothermic so the core collapses by the photodisintergration of 
iron. The SN event is complex and includes much extreme physics (Arnett, 1994; Mueller E & Janka, 1994; 
Mezzacappa & Bruenn, 2000; Warren, 2002). Collapse leads to the formation of a neutron star or black 
hole and releases a vast amount of energy in neutrinos that leads to the ejection of the mass surrounding the 
compact remnant. This gives rise to a type II SN. If little of the hydrogen envelope has been lost, a type IIP 
SN. These are SN that have a long plateau phase where the SN luminosity remains constant. This is because 
the photosphere which is located at the point of hydrogen reionisation remains at a constant radius. This is 
due to the ejecta expanding outward while the ionisation front moves inwards in mass. 

Observations of the progenitor to SN2003gd have found it to be a red supergiant (Smartt el al., 2004). 
This and other data on the progenitors of type IIP SN from Smartt et al (2003) confirm that these normal 
progenitors do exist and maybe quite common. 



3.2.2 Super- AGB Progenitors 

There are two paths of evolution to become an AGB star. For lower masses the helium burning shell catches 
up with the hydrogen burning shell, while for more massive AGB stars (above 3M ) the convective envelope 
penetrates down to the CO core in second dredge up. This puts the hydrogen and helium burning shells 
in close proximity, an unstable arrangement that leads to thermal pulses rather than steady shell burning. 
Furthermore because the hydrogen burning shell is now at a higher temperature near the helium burning 
shell hydrogen burning occurs as a much greater rate and the luminosity of the star increases by at least a 
factor of 10. 

Most AGB stars have CO cores and do not ignite carbon. These go on to form CO White Dwarfs as the 
envelope is removed by a stellar wind during the thermally pulsing phase. If the envelope was not removed 
these stars would eventually form a degenerate CO core approaching Mch and carbon would ignite degener- 
ately in the centre of the star leading to a thermonuclear SNe. The resultant SN would have observational 
characteristics of a type la SN with a type II mask. They are sometimes referred to as a type 1 ^ as suggested 
byIben&Renzini(1983). 

However stars above about 7M Q experience carbon burning as detailed in Garcia-Berro, Ritossa & Iben 
(1997), Iben, Ritossa & Garcia-Berro (1997) and Ritossa, Garcia-Berro & Iben (1999). The lowest mass 
stars in this range ignite carbon in a shell around the more degenerate central regions that have been cooled 
by neutrino losses. The burning shell then slowly burns inwards, raising the degeneracy of the inner shell 
before burning takes place. In more massive stars carbon ignites in the core and then burns outward in brief 
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shell flashes rather than in a steady burning shell. This is due to the mild degeneracy of the CO core when the 
burning begins. Carbon burning leads to the formation of an ONe core. With the new core the star becomes 
a Super-AGB star. The carbon burning continues outwards until it reaches the base of the helium burning 
shell at which point it extinguishes. However the hydrogen and helium shells still continue burning outwards 
forming a CO mantle above the ONe core. It is likely that mass loss will lead to envelope loss and formation 
of an ONeWD. 

In the cases with the most massive cores it is possible that the envelope is not lost before the core reaches 
Mch- This can occur in the most massive stars when the core is close to Mch after second dredge-up. At this 
point the core collapses, however now there is no carbon to ignite and prevent collapse. Densities can be 
reached of around log (p/gcm~ 3 ) =9.6 — 9.8 where electron capture on to 24 Mg can occur. This accelerates 
the core-collapse and formation of a neutron star and a type II SN. The only difference from other SN is that 
the core-collapse is via electron capture rather than iron photodisintigration. 

We find because of our new nuclear reaction rates and opacity tables the mass ranges for different behaviour 
have shifted relative to previous studies. It is possible to summarise the behaviour as follows: 

• Mioos < 7M , Mos < 5Mq, stars undergo second dredge-up and thermal pulses with a central CO 
core as a thermally pulsing AGB star. These lose their envelope and leave CO white dwarfs before 
their cores reach Mch- 

• Mioos ~ 8M , 6M q < Mos < 7M Q , during or after second dredge-up carbon ignition occurs in a shell 
because the core is degenerate and has a temperature inversion caused by neutrino losses. The star 
then undergoes thermal pulses with an ONe core. Since the core mass after dredge-up is less than Mch 
these stars can lose their envelopes and form ONeWDs. These are Super-AGB stars. 

• Mioos ~ 9M ra , Mos ~ 7.5M©, carbon ignites before second dredge-up (in a shell if the centre is 
degenerate). Thus at dredge-up there is a growing ONe core. If this can reach Mch before the envelope 
is lost then the star undergoes a SN. The outcome depends on the nature of the thermal pulses. 

• M no os ~ 1OM , Mos ~ 8M , the CO core is greater than Mch before dredge-up. However shell 
carbon burning (enhanced by a thin neon burning shell in the most massive stars of this type) drives a 
convection zone that reduces the size of the CO core to Mch so dredge-up can occur. This CO material 
is mixed with the envelope and increases the CO abundance at the surface during second dredge-up. 
After dredge-up the star has an ONe core of Mch which progresses to a SN by electron capture on to 
24 Mg. These are Hyper-AGB stars. 

• M 10 os — 1 LM©> ^os > 9M Q , the helium core or CO core masses are too great for dredge-up to occur. 
The limiting mass for the helium core is 3M and for the CO core 1.5M©. Nuclear burning in these 
stars progresses until it cannot support the core and a SN occurs. 
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Where Mos and M no os are the mass ranges with and without convective overshooting included. The lu- 
minosity of SN progenitors determines how many we can expect to find in current observational searches 
for progenitors. If Super-AGB stars do undergo SN then the luminosity of the lowest mass SN progeni- 
tors depends upon the occurrence of second dredge-up in the late stages of evolution (Smartt et al., 2002). 
In the mass range of interest second dredge-up can increase the final luminosity from log (L/L ) < 4.6 to 
log(L/L ) > 5.2. Therefore these stars, if they occur in nature, should provide a large population of luminous 
red supergiant progenitors. 

Most Super-AGB stars probably do not undergo SNe but lose their envelopes and form ONe white dwarfs 
before core-collapse can occur, as in the work of Garcia-Berro, Ritossa & Iben (1997) and Iben, Ritossa & 
Garcia-Berro (1997). Also the much higher luminosity increases mass loss. From our models it can take 
longer than 10 5 years for the burning shells to advance and for conditions in the core to reach those for core- 
collapse. To lose their envelopes during this time a mass-loss rate of less than 10 M Q yr is required. This 
is not unreasonable for a luminous red supergiant. Also if the helium burning shell is unstable and thermal 
pulses do occur, the burning shell's advance is slowed and there is more time for mass loss to occur. 

However there is a class of stars, extreme Super-AGB or Hyper-AGB stars, that have CO cores greater than 
Mch when second dredge-up occurs and a convection zone between the hydrogen and helium burning shell 
forms to reduce the CO core mass to Mch- Figure 3. 1 shows examples of such stars of different masses 
at solar metallicity without convective overshooting. The 9M G model undergoes second dredge-up and the 
core mass grows outwards and collapse occurs at just below Mch- For the 9.5M Q model the CO core mass is 
just below Mch at second dredge-up and a small intershell convection zone forms as the convective envelope 
penetrates downwards. At 10 and 10.5M e the CO cores are now greater than Mch as second dredge-up 
occurs and the intershell convection zone that penetrates into the CO core reducing its mass is now much 
larger and eventually merges with the hydrogen envelope at which point helium burning products such as 
carbon and oxygen are mixed to the surface increasing the observed abundance. 

The carbon burning in these stars can reach quite high temperatures, above 10 9 K. This is high enough for 
some neon burning to occur so the cores acquire a larger fraction of magnesium. These stars do not have to 
wait for the hydrogen and helium shells to burn outward to reach Mch- Neutrino losses cool the core and 
increase the central density to the point where electron capture can occur, we take this to occur at a central 
density of log(p/gcm~ 3 ) = 9.8 (Miyaji et al., 1980; Nomoto, 1987). We define this as the point when the 
SN occurs and we end the evolution of the model. 

Extreme Super-AGB stars are difficult to model. The carbon burning shell flashes cause the timestep re- 
quired to change unpredictably. In the later stages of the most massive extreme Super-AGB stars a neon 
burning shell ignites in the degenerate ONe core before collapse occurs and causes numerical instabilities 
that make further evolution troublesome. To experiment with this important range of masses we make some 
simplified models. We artificially limit the carbon and neon burning rates. This means they do not runaway 
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Figure 3.1: The internal structure of 9, 9.5, 10 and 1O.5M stars with solar metallicity and no convective 
overshooting. The are shown above, top left, top right, bottom left, bottom right. Evolution continues up 
to the point the centres experience conditions for electron capture collapse. The thin black lines indicate 
convective regions, the red line indicates the hydrogen burning shell, the yellow line the helium burning 
shell and magenta line the carbon burning shell. 
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Figure 3.2: The luminosity of solar metallicity SN progenitor models before the central density reaches the 
conditions required for core-collapse. The upper plot is without convective overshooting, the lower plot with 
convective overshooting. All units are in solar units. 
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in degenerate regions where there is no thermostatic control. They proceed to burn but more slowly than 
might be the case in reality. The effect can be see in figure 3.1 where the carbon burning shell moves out in 
a stable convective shell rather than the separate convective shells that occur when the burning proceeds via 
brief flashes of burning. These models provide an upper limit to the maximum mass of extreme Super-AGB 
stars. Neon burning in the ONe core is likely to have a small effect on preventing second dredge-up and at 
the highest masses the burning may progress all the way to the formation of a small iron core and collapse 
might not occur via electron capture. 

In figure 3.2 we plot the luminosity of extreme Super-AGB models with and without convective overshooting 
against the time before the central density is reached when electron capture reactions would initiate. In these 
models we have not included mass loss. It has little effect on the final outcome of these models. The lowest 
mass models are on the limit of Super-AGB to extreme Super-AGB stars and may eject their envelope to 
form an ONe WD. The more massive models increase in luminosity 100 years to 1 year before SN. This is 
because before second dredge-up the stars have more massive cores which leads to a denser core after second 
dredge-up. Therefore less time is required to achieve the conditions for core collapse. 

The luminosity rise from second dredge-up is substantial in all cases. Interestingly in figure 3.2 just before 
the increase there is a slight fall in the lightcurve. This is due to the hydrogen burning moving as the 
convective envelope penetrates into the helium core. The most massive star included in figure 3.2 does not 
experience second dredge-up and is the lowest luminosity progenitor that is possible. The only red supergiant 
progenitor observed, SN2003gd, had a pre-SN luminosity of log(L/L ) = 4.3 ± 0.3. The lowest luminosity 
progenitors are within 1 a of this value. We can be sure that the progenitor did not go through second dredge- 
up because the observations of the progenitor were taken less than one year before the SN. Progenitors with 
a lower luminosity than this model might exist as an outcome of binary evolution or we might consider how 
to prevent second dredge-up. One way to do this would be by making convection less efficient so the process 
of second dredge-up is slower. 

After observing a number of SN progenitors, we shall be able to determine more about the low mass end of 
stars that go SNe. If we make observations of fields containing prospective SN progenitors we shall observe 
that between 9-17% of all SN progenitors will undergo second dredge-up before SN, with overshooting the 
percentage is 10-20% with the IMF of Kroupa, Tout and Gilmore (1993). The exact values depend on the 
range of masses for S-AGB stars to occurs. The range of Super-AGB stars are sensitive to the strength of 
convection used in our models. If we inhibit convection we find second dredge-up does not occur before 
core collapse. 

This is a promising observational test of whether electron capture SNe, and S-AGB progenitors exist. The 
only other possibility is small differences in the nucleosynthesis products in these SN. In normal SN progeni- 
tors the collapsing core is surrounded by shells of silicon, oxygen, neon, carbon and helium while in electron 
capture SN there is only the hydrogen envelope so, important differences may exist, although it appears that 
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a similar amount of 56 Ni is produced in the SN of Super-AGB stars (Swartz, Wheeler and Harkness, 1991). 



3.3 High mass progenitors 



In stars of M > 2OM black holes can be formed in SNe rather than neutron stars and mass loss becomes 
severe enough to effect the evolution of the star and remove a large fraction of hydrogen or in the most 
extreme cases all hydrogen is lost from the star. This loss leaves the core of the star exposed as a naked 
helium star. These massive stars which are loosing mass in a strong wind are called Wolf-Rayet stars. 
These are the last type of SN progenitor. Wolf-Rayet stars are classified by their spectra and there are three 
main classes WN, WC and WO in which different elements are observed. They represent the products of 
different stages of nuclear burning. WN stars still have some hydrogen and nitrogen is the dominant feature 
of the spectrum. WC stars start exposing the products of low-temperature helium burning and carbon begins 
to dominate. Finally WO stars expose the higher temperature helium burning where oxygen becomes the 
dominant element. 

For single stars WR stars are probably formed at M > 3OM at solar metallicity. However there is a com- 
plication in that binary interactions can lead to mass-loss events such as Roche Lobe overflow and common 
envelope evolution so that stars lower than this limit can lose their hydrogen and become Wolf-Rayet stars. 
We must be careful when looking at mass-loss rates not to pick a prescription that reproduces the lowest 
mass WR stars in the sky. These are likely to be in binaries so our single star limit must be higher. 

The stars that retain hydrogen in the progenitor will give rise to type II SN. However they are more likely 
to be type ILL, Iln or lib because they do not have extended hydrogen envelopes or a dense circumstellar 
environment with which the ejecta can interact. If there is enough hydrogen for a type IIP SN it may be very 
different if a black hole is formed at the centre of the collapsing core rather than a neutron star as suggested 
by Zampieri et al. (2003). Although an alternative scenario proposed by Chugai & Utrobin (2000) is that 
these low luminosity IIP SN could also arise from lower mass progenitors with initial mass between 6 and 
1OM . 

Stars that have lost all their hydrogen give rise to type I SNe, either type lb or Ic since la is a thermonuclear 
event. SN lb have helium in their spectra while Ic SN do not. There is some uncertainty in what difference 
in progenitors give rise to the difference between lb or Ic SNe and since most progenitor models still retain 
some helium there is probably a continuum between the two objects but we shall return to this dilemma later. 
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3.4 Mass-loss rates 

Mass loss is the biggest factor in determining the mass at which type I SN occur rather than type II. In 
chapter 4 we investigate different mass-loss prescriptions and their effect on supernovae progenitors over a 
wide range of metallicities and initial stellar masses. For the models below we use our preferred mass-loss 
rates for local metallicities (10~ 3 < Z < 0.05) and the inclusion of extra mixing in the form of convective 
overshooting. We shall examine 3 metallicites, solar, LMC and SMC (Z = 0.02,0.008 and 0.004). 

When deciding which rates to use it is sensible to opt for empirical mass-loss rates if we aim to model nature. 
The mechanism for mass loss in OB stars is radiatively driven winds and theoretical predictions agree with 
observed values (Vink, de Koter & Lamers, 2001). WR theoretical rates do exist but the WR star mass-loss 
mechanism is not known for certain. It could be driven by radiation driving optically thick winds, pulsations 
or opacity (Heger & Langer, 1996; Nugis & Lamers, 2002). Because of this uncertainty it is currently best 
to use empirical rates for WR stars. 

We split our mass-loss rates into pre-WR and WR because the nature of mass loss changes once the hydrogen 
envelope has been removed and observations show quite different rates. WR stars have greater mass loss than 
an OB star of the same luminosity. Our preferred rates are essentially the NL rates from Dray & Tout (2003) 
but with the inclusion of the rates of Vink, de Koter & Lamers (2001) for the OB stars. 

Pre-WR Evolution: We use primarily the mass-loss rates of de Jager, Nieuwenhuijzen & van der Hucht 
(1988), here after JNH. The JNH rates are old and complex but considered to be the most accurate (Crowther, 
2001). The rates are dependent on surface luminosity and temperature and are derived from a large pool of 
observations. We apply the commonly adopted scaling with metallicity of M(Z) = M(Z Q ) x (Z/Z ) 5 as 
in Heger et al. (2003) and Dray (2003). This scaling arises from the assumption that stellar winds are line 
driven and with lower surface opacity at lower metallicity there are weaker winds. However while there 
is agreement that mass loss scales in this form there is a range of suggested values for the exponent. We 
investigate the effect of this scaling in chapter 4. We supplement these rates by with the rates of Vink, de 
Koter & Lamers (2001) for OB stars. These include a theoretically predicted metallicity scaling. 

WR Evolution. WR evolution begins just before the hydrogen envelope is completely removed from the 
star. The definition of when a model enters each stage is arbitrary. We use those denned by Dray (2003). 

• WNL: When X < 0.4 and r eff > 10 4 . 

• WNE: When X < 0.001 and (X c +X ) /Y < 0.03. 

• WC: When X < 0.001 and 0.03 < (X c +X ) /Y < 1.0. 
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• WO: When X < 0.001 and (X c +X Q ) /Y > 1.0. 



When X is the surface hydrogen abundance by mass, Y, Xc and Xo are the surface abundance by number of 
helium, carbon and oxygen. During WN and WC evolution we use the rates of Nugis & Lamers (2000) that 
depend on luminosity and surface abundance. They find that these depend strongly on composition as well as 
luminosity. During WO evolution we use a constant rate as in Dray & Tout (2003). One straightforward test 
of this mass-loss prescription is to compare the number ratios of WR stars to O stars, and the WR subtype 
ratios. Dray & Tout (2003) performed this test in detail and found good agreement with observed ratios. 
The one exception was the WN/WR ratio for which this theoretical value is too low. There are a number 
of plausible reasons for this such as the effect of binaries on mass-loss or misclassification of some WN 
stars because spectral features change with metallicity. There is evidence that WR mass-loss rates should be 
scaled with the initial metallicity (Crowther et al., 2002), we return to this in chapter 4 but here we include 
the scaling with metallicity as for pre-WR evolution of M{Z) = M(Z ) x (Z/Z ) 5 . 



3.5 The Evolution of SN Progenitors 



We can now closely examine the variation between SN progenitors over a range of masses. We look in this 
case at solar, LMC and SMC metallicities and masses from 10 — 2OOM . We do not go to lower masses 
to avoid the electron capture SN region. We have included convective overshooting in these models. The 
compiled HR diagrams for the stars are shown in figures 3.3, 3.4 and 3.5. After evolution on the main 
sequence stars move over the Hertzprung gap before becoming a red supergiant. Higher mass stars then try 
to go up the supergiant branch but mass loss drags them back across the HR diagram to the blue side and 
they become Wolf-Rayet stars evolving along the WR sequence. At the upper end of the mass scale the stars 
become WR stars straight from the main sequence. For the LMC and SMC samples these transitions happen 
at higher initial masses than in the solar case. 

Figures 3.6, 3.7 and 3.8 display various details of these stars after core carbon burning or neon burning 
in some cases. There are many trends with mass and metallicity. Mass ranges of certain features tend to 
become spaced out over wider regions at lower metallicity. For example in figure 3.6 at solar metallicity 
there is a sharp changeover in radius, there are supergiants and then with the increment of 1M the stellar 
radius drops to around that of the Sun, while for the LMC and SMC there is a smoother transition. Figure 
3.7 displays the final mass and the final mass-loss rate of the stars. These show that in the most massive 
stars with hydrogen envelopes there are extreme winds. Higher final masses for stars only become possible 
at lower metallicities owing to the reduced mass-loss rates. The maximum final mass-loss rates increases 
with decreasing metallicity. This is because with less mass lost on the main sequence the stars spend longer 
ascending the supergiant branch and become more luminous which increases the mass loss from these stars. 
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Figure 3.3: Solar metallicity HR diagram. On these diagrams we plot all our models. Selected models have 
their mass indicated the lighter grey lines are the oth^Qmodels. 
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Figure 3.5: SMC metallicity HR diagram. On these diagrams we plot all our models. Selected models have 
their mass indicated the lighter grey lines are the oth^2 m °dels. 



3.6 CHARTING THE CHANGE BETWEEN SN TYPES: OBSERVATION 



Figure 3.8 looks at the internal composition of the different stars. What we must consider carefully here 
is when the changeover between type lb and Ic might occur. The helium rich progenitors become more 
prevalent at low metallicity and there seems to be a very clear definition between the two. Although at solar 
metallicity even with models at every integer mass it is difficult to resolve the helium rich progenitors that 
might give rise to the lb SN, from those deficient in helium that give rise to the Ic. However to complete the 
type I progenitor story we must look at binaries in chapter 5. 

Figure 3.9 shows estimates of the remnant and ejected masses if the stars were to go SN. These are obtained 
with the binding energy of the star. We assume a neutron star is formed at the centre of the star after core 
collapse of mass Mch = 1.44M . This produces about 10 46 J of energy from the release of gravitational 
binding energy in neutron star formation. We then assume a hundredth of this energy is transferred into the 
envelope by some unknown mechanism. The current suggestion is the transfer occurs via neutrinos released 
when the proto-neutron star forms that are thermalised within the envelope and the mantle of the core. 

We integrate the binding energy of the star from the surface towards the centre until we reach 10 44 /. The 
envelope outside this region is ejected with the remaining amount forming the remnant. 



If we have a remnant mass M rem > 2M the remnant is a black hole. This is a very rapid method to estimate 
the remnant and ejected masses quickly. It must be remembered however that the remnant masses are at best 
a lower bound as we use a very simple model for a complex event. From the figures it is possible to see that 
there are only very large remnants possible at low metallicities. Otherwise commonly there are black holes 
that are a few solar masses and neutron stars. Interestingly there is a maximum possible ejected mass in SN. 
This is important to compare to the calculated ejected masses from SN observations such as those of Hamuy 



With these models we can choose which features we need to consider in deciding the SN type a certain 
progenitor produces. But first we must investigate how SNe are separated into their various types. 

3.6 Charting the change between SN types: Observation 

Type I SNe do not have hydrogen in their spectra and are split into types la, lb and Ic. Type la are ther- 
monuclear events from low-mass stars (Mjuitiai < 8M ). They are identified by the presence of silicon lines 
but they are not core-collapse events so are not of interest to us here. Types lb and Ic are SN with WR-star 
progenitors (Ensman & Woosley, 1988). Comparing their spectra we find lb SN to have helium while Ic do 




(3.1) 



(2003). 
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Figure 3.6: The comparison between the different metallicities of the final radius, surface temperature and 
luminosity of our models versus initial mass. 
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Figure 3.7: The comparison between the different metallicities of the final mass and mass-loss rate versus 
initial mass 
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Figure 3.9: Estimates of the remnant and ejected masses from our models. 
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not and this is how they are distinguished. An example of a type lb SN is SN 19841 (Leibundgut, Phillips & 
Graham, 1990) while an example type Ic is SN1987M (Filippenko, Porter & Sargent, 1990). From the pres- 
ence or absence of helium in the spectrum we could jump to the assumption that lb progenitors must have 
helium while Ic progenitors do not. However it's not that simple. In lb SNe the helium lines do not become 
prevalent until after a few weeks from maximum light and simulations show that helium can be difficult to 
see if it is only present in trace abundances (Wheeler et al., 1987). There is also some evidence that Ic SNe 
do have helium present in the progenitor (Clocchiatta et al., 1996). We also know that the light curve and 
spectra of a SN depends on the size of the star. Helium giants will give rise to quite different SN to those 
from helium dwarfs. In all the models presented we find a small proportion of helium. So unless all Ic SNe 
are the results of a helium giants after a common envelope phase where they lose their helium envelope, we 
must determine which stars become lb or Ic SNe. Binary stars models are required for a full understanding. 

Type I SN types are poorly defined. They are defined by what they lack rather than what they have. For 
both they do not have hydrogen or silicon lines while Ic also has no helium. Indeed observations include a 
large number of SNe that supposedly fit into this category, including the newly identified hypernovae with 
energies ten times that of normal SNe, SNe that have been associated with Gamma-ray bursts and also fainter 
SNe events. The main reason for this is that there are not enough high-quality observations of type lb and 
Ic events to recognise the different types of SNe. However there are groups around the world that will be 
collecting vast samples of data on core-collapse events to improve our observational understanding (Gal- Yam 
et al., 2004; Kaiser et al., 2004). In this work we can tackle this problem from the theoretical end, describing 
what is possible from our current understanding of stellar evolution. This should help classify and recognise 
new SN categories in future. 

Type II SNe are better defined by what they do have. Most importantly they do have hydrogen in their 
spectra. The main types are IIP, IIL, lib, Iln and Ilpec. Type IIP have a long plateau phase where after the 
peak luminosity the SN luminosity remains constant owing to the photosphere remaining at a constant radius. 
A well observed SN of this type was SN1999em Leonard et al. (2002). The plateau occurs because while the 
ejecta expand outward the ionisation front moves inwards in mass. The photosphere is located at the point 
of hydrogen reionisation. A massive hydrogen envelope is required for the progenitor. Red supergiants fit 
the requirement well. 

IIL SNe have a linear decay in their light curves after maximum light mainly powered by radioactive decay 
of 56 Co. An example is SN1980K (Buta, 1982). Single-star progenitors of this type of SN are probably 
initially more massive than the IIP progenitors, close to the boundary of the type I SNe and have undergone 
more extreme mass-loss. We should note that IIL SNe are also possible from binary evolution. Type lib are 
a step closer to type I SNe. These objects were first discovered looking like a type II SN. However, after a 
few weeks, they change to type lb objects. Therefore the objects have only very small hydrogen envelopes. 
The classic example of this type is SN1993J (Filippenko, Matheson & Ho, 1993), one of the three SN with 
an identified progenitor. 
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Type Iln are distinguished by narrow lines in their spectra, SN1998S was a SN of this type (Fassia et al., 
2001). The naiTow lines indicate slow expansion velocities. Either this is due to a low explosion energy or 
a dense circumstellar environment where the ejecta is tamed by its interactions with circumstellar material. 
Either this is material ejected from a common envelope phase or a strong wind from the progenitor star. 
From radio observations of SNe it is possible to put limits on the strength of the wind required to be between 
10 -5 and 10 -3 M Q yr _1 . There are indications that Iln SNe are related to other SN types with features that 
are common to the other type II SN such as a plateau or linear lightcurve. It is also possible to split Iln SN 
into their two possible progenitor systems. For example a IIP SN with a low explosion energy would have 
a plateau phase with narrow lines which would have a P-Cygni profile which is evidence of an expanding 
hydrogen envelope. However if the ejecta interacts with a dense circumstellar environment the lines will not 
have the P-Cygni profile. 

We should note there are also impostor type Iln events. One class of these are luminous outbursts from 
luminous-blue variable stars or WR stars, examples of such events are 1961V and 1997bs. These are not SN. 
The other class is type la supernovae in dense circumstellar environments. An example is SN2002ic where 
a type la SN occurred in a hydrogen rich environment. 

Finally we have Ilpec, where we put the type II SN flotsam and jetsam. SN1987A is the most common 
example and can give some idea of what makes a SN peculiar. The progenitor star, Sk -69 202, was a blue 
supergiant, there is evidence for various mass loss events over the time of evolution and there is now a large 
amount of evidence that points towards a binary merger of two stars (Podsiadlowski, 1992). The surrounding 
environment is highly complex and the star was relatively compact and massive. Also the core structure may 
be quite unique due to the merger. 



3.7 Charting the change between SN types: Theory 



From stellar models we can draw some conclusions about SNe. However the only way to know for sure what 
type of SN a model produces would be to model the SN itself, the expansion of the envelope, hydrodynamics 
of the core, detailed nucleosynthesis and interactions with circumstellar material. We would have to simulate 
the light curve and spectra. Such models exist but there are a number of parameters, such as the amount of 
mixing that can be varied to produce different results. There are some analytical models but care must be 
taken in their use. For example Hamuy (2003) and Nadyozhin (2003) found very large ejected masses that 
are at odds with our models and the direct limits deduced by Smartt et al. (2003,2004). 

To predict the SN a model gives rise to we use analytic methods detailed below and specific details about 
the models. We also use the ejected mass and remnants mass estimated from equation . Then we compare 
to observations. The best observations are of SN which have identified progenitors. There are only three 
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observed progenitors to date and there are another six where observations exist to put limits on the properties 
of the progenitor stars (S mailt et al, 2003). 

We do not produce a theoretical scheme to name SN. Such a scheme would have limited use and it is 
preferable to estimate what we will find in observations so we use the observational scheme as a starting 
point. To estimate what we might see we need a number of details from our models. They are mostly those 
that have already been plotted in figures 3.6, 3.7, 3.8 and 3.9. But there are some details we need that are not 
directly available from our models so we must make some analytical estimate of their value. 

One of these is the remnant and ejected masses as discussed earlier. However to determine whether a type II 
SN has a plateau phase or not we use the following expressions from Popov (1993), 



A = 11.74 



4 34 M 



3/2 T 4 
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In these equations K0.34 is the opacity in multiples of 0.34gcm~ 2 , 7?o,soo the radius in multiples of 5OO/? , 
7ion,5054 the hydrogen ionisation temperature in multiples of 5054K and E44 the explosion energy in multiples 
of 10 44 J. The mass, M 10 , is the mass ejected in the SN in multiples of 10M Q . Our calculated values for 
plateau lengths agree with values from observations. In all our calculations we set K0.34 = £44 = 7i on ,5054 = 1- 
It is thought that the explosion energy is not a constant between supernovae so there will be some scatter 
around our values. However A has the greatest dependence on the energy and not the plateau length. For 
a plateau phase to occur A»l. Comparing models to these limits we set the limiting value A > 2. The 
plateau length must be sufficient that it can be observed. We take this limit to be 30 days. A is a variable that 
defines the strength of the plateau phase and it gives the limit that there must be enough mass in the envelope 
compared to the radius of the envelope. 

We list the information we use for SN type determination in table 3.7. If we think a model might give rise to 
a IIP SN we denote it as HP8o,3.o- This means that it will be a type IIP SN with a plateau length of 80 days 
and the second number is the value of A as defined above. Or for type Ibc we denote it as Ib3Co 95 or Ibo.iCo.2 
with the subscript on the b equal to the mass of helium present and the c subscript the helium mass fraction. 

From our example models we have calculated these factors and listed them in tables 3.2, 3.3 and 3.4. For 
each model we have listed the values relevant for all subtypes. By doing this it is easier to decide over which 
mass ranges different SNe might occur. We start by calculating whether the star is type IIP. If it is not then 
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Table 3.1: Theoretical SN Classification Scheme. 



Label 


Indication 


I 


No hydrogen present. 


II 


T-TvHttictpti nrp*iPTit 
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Thermonuclear explosion mechanism. 
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If type I, x is the helium mass fraction. 


n x,y 


If the progenitor model has a strong pre-SN wind, x is the mass-loss rate in lO _5 M yr _1 
and y is the total mass lost during evolution in M & . 


(i) 


Pre-SN radius, R < RQ. 


(ii) 


R Q <R< 10R q . 


(hi) 


10R q <R< lOOfl . 


(iv) 


lOO/?0 <R< lOOO/?0. 


(v) 


R > lOOO/?0. 



we describe it as IIL this seems to occur when the hydrogen envelope mass is less than about 6M when 
we apply the limit from the value of A. This is a little larger than the 2M Q used by Heger et al. (2003) but 
is similar to the upper limit of 5M determined by the lightcurve models of Swartz, Wheeler and Harkness 
(1991). However the plateau length is still quite large therefore these SNe could be the low luminosity IIP 
SNe that have been observed by Zampieri et al. (2003) rather than IIL SNe. These SNe also form black holes 
in their centres so the SN might be different. For type lib SN we will say these occur when the hydrogen 
envelope mass is less than 1M . Simulations of lib SNe indicate that while there is some uncertainty the 
value is of this order (Woosley et al., 1994; Shigeyama et al., 1994; Utrobin, 1994). 

Progenitors of type II SNe eject most of their mass as the star ascends the red supergiant branch. Because 
of this, during the SN, there could be substantial mass in the circumstellar environment in the most massive 
type II progenitors. 

While we do not specifically identify any model as a Iln SN we list the final mass-loss rate and total mass lost 
in the winds during the evolution if the pre-SN wind is stronger than lO _5 M yr _1 . This enables evaluation 
as to whether these stars will be observable by radio observations or the strength of any Iln signature in the 
SNe. SNe Iln probably have an underlying SN of another type. For now we separate type lb and Ic by 
choosing so that if y sur face > 0.5 we have a lb SN, otherwise we have a Ic. 
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Table 3.2: Theoretical SN details for solar metallicity models with overshooting. 



M 


Solar 


/M 




10 


nP 9 l, 4 .3L6.2,7.7(iv) 


11 


nP 98 , 4 . 5 L 6 . 8 , 8 . 6 (iv) 


12 


HPl03,4.5L7.2,9.3(rv) 


13 


nPi09,4.6L7.6,10.2(iv) 


14 


nPn 5 ,4. 7 L8.i,ii.o(iv) 


15 


IIPl20,4.2L8.1,11.7(v) 


16 


HPl25,4.2L8.4,12.4(v) 


17 


HPl30,3.9L8.4,13.o(v) 


18 


nPl34,3.6L8.3,13.5(v) 


19 


nPl37,3.3L8.2,13.90) 


20 


IIPl40,3. 1^8.1, 14.3( v ) 


21 


HPl42,2.7L7.7,14.5(v) 


22 


IIPl43,2.3L7.1,14.3 n 1.3,6.o(v) 


23 


HPl39,1.9L6.4,13.3 n 1.7,7.2(v) 


24 


IIPl33,1.4L5.3,12.in2.2,8.8(v) 


25 


IIPl25.0.9L4.0,10.6n3.o,10.5(v) 


26 


nPll5,o.3Ll.8,9.3bo.81 n 4.2,13.2(v) 


27 


Ib().96Co.98 n 3.4,15.9(ii) 


28 


Ibl.l5C0.96 n 2.6,16.9(ii) 


29 


Ibo.89Co.98 n 2.4,17.9(ii) 


30 


Ibo.46Co.89 n 2.4,19.2(ii) 


31 


It>0.16C0.30 n l. 2,20.7(0 


32 


Ibo.l9Co.24(i) 


33 


Ibo.23Co.24(i) 


34 


Ibo.26Co.25(i) 


35 


Ibo.28Co.26(i) 


36 


Ibo.28Co.26(i) 


37 


Ibo.29Co.26(i) 


38 


Ibo.29Co.26(i) 


39 


Ibo.29Co.26(i) 


40 


Ibo.28Co.26(i) 


45 


Ibo.28Co.27(h) 


50 


Ibo.28Co.27(i) 


55 


Ibo.29C0.27 n l.l,44.3(i) 


60 


Ibo.29C0.27lli. 1,48.7(0 


70 


Ibo.25Co.26lll.3,56.5(ii) 


80 


Ibo.43Co.25ni.8,63.o(ii) 


100 


Ibo.25Co.26ni.3,86.6(ii) 


120 


Ibo.29Co.27Hl. 1,108.4(0 


150 


Ib().30Co.27 n 1.0,139.7(i) 


200 


Ibo.28Co.28(ii) 
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Table 3.3: Theoretical SN details for LMC metallicity models with overshooting. 





LMC 


/M Q 




10 


nP 9 2,4.7L6.3,8.l(iv) 


11 


nP 9 9,5.2L 7 .0,9.l(iv) 


12 


HPio5,5.4L 7 .6,io.o(iv) 


13 


HPlll,5.7L 8 . 2 ,10.9(iv) 


14 


nPii5,5.5Lg.3,11.5(iv) 


15 


nPl21,5.2L 8 .7,12.3(iv) 


16 


HPl25,4.9L8.6,12.9(i v ) 


17 


IIPl30,4.7L8.8,13.6(v) 


18 


IIPl34,4.5L8.9,14.2(v) 


19 


IIPl38,4.3L8.9,14.8(v) 


20 


HPl42,4.lL8.9,15.2(v) 


21 


HPl44,3.7L8.6,15.6(v) 


22 


HPl45,3.4L8.4,15.4ni.0,4.7( v ) 


23 


HPl42,3.lL8.0,14.8 n 1.3,5.6( v ) 


24 


IIPi39 ! 2.5L7.2,13.8 n 1.7,6.9(v) 


25 


HPl35,2.oL6.3,12.9 n 2.2,8.2(v) 


26 


IIPl29, 1.41-5.2,1 1.7n 2 .9,9.9(v) 


27 


IIPll9,o.9L3.7,10.on3.7,11.9(v) 


28 


IIPl05,0.3Li.8,8.lbo.60 n 4.8,14.3(v) 


29 


IIP 9 5,0.lLo.7 ; 7.obo.l9n5.8,16.o( v ) 


30 


IIP 9 4,0.lLo.6,6.9bo.l6 n 7.3,16.7(v) 


31 


HP69,0.5Lo.6,7.ot>0.1 1115.7,17.3 (iv) 


32 


nP59,0.3Lo.2,7.obo.06 n 4.0,18.2(iii) 


33 


nb .oin 2 .5,i9.i(ii) 


34 


Ibl.66C0.99 113.0.20.3(0 


35 


Ibl. 83 Co.99 112.9,20.7® 


36 


Ibl.48C0.99 n 2.6,21.7(i) 


37 


Ibo.6lCo.93 n 2.2,23.4(ii) 


38 


Ibo.47Co.84lll.9,24.3(ii) 


39 


Ibo.3lCo.66 n 2.6,25.2(i) 


40 


Ibo.2lCo.47H2.1,26.2(i) 


45 


Ibo.23Co.27Hl.O,30.l(i) 


50 


Ibo.3lCo.28 n 1.2,34.2(i) 


55 


Ibo.3lCo.31 n 1.4,38.7(i) 


60 


Ibo.3lCo.30 n l. 4,43.6(0 


70 


Ibo.30Co.30 n 1.6,51.2(i) 


80 


Ib().30C0.29 n 1.8,58.7(ii) 


100 


Ib().68C0.28ll2.4,71.3(ii) 


120 


Ibo.80Co.29 n 3.1,83.l(ii) 


150 


Ibl.l5C0.34ll5.0,99.8(iii) 


200 


Ibl.85C0.48ll5.0,130.8(iv) 



93 



3. 7 CHARTING THE CHANGE BETWEEN SN TYPES: THEORY 



Table 3.4: Theoretical SN details for SMC metallicity models with overshooting. 



M 


SMC 


/M 




10 


nP 9 2,4.9L6.3,8.2(iv) 


11 


nP99,5.7L 7 .i,9.2(iv) 


12 


HPl05,6.()L7.7,10.2(iv) 


13 


HPl09,6.2L8.1,10.9(rv) 


14 


IIPll4,6.oL8.4,11.7(iv) 


15 


nPl21,5.7L 8 .7,12.5(iv) 


16 


HP 125 ,5.81-9.0,13. 3(iv) 


17 


HPl30,5.6L9.2,14.l(rv) 


18 


HPl34,5.5L9.4,14.8(iv) 


19 


nP 13 9,5.3L9. 5 ,i 5 . 4 (v) 


20 


HPl42,5.lL9.6,16.o(v) 


21 


HPl46,4.8L9.5,16.4(v) 


22 


HPl46,4.5L9.4,16.3( v ) 


23 


HP 144,4.()L9.0, 15.7 n l . 1 ,4.5 ( v ) 


24 


IIPl43,3.6L8.5,15.2ni.4,5.5(v) 


25 


HPl41,3.3L8.2,14.7 n 1.7,6.4(v) 


26 


IIPl37.2.5L7.1.13.6n2.3,7.9(v) 


27 


nPi33 J 2.oL6.1,12.7 n 2.9,9.4(v) 


28 


nPi22,1.3L4.5,10.8 n 3.8,11.5(v) 


29 


IIPll4,0.8L3.3,9.6n4.8,13.2(v) 


30 


HPlOO,0.3Ll.5,7.6bo.43 n 5.9,15.6(v) 


31 


IIP99,0.2Ll.3,7.4bo.38 n 7.8,16.3(v) 


32 


IIP99,0.2Ll.3,7.4bo.35 n 9.6,16.8( v ) 


33 


nPg 8 ,0.2Li.2,7.2bo.32ni 1.6,17.5 (v) 


34 


IIP97,0.2Ll.l,7.lbo.29 n 13.7,18.l(v) 


35 


IIP96,0.lLl.o,6.8bo.26 n 14.8,18.7(v) 


36 


IIP96,0. 1 Lo.9,6.9bo.24ril7.0,19.3 (v) 


37 


HP95,0.lLo.8,6.7bo.2ini9.8,19.9(v) 


38 


IIP94,0.lLo.8,6.5bo.l9 n 22.8,20.5(v) 


39 


HP93,0.lLo.6,6.4bo.l5 n 24.8,21.2(v) 


40 


IIP91 ,0.oLo.5,6.3bo. 1 1 n26.8,21 .9 (v) 


45 


IIP41 ,o.oLo.o,6.obo.oi 04.0,25.4(111) 


50 


Ibl.46C0.99 n 5. 3,29.7(0 


55 


lb 1.00C0.94111. 0,33.3(11) 


60 


Ibo.90Co.89 n l. 9,37.3(ii) 


70 


Ibo.42Co.45(ii) 


80 


Ibo.76Co.43(ii) 


100 


Ibo.92 Co.33n3.2,59.5(i) 


120 


Ibo.92Co.33 n 3.2,79.5(i) 


150 


Ibl.28C0.37 n 4.6,97.0(il) 


200 


Ib2.59 c 0.54ll5.0, 1 1 9.9 (ill) 
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3.8 NAMING THE SN FROM THE PROGENITORS 



Table 3.5: The Plateau lengths from S-AGB and Extreme S-AGB Progenitors. The first set are models 
without convective overshooting, the second set are models with convective overshooting. 



Initial 




Initial 




Mass 


Without 


Mass 


With 


/M Q 


Overshooting 


/M Q 


Overshooting 


9.5 


HPl05,3.oL8.1 


7.5 


IIP90,2.lL6.1 


10 


HPl04,3.oL8.6 


7.7 


IIP92,2.2L6.3 


10.1 


HPl05,3.oL8.7 


7.9 


HP93,2.4L6.5 


10.2 


HPl06,3.2L8.8 


8.1 


IIP903.2L6.7 


10.3 


IIPl07,3.3L8.9 


8.3 


IIP913.4L6.9 


10.4 


HPl07,3.4L9.o 


8.4 


IIP9L3.6L7.0 


10.5 


IIPl05,3.8L9.1 


8.5 


HP9OA1L7.I 


10.6 


IIP103A5L9.2 


8.6 


HP89.4.5L7.2 


10.7 


HP99.10.7L7.3 


8.7 


nPg6,6.5L6.1 



3.8 Naming the SN from the progenitors 



We list the details of our models in terms of the theoretical SN scheme in tables 3.2, 3.3 and 3.4. From these 
we list the ranges of SN types with metallicity in table 3.7. Comparison with observed ratios is not possible 
because we have not included binaries although the ratio of type II to type Ibc is in the correct region for 
spiral galaxies. However we should first notice that there are not enough IIL SNe to fit the observed value of 
10-50% of II SNe being IIL (Gaskell, 1992; Cappellaro et al., 1997) observed today. To resolve this without 
binaries we would have to lower the mass range for IIL SNe or include some He SNe in the IIL group. There 
are also candidates for the type Iln SN. In the range of stars we have labelled as IIL there are stars that have 
extremely high mass-loss rates before the SN. These stars provide a very dense circumstellar medium that 
absorbs some of the energy in the expanding envelope. Although, as with all the determinations, there is a 
continuum of objects. 

Looking at the type I SN we see that for the models all hydrogen deficient SN retain some helium. However 
the surface abundance does vary. Therefore it is possible that this is the source of the difference between the 
two SN types. All the type I progenitors are very small and therefore tightly bound so only a little mass will 
be ejected to provide a little display so some of the more massive type lb stars may not be observed at all. 

In table 3.6 we use the values from tables 3.2 and 3.3 to constrain the mass ranges or the progenitors. Note 
that from this list we have removed SN1987A and SN1993J due to these systems being confirm as in binaries. 
SN1980K we discuss below due to the disparity between out models and the observational limits. With the 
remaining SNe we are able to fit them reasonably well with out models. Assuming our LMC models are at 
roughly the same as half solar metallicity. 
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3.8 NAMING THE SN FROM THE PROGENITORS 



Table 3.6: Comparing values from models to observations of progenitors. Adapted from the table in Smartt 
et al (2003) and details from Smartt el al. (2004). 





Observed 




Observed 




Observed 


Theoretical 


Theory 




SNe 


Observed 


Stellar 


Initial 


Ejected 


Mass Range 


Ejected 


SN 


Type 


Z/Z 


Type 


Mass /Mq 


Mass/M Q 


Mq 


Mass /M Q 


2002ap 


Ic 


0.5 


WR? 


<40 


2.5-5 


38-40 


<6 


1980K 


IIL 


0.5 


9 


<20 




25-28 


6-12 


2003gd 


IIP 


0.5 


M 


8 




8 


6.5 


2001du 


IIP 


1 


G-M 


< 15 




8.7-15 


6-12 


1999em 


IIP 


1-2 


K-M 


< 15 


5-18 


8.7-12 


6-10 


1999gi 


IIP 


2 


G-M 


< 12 


10-30 


8.7-12 


6-10 



Table 3.7: Relative populations for SN types from theoretical models with convective overshooting. The 
figures are obtained from tables 3.2, 3.3 and 3.4 applying our determination on SN types. 



SNe Type 


Solar 
Mass Range 

/M Q 


Population 


LMC 
Mass Range 

/Mq 


Population 


SMC 
Mass Range 

/Mq 


Population 


He 


7.5-8.6 


21.6% 


6.5-7.5 


22.5% 


6-7 


24.0% 


IIP 


8.7-22 


63.3% 


7.6-25 


67.9% 


7.1-27 


68.7% 


IIL 


23-26 


3.4% 


26-28 


1.7% 


28-35 


2.7% 


lib 




0% 


29-32 


1.6% 


36-40 


1.3% 


lb 


27-30 


2.5% 


33-39 


1.8% 


45-200 


2.8% 


Ic 


31-200 


8.9% 


40-200 


4.4% 




0% 



Returning to the lower end of the mass range we take the minimum for the IIP region as the upper mass from 
the He SNe from extreme Super-AGB stars. We separate out the He SN for the specific reason that these 
stars may have strange SN properties. Table 3.5 shows that they have similar properties to our IIP models. 
However with overshooting the case for these having a plateau weakens as the models with the longest period 
to SN after second dredge-up will lose a large amount of mass, reducing the envelope mass. These models 
could in fact undergo a type IIL SN, as suggested by Swartz, Wheeler and Harkness (1991). 

If we include overshooting then S-AGB stars have smaller hydrogen envelopes. Making it is easier for 
mass-loss to remove the envelope to the point where a type IIL SNe occurs. Table 3.5 also provides further 
evidence with the lower values of A for the overshooting models. Therefore to search for the progenitors of 
electron-capture SNe we should search for IIL SNe as well as IIP. 

There is one SN that provides a tentative hint that electron captures do give rise to IIL SN, 1980K. This was 
a type IIL SN and even though the progenitor was not detected the limit placed on its mass was M < 2OM . 
There are three possible interpretations. First that the stellar winds must be high to remove the envelope of 
a star between 15 and 2OM , enough for a IIL to occur, second that a binary interaction occurred removing 
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the mass or third that the initial mass was in fact quite low and we have an electron capture event. 

There is a further piece of information to consider. Montes et al. (1998) found evidence for a change in the 
mass loss 10,000 years before SN over a period of about 4000 years. This information fits most comfortably 
with the binary or electron capture models. From figure 3.2 we can see that a 9.5M star (7.5M with 
overshooting) fits the required change in mass loss history. A strong wind over 10,000 years would have 
reduced the envelope mass to a great degree and a type IIL SN would be the most likely outcome. The 
derived mass-loss rate is about 2 x 10 M Q yr _1 . This would be enough for both types of star to remove 
sufficient envelope for the SN to be IIL in the required time. 

The binary model is the only other viable alternative. The problems with the binary model is that the mass 
loss must be relatively constant and must occur at a very specific time close to the time of SN. Thus we would 
require a stable non-conservative mass transfer phase, something that may be difficult but not impossible to 
achieve. 



3.9 Conclusion 

In this chapter we have shown that single stars provide a restricted number of possibilities for types of SNe. 
It is also apparent that the definition of which progenitors give rise to which SNe is rather arbitrary. This 
problem arises because there are not clear cut changes but a continuum. At high metallicities this continuum 
is compressed in the range of masses over which it occurs and is only resolved at lower metallicities. 

We have also described electron capture SN progenitors that evolve rapidly just before the SN owing to 
the event of second dredge-up. There may also be some current observational evidence that these electron 
capture events do occur. Only observational searches for progenitors will find these stars. 

Stellar wind mass-loss rates are probably not as large as commonly assumed. It is likely that the require- 
ments for a IIL SN are not as restrictive as previously thought or that the transition objects are currently 
misidentified as Ilpec or Iln before the IIL behaviour dominates. 

Wolf-Rayet stars are quite similar at solar metallicity. Most end up as WC stars. With lower metallicity the 
reduced mass loss slows down WR evolution through the WR types. This means the fate of the most massive 
stars changes with metallicity. 

With this detailed study we must next move on to study how changing the mass-loss scheme changes the 
ratios we derive and how these change over a wider metallicity parameter space. 
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Chapter 4 



Mass-Loss Prescriptions & Progenitor 
Populations 

"The thing's hollow - it goes on for ever - and - oh my God - it's full of stars!" 
Dave Bowman, 2001: A Space Odyssey, by Arthur C. Clarke. 



4.1 Introduction 

There are many examples of similar studies to that in the previous chapter, Woosley, Heger & Weaver (2002) 
or de Donder & Vanbeveren (2003) for example. However each one uses a different mass-loss prescription 
and no one has investigated the difference in models between different prescriptions using the same stellar 
evolution code. Because of this and the basic uncertainty in the mass-loss rates from massive stars a detailed 
comparison of how the evolution of massive stars depends on the mass-loss prescription is required. There- 
fore in this chapter we will compare a number of different prescriptions to determine the effect on the SN 
progenitor population. We will use observations to constrain which mass loss prescription we should use. 

The mass-loss rates we study are broken down into pre-WR and WR rates. We also study the effect of scaling 
the mass-loss rates with metallicity by a factor of (Z/Z )' ! for various n rather than the usually assumed 0.5 
taken from the work of Kudritzki, Pauldrach & Puis (1987) and Kudritzki et al. (1991). This scaling arises 
from the assumption that stellar winds are line driven and with lower surface opacity at lower metallicity 
there are weaker winds. However while there is agreement that mass-loss scales in this form there is a 
range of suggested values for the exponent. We have tested the sensitivity to changes in this scaling by 
producing a grid with n = 0.7, the exponent suggested by Vink et al. As expected a larger exponent increases 
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variation with metallicity: the type II/Ibc SN boundary varies with metallicity to a greater degree. Previous 
authors have commonly assumed that this scaling applies at all stages of evolution even though there is some 
evidence that it varies with spectral type. This seems likely because the exponent assumes the winds are line 
driven while mass loss from giants is driven by some quite different, unestablished mechanism. If we do not 
scale the mass loss for giant stars (log(r e ff/K) < 3.7) the result is that the minimum initial mass for WR star 
formation decreases with decreasing metallicity. This is contrary to observations of WR stars in the LMC 
and SMC (Massey, Waterhouse & DeGioia-Eastwood, 2000). This provides evidence that mass-loss from 
giants must scale with metallicity. 

We also investigate whether the scaling with initial metallicity should be included for WR stars and the 
effect this has on our results. We should note that in this thesis we study a broad range of metallicity, high 
(1 x 10~ 3 < Z < 0.05) and low (10~ 8 < Z < 10" 3 ). We include an overlap to see how different schemes in 
the two regimes match. 

First we present the various mass-loss rates we use. We next discuss the effect of including the metallicity 
scaling for WR mass-loss rates and then compare a large number of mass loss rates at high and low metal- 
licity. Finally we present the results from our preferred mass loss rates and compare the differences caused 
by removing convective overshooting. 



4.2 The Mass-loss Rates 

In order to produce the various mass-loss prescriptions we combine various rates. Some of these are designed 
for low metallicity only and this requires us to be careful in which regions we apply the rates. 

The pre-WR mass-loss rates we use are 

• JNH: de Jager, Nieuwenhuijzen & van der Hucht (1988). Despite being old these are accurate and 
detailed (Crowther, 2001). They give a mass-loss rate for a specific place in the HR diagram. They are 
quite complex with over 20 terms in the equation but the extra detail increases the accuracy. 

• NJ: Nieuwenhuijzen & de Jager (1990). These rates are based on the same data as for JNH but they 
are simpler with only 3 terms in the mass-loss equation. This does provide some advantages. They are 
much easier to extrapolate from than the JNH rates and thus do not have such a limited coverage of 
the HR diagram. This rate generally leads to greater mass loss than the JNH rates. 

• Vink: Vink, de Koter & Lamers (2000) and (2001). These are theoretical rates that match well to 
observed values listed in de Jager, Nieuwenhuijzen & van der Hucht (1988). The model is only valid 
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over a metallicity range of 4, < (Z/Z ) < 3. 

• Krki: Kudritzki (2002). These are theoretical rates for the low metallicity regime. They seem to 
indicate a break in the mass-loss rates at around Z = 10 5 . We have only used this rate in our low 
metallicity regime. 

The Vink and Krki rates are only valid for OB stars so they are combined with the NJ and JNH rates. 
When X sur f ace < 0.4 and r e ff > 10 4 K a star becomes a WR star. The WR mass-loss rates we use are 

• NL: Nugis & Lamers (2000). These rates are the most recently derived WR mass-loss rates from 
observations. There are some theoretical rates for WR stars but they do not agree well. These also 
include the rate for WO stars used by Dray (2003). 

• Vbev: Vanbeveren, De Loore & Van Rensbergen (1998a) and Vanbeveren (2001). This rate is ex- 
tremely simple depending only on luminosity and does not consider the change of mass loss with WR 
stellar type. 

• WL: Wellstein & Langer (1999). These are effectively the rates from Langer (1989) but lowered by a 
factor of 3 to account for clumping of the winds. 

The analytical expressions used for these mass-loss rates are detailed in appendix A. We combine them in 
various ways to produce results similar to those in chapter 3 but with coarser resolution in mass because we 
here have a wider range of metallicity. The combinations we use at high metallicity are 

• HG: NJ+WL, from Heger et al. (2003). 

• VB: NJ+Vbev, based on the work of Vanbeveren (2001). 

• NJ: NJ+NL, for comparison with the JNH rates. 

• NJV: NJ+Vink+NL, to observe the effect of the Vink rates. 

• JNH: JNH+NL, the rates of Dray (2003). 

• JNHV: JNH+Vink+NL, to observe the difference of the Vink rates. 
At low metallicity we also use 

• NJK: NJ+Krki+NL, to observe the effect of the Krki rates. 

• JNHK: JNH+Krki+NL, to observe the effect of the Krki rates. 
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JNH Mass-loss Rates / log(M G /yr) NJ M a s s — I :) :■; s Nine:-: ,/" k;g: H. ; ,/./i ; 




3.8 4.0 4.2 4.4 4.6 4.8 5.0 3.6 3.8 4.0 4.2 4.4 4.6 4.8 5.0 

Temperature / log (T/K) Temperature / log (T/K) 



Figure 4. 1 : The JNH (left) and NJ (right) mass-loss rates on the HR diagram. The NJ rates are calculated for 
a 1OM star. Changing the mass of the stars for which the rates are calculated has only a small effect on the 
rates. 

4.3 Models details 

In this work we use 15 zero-age main-sequence models that have masses from 10 to 2OOM with a uniform 
composition X = 0.75 — 2.5Z and Y = 0.25 + 1.5Z, where X is the mass fraction of hydrogen, Y that of 
helium and Z is the initial metallicity. We have also created models with alternate compositions: X = 0.7, 
Y = 0.3 — Z and X = 0.76 — 3Z, Y = 0.24 + 2Z. With constant hydrogen abundance the variation with 
metallicity is accentuated. For example 2nd dredge-up occurs at higher masses with increasing metallicity. 
Changing the primordial hydrogen abundance from 0.75 to 0.76 as in the second case only produces small 
changes in our results. 

Unless stated otherwise all models are evolved with convective overshooting included as described in chapter 
2 and are evolved with 499 mesh points. Some problems were encountered in modelling the highest mass 
stars at high metallicity. There is a small range of mass where large helium stars with high mass-loss rates 
are formed and these cannot be evolved. The physical parameters move outside of the region covered by 
the opacity tables, to high temperature and low density. Linear extrapolation of the tables does not work 
and the density continues to decrease. The star is expanding and becoming a high mass helium giant. For 
most of the sets of models we have interpolated over mass to complete the grid of models. If the models 
are at the extremities of our range we have not extrapolated into that region. The method to produce models 
for this regime is to limit the mass-loss rate during the WR phase to a constant value that gives models of 
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similar mass and only slightly more massive than neighbouring models and similar to the expectations from 
interpolation. This indicates the interpolation is an acceptable compromise. The range of initial masses is 
around 80— 1OOM at solar metallicity. This range decreases with increasing metallicity and increases with 
lower metallicity. When Z < 0.001 the problem region moves above the range of mass studied. These stars 
indicate that there may be a maximum WR mass-loss rate because the mass loss for these stars is so severe 
it causes them to almost evaporate. Alternatively it could mean that the mass loss should be greater and that 
the low-density regions actually should be ejected. To properly investigate this theoretical rates for WR stars 
are required to better understand this feature and discover if it is real or not. 



4.4 Supernovae Type and Remnant Determination 



We use similar methods to those in chapter 3 but for comparison we use simpler schemes. For a star with 
a model at the end of core carbon burning we say it will become a type II SNe if there is any hydrogen 
in its envelope. Furthermore we know observationally that type IIP SNe have retained a good proportion 
of their original hydrogen. This hydrogen causes the light curve to follow a plateau powered by a moving 
hydrogen ionisation front. Type IIL SNe only retain a small fraction of their hydrogen and the light curve 
decays linearly. We adopt the scheme of Heger et al. (2003) with type IIP SNe until the mass of hydrogen in 
the star drops below 2.OM when type IIL SNe occur 

When there is no hydrogen the SN will be type I. As already detailed it is difficult to discriminate between 
type lb and Ic. Therefore we adopt a scheme to use the surface helium abundance for the moment. If 
^surface > 0.7 a type lb SN occurs; if 0.3 < F sur f ace < 0.7 we denote it as a type Ibc due to uncertainty and 
finally if y sur f ace < 0.3 type Ic. These are arbitrary boundaries but are simple enough that we can gain some 
understanding of the relative populations. 

We adopt from Heger et al. (2003) a method of determining the strength of a type Ibc SN to quantify the 
number of SNe that have no observable display and are thus unseen. This allows us to estimate ratios of type 
II to type Ibc SNe. The largest progenitors are thought to give no display because the core is so massive that, 
even with a large explosion energy, nothing escapes the forming black hole. The exception would be if a 
jet driven SNe occurs that makes a black hole and produces an observable display (MacFadyen, Woosley & 
Heger, 2001). We take the ranges from Ensman & Woosley (1988) for the strength of the display: 

• ^Core(He) > 15M , no display. 

• 15M Q > Mcore(He) > 8M , faint SN. 

• 8M Q > Mcore(He) > 5M G , possibly faint SN. 

103 



4.4 SUPERNOVAE TYPE AND REMNANT DETERMINATION 



• Mcore(He) < 5M , bright SN. 



Determination of the remnant formed at the heart of a SN is a black art. The physics is extreme and poorly 
understood. However there are many prescriptions which rely on the conjecture that a more massive core 
leads to a more massive remnant. We consider here three types of stellar remnant, white dwarfs, neutron 
stars and black holes. 

White dwarfs come in three flavours. They can be made of helium (a He WD), a mixture of carbon and oxygen 
(a COWD) or a mixture of oxygen and neon with some sodium and magnesium (an ONeWD). HeWDs are 
formed by very low mass stars or more probably in binary stars and are not of interest to us here. COWDs 
are also not formed in SNe but are the remnant cores of stars that lose their envelopes on the asymptotic 
giant branch (AGB). At solar metallicity stars in the mass range 0.8 to 7M (no convective overshooting) 
produce a COWD while larger stars, up to 1O.5M develop an ONe core during the early AGB. Some of 
these lose their envelope and form an ONeWD. The most massive undergo a SN driven by electron capture 
on to magnesium. The dividing line between these two is of great interest but is somewhat uncertain and 
depends on the mass-loss rate and the details of the thermal pulses in AGB stars as discussed in the previous 
chapter. 

In our selected mass range (10 — 2OOM ) all stars undergo a SN and form a neutron star or black hole. To 
distinguish between neutron stars and black holes there are a number of different prescriptions. We compare 
that of Heger et al. (2003) with our own. 

Heger et al. (2003) use a simple system based on the helium core mass at the end stage of evolution. This 
is useful because the helium burning life-time is slower than for the later stages of burning so the core mass 
increases by only a small amount in the late burning stages. However the method relies on setting the values 
for different remnants by using simulations of SNe. The values used are 

• Neutron Star: M Cor e(He) < 8M . 

• Black Hole (by fall back): 8M < M Cor e(He) < 15M . 

• Black Hole (directly): M Cor e(He) > 15M . 



These values are based on the SN models of Fryer (1999). By using this scheme we can compare directly 
with Heger et al. (2003). We should also note that SNe that form a black hole directly also have no display, 
while those SN forming black holes by fall back will be faint unless a jet driven SNe occurs producing a 
black hole and an observable display. 
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4.5 WR Mass-loss rates: scaling with initial metallicity 

The first matter to investigate is that of the scaling of WR mass-loss rates with metallicity. WR stars have 
higher mass-loss rates than OB stars of the same mass. The main difference between OB and WR stars is 
that the WR stars have very little hydrogen or no hydrogen at all with the surface dominated by helium or 
in later WR types carbon and oxygen. This affects the number of lines for radiatively driven winds. This 
provides evidence that the composition affects mass-loss. 

It is known that mass loss is affected by metallicity in OB stars. Our question is, 'if OB winds depend on 
initial metallicity do WR winds too?' This would mean that the metals in the initial metallicity have a strong 
effect on the mass-loss mechanism of WR stars, perhaps providing evidence that WR winds are still linked 
to radiative driving. There is observational evidence for this too, Crowther et al. (2002) see a scaling of 
(Z/Z ) - 5 from observations of WR stars in the LMC. From these tentative observations it is necessary to 
determine how this affects the evolution and final end states of massive stars. For this study we use the rates 
of JNH, Vink and NL. We then vary the metallicity scaling of the NL rates, setting the exponent to 0, 0.5 and 
0.7. 

Massey, Waterhouse & DeGioia-Eastwood (2000) determined the limits for WR stars to form at differ- 
ent metallicities by studying the turn-off ages of clusters with WR stars. The limits are M m j n > 25M , 
Mnin,LMC > 3OM and M m j n ,sMC > 7OM . They also derive lower limits for the formation of WC stars but 
these are all limited by small number statistics. The WC minimum masses are higher than the WR minimum 
masses and increase at lower metallicity. Our models must therefore agree with this result. But we must con- 
sider the fact that the observations are for a very small sample of stars and thus have large errors associated 
with them. Some of the observed stars are also likely to be binary stars. 

In figure 4.2 we present the mass just before the SN for the stars. Figures 4.3 and 4.4 show the hydrogen 
and helium mass fractions and figure 4.5 shows the final helium core mass that we can use to determine the 
remnant type. One thing that is evident from these figures is that the effect of scaling the mass loss follows 
some general trends but can also be non-linear in some places. 

The final masses of the pre-SN provide our main limits on the mass loss at metallicities above solar. In 
general without scaling the remnant masses are higher and all WR stars that go SN are more massive than 
5M . Therefore, according to the above prescription, this means that we see no bright Ibc from single stars. 
When the scaling exponent is 0.7 rather than 0.5 we find that the resultant masses are slightly lower than that 
with 0.5. 

When we move our attention to metallicities lower than solar the picture becomes quite different. The 
remnant masses without the scaling are consistently lower than when it is included. There is little difference 
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Figure 4.2: The final mass versus initial mass for different WR mass loss scaling laws at different metal- 
licities. The numbers in the legend are the value of the exponent, n, in the initial metallicity scaling 

M(Z)=M(Z & )(Z/Z y. 
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between the results achieved when using different values for the scaling exponent. At the lowest metallicities 
the effects are weak due to the later on set of the WR phase because pre-WR mass loss is also lower at lower 
metallicity and gives rise to initially more massive WR stars so there is more mass for the WR mass loss to 
remove. 

A similar tale is told by the mass fractions of helium and hydrogen. A larger range of WR stars at lower 
metallicity retain more hydrogen than their higher metallicity cousins with the same mass. This indicates 
quite strong evolution of WR type ratios with metallicity especially if the mass-loss rates do scale with initial 
metallicity. Indeed below SMC metallicity there will be very few WC stars. At solar and above nearly all WR 
stars end their evolution as WC or WO stars with low surface helium while at lower metallicity a growing 
fraction retain more helium at a high mass fraction, Y > 0.9. The graphs also show strong evidence for two 
populations of WR SN progenitors. This could be the type lb and Ic division. 

The helium core mass in figure 4.5 is important for determining remnants. At the highest metallicities, first 
we can see with the scaling there are no black holes although without the scaling there are black holes formed 
by fall back. Also we can see black holes cannot form directly until below solar metallicity. Then black holes 
formed become more massive at lower metallicity if the scaling exists. Therefore limits on this scaling could 
be possible from measuring the black hole initial mass function and how it depends on metallicity. However 
this is likely to be a very difficult to measure. 

We have compared the ratios of WR types from these models against observations. The models with no 
scaling give a constant ratio of WC/WR and WC/WN which is at odds with observations. The models with 
scaling do roughly match the trend as can be seen in figure 4.23. Therefore observation leads us to include 
the extra scaling for WR stars. Models with the extra scaling of mass loss with initial metallicity do agree 
with the observations of Massey, Waterhouse & DeGioia-Eastwood (2000) that the limit for WR and WC 
stars increases with lower metallicity. However the direct measurement of WR mass loss by Crowther et al. 
(2002) provides the most convincing argument. A measure will also be possible from the ratio of type I/II 
SNe and the ratio of type Ib/Ic against composition. This requires more observations. In all the mass-loss 
schemes we now include the scaling of WR mass loss with metallicity using an exponent of 0.5. 



4.6 Comparing Prescriptions at High Metallicity 



We now compare a number of commonly used mass-loss schemes. We also include a set of models that have 
zero mass loss. We also show results with JNH rates with different scaling of the mass-loss with metallicity. 
JNH5 has a scaling of (Z/Z ) a5 while JNH7 has a scaling of (Z/Z e )° 7 . 

Two of the schemes we study are similar to those used by Heger et al. (2003) and Vanbeveren, De Loore 
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Figure 4.3: The final hydrogen mass fraction on the surface versus initial mass for different WR mass loss 
scaling laws at different metallicities. The numbers in the legend are the value of the exponent, in the 
initial metallicity scaling M(Z) =M(Z e )(Z/Z Q ) n . 
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Figure 4.4: The final helium mass fraction on the surface versus initial mass for different WR mass loss 
scaling laws at different metallicities. The numbers in the legend are the value of the exponent, in the 
initial metallicity scaling M(Z) = M(Z Q )(Z/Z Q ) n . 
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Figure 4.5: The final helium core mass versus initial mass for different WR mass loss scaling laws at different 
metallicities. The numbers in the legend are the value of the exponent, n, in the initial metallicity scaling 

M(z)=M(z G )(z/z & y. iiq 
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gure 4.6: The final mass versus initial mass for different mass-loss schemes at different metallicities. 
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& Van Rensbergen (1998a). For the HG mass-loss scheme we use the NJ rates for pre-WR and the rates of 
Wellstein & Langer (1999) for the WR phase, while for the VNB rates we use the NJ rates which are similar 
to the VNB rates for OB stars. We have experimented with their full mass-loss scheme but find it non-trivial 
to implement due to the large disparity in their mass-loss rates. The NJ rates do provide a similar magnitude 
of mass loss and prevent the occurrence of giant stars at solar metallicity with Mt, i > —9.5. This is argued 
to be necessary because red giants of this luminosity are not observed in the solar neighbourhood. We do 
use their WR mass loss rate despite its extreme simplicity because it does not consider the effect of varying 
surface composition or WR type indicated by Nugis & Lamers (2000). For all WR mass loss we scale the 
rate with the initial metallicity such that M WR (Z) = M(Z Q )(Z/Z Q ) - 5 . 

We shall analyse the following details from our models, the final mass in figure 4.6, the surface hydrogen 
and helium mass fractions in figures 4.7 and 4.8, the final helium core mass in 4.9 and the length of the type 
II plateau phase in figure 4. 10. From these we can draw conclusions on the differences between the mass 
loss rates. 

The final mass data shows interesting agreements between the different stars. In all the diagrams there is 
a clear split between those with NJ and those with JNH rates. Around solar metallicity and above this 
behaviour is clearly seen when Minimal < 3OM . The JNH rates give rise to higher mass progenitors than 
the NJ rates for type II SNe. As the metallicity drops the difference alters so that the NJ rates give rise 
to more massive progenitors. The reason for these differences is the extra structure in the JNH rates that 
effectively leads to lower mass loss for stars than the NJ rates. However, at lower metallicity, the stars are 
more luminous and this can lead to more mass loss for the most massive stars because they end up in the 
region with greatest mass loss as shown in figure 4.1. The mass loss becomes severe despite the metallicity 
scaling. This demonstrates the only problem of the JNH rates which were measured at solar metallicity. At 
different metallicities the position of the evolution tracks of the star shuffles due to different opacities and 
different initial hydrogen and helium composition. Near to solar this effect is small but it escalates as the 
metallicity becomes more extreme. The NJ rates, because of their simplicity, suffer less from this problem 
although it is their simplicity that makes them less accurate at near solar metallicities where mass loss is too 
great. The only way to resolve these problems is further observations of the mass loss from stars at different 
metallicities and theoretical mass-loss rates for all stars. 

The second region of interest in these diagrams is the region when Minimal > 30M© at around solar metallicity 
and above. They are all WR despite having different pre-WR mass-loss rates. At the highest metallicity there 
are two populations. Those using the NL rates have higher final WR masses than those using alternative WR 
mass loss rates. As the metallicity moves toward solar we find that the VNB rates move into agreement with 
results from the others while the HG rates continue to have lower pre-SN masses. At the lowest metallicity 
the difference becomes small due to the metallicity scaling. The most discrepant result seems to be that 
of HG because of the estimate of lower WR masses. The mass-loss rate they may be using has not been 
reduced enough from the rates of Langer (1989) however the result is not discrepant to more than a factor of 
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2 at solar metallicity. The final point to note is that by comparing the JNH results, we can also see that the 
larger exponent also has its greatest effect at the lowest metallicities. The addition of the Vink rates to the 
NJ and JNH rates makes only small changes. 

In 4.7 and 4.8 it is easier to see over what ranges we see the WR stars. The mass of hydrogen and helium has 
a tight correlation to the mass fraction which means we can deduce the amount of hydrogen and helium in 
the models from these plots. From the hydrogen plot at solar metallicity and above there is a clear separation 
between the NJ and JNH rates with different minimum masses at which all the hydrogen is removed from a 
star. Immediately below solar metallicity the results start to agree more closely until the relation is inverted 
at the lowest metallicities where JNH produces more hydrogen deficient stars than NJ. From studying the 
helium graphs at solar and above once the hydrogen has been removed all the models agree to a large degree 
but with the HG and VNB results showing the lowest helium fraction due to greater mass loss. Below solar 
metallicity the difference in the helium fraction grows until below SMC metallicity when the JNH rates all 
have more helium than all the NJ based rates. This is because the JNH rates remove more hydrogen from 
the star than the NJ rates so in the latter much of the surface mass is still hydrogen leading to lower helium 
abundance. 

The helium core mass details in figure 4.9 confirm the above conclusions. It is apparent that there are different 
maximum metallicities at which black holes will form directly, with the HG rates needing a metallicity of 
below Z = 0.01 for this to occur. They also require a lower metallicity, lower than solar, to form black holes 
by fall back. The other schemes have direct black holes or fall back black holes at metallicities higher than 
these rates. This again indicates that the HG rates are higher on average than others. For all rates very 
massive helium cores are not possible until SMC metallicity and below so in these regions we find the most 
massive black holes. In view of the lowest helium core masses only the HG and VNB rates have WR stars 
below 5M Q that will give rise to bright lb events. 

Our final figure 4.10 of plateau duration for IIP SN splits the JNH and NJ rates again with a clear distinction 
between the two. The JNH rates give rise to a maximum plateau length of 140 days while the NJ rates give 
rise to a maximum of 1 10 days at solar and above. These increase below this range due to the greater amount 
of mass left in the stars as the stellar winds weaken. The longest plateaus are achieved by the JNH7 rates. 
The remaining rates all tend to agree at the lowest metallicities. In table 4. 1 we compare these values with 
those from our two main contending schemes. 

Both of these agree with the results but do favour JNH with the longer possible plateau lengths. We should 
bear in mind that more massive envelopes lead to more luminous plateaus so there is also a selection effect 
on these statistics and they are for a small sample of SN from Hamuy (2003). Also in calculating the plateau 
lengths we have used a constant explosion energy. Despite this the results for the JNH models agree best with 
these observations, apart from JNH7 at low metallicity. While NJ rates only agree at the SMC metallicity and 
below, they lack plateau lengths of around 140 days. We really need further observations to provide details 
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Figure 4.7: The final surface hydrogen mass fraction versus initial mass for different mass-loss schemes 
different metallicities. 
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Figure 4.8: The final surface helium mass fraction versus initial mass for different mass-loss schemes at 
different metallicities. 
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Figure 4.9: The final helium core mass versus initial mass for different mass-loss schemes at different metal- 
licities. 
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Figure 4.10: The plateau length duration for SN with IIP light curves versus initial mass for different mass- 
loss schemes at different metallicities. 
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Table 4. 1 : The average values for the plateau duration of the IIP SN light curve. 



Source 


Average Value 


Shortest 


Longest 




/days 


/days 


/days 


Nadyozhin 


105 ±15 


80 


140 


Hamuy 


131 ±21 


105 


1 HI 

171 


JNHV Solar 


104 


90 


140 


JNHV LMC 


105 


90 


140 


JNHV SMC 


104 


90 


140 


NJV Solar 


97 


90 


112 


NJV LMC 


98 


90 


120 


NJV SMC 


104 


90 


132 



of how the plateau length changes with metallicity to choose one of these schemes over the other based on 
these details, especially since they will also have strong winds that might give a type Iln appearance if they 
are in the early WR phase. 

In summary there are some important discrepancies between mass-loss schemes and there are some important 
similarities. Which prescription is best is testable by SN and WR observations. The rates which seem 
to provide the best results are the JNH5 and JNHV rates. The latter are our preferred rates at the higher 
metallicities covered here and we shall investigate them further below. 



4.7 Comparing Prescriptions at Low Metallicity 



We now move our analysis to the lowest metallicities. These are closer to population III stars where Z = 
however they still have evolution similar to the stars at higher metallicities. The biggest divergence of 
evolution of low metallicity stars occurs around Z = 10~ 12 . At this level the CNO abundance is too small 
for the CNO bi-cycle to support the star from collapse so the star continues to contract until helium ignites 
and forms enough carbon that is then processed to nitrogen and oxygen by hydrogen burning which can then 
support the star. The mass-loss rates for these stars however remain largely unknown. For the hot OB stars 
we can use the Vink rates although we are far below the regime they were derived for. The only rates that 
do exist are those of Kudritzki (2002). These extend far down into the extremely low metallicity regime. 
At these low metallicities there are few WR stars so their rates are not as important here. In this section 
therefore we must compare the NJ and JNH rates combined with the theoretical rates for OB stars. 

To compare we use the final masses and final amount of hydrogen present as in figures 4. 1 1 and 4. 12. In 
the final masses we see that when Z = 0.001 the results all agree well but then upon lowering metallicity we 
find they exhibit peculiar behaviour. The 2OOM JNH models retain most of their mass as do stars below 
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about 70M©. The stars in between however lose a large proportion of their initial mass and hydrogen down 
to the lowest metallicities. This is because the star's position on the HR diagram is moving into either the 
region where the mass-loss rates are greatest or outside the valid region of the JNH rates. This means that 
the results are inaccurate because we are not applying the rates in the correct way. Either we could move the 
rates so a star starts at the same position on the rates diagram or just use the NJ rates that provide similar 
results for lower masses. Interestingly there are only minimal differences between the Krki and Vink rates. 

Over this range of metallicity we use the NJ rates supplemented with the Krki rates. The NJ rates avoid 
the problems of the JNH rates and the Krki rates are derived specifically for this range of metallicity. Un- 
fortunately it is difficult to obtain observational constraints so we are reduced to using our preferred rates. 
Note however that we have shown that the differences are minor. Therefore we are now able to combine our 
preferred mass-loss rates from the two grids for coverage over the entire range from Z = 10~ 8 to 0.05. 



4.8 Our preferred mass loss prescriptions 

Combining the above rates we can produce grids that represent the evolution of single SN progenitors over a 
large range of metallicity. While we could plot any details of the progenitors we concentrate on a small group 
and briefly examine their shape on the grids. We display the final mass, luminosity, surface temperature, 
radius, plateau length, helium core mass, SN type, estimated ejected and remnant masses, hydrogen envelope 
mass, mass of hydrogen and helium and discuss each result in detail and compare the overshooting and no 
overshooting cases. After presenting the results of this test we describe some tests of our models from SN 
type ratios, SN progenitors and WR ratios. In appendix ?? we show more details of these models and in 
appendix ?? we show the same graphs for the NJV scheme at high metallicity, the second scheme on our list. 

Our maps tend to broadly agree with those of Heger et al. (2003). Owing to our specific coverage in metal- 
licity, we obtain more structure within the map. All of the structure is due to the mass-loss prescriptions we 
have used in producing the models. It is important to note however that while the broad nature of these maps 
are correct some of the fine structure would in practice be blurred out by different helium content, rotation, 
magnetic fields and other mixing processes. We assess our maps for each variable individually. 

4.8.1 Luminosity, Figure 4.13 

In the no overshooting diagram we can see in the upper left corner stars that undergo 2nd dredge-up. For 
models at lower metallicity and for models with overshooting we do not go to low enough masses to resolve 
this region. We will cover this in more detail later. The main feature in the diagram is the structure in the 
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upper right corner at high metallicity. As mass increases we see the last few red giant progenitors are quite 
bright but are then replaced by WR stars and the pre-SN luminosity drops abruptly. The luminosity increases 
again as higher mass WR stars form. At the highest initial masses these WR stars once more become less 
massive and therefore less luminous. This structure becomes less pronounced at lower metallicities owing to 
the lower mass-loss rates. At low metallicity there is limited mass-loss and the main structure at the lowest 
metallicities is due to changes in the nature of burning reactions at the lowest temperatures. 



4.8.2 Spectral type, Figure 4.14 



High metallicity results show a strong difference between giant stars and the hotter WR progenitors. As 
metallicity decreases this change over becomes more gradual. Towards the middle right section there are 
slightly cool WR stars. These are those WR stars for which to evolve the model, we have had to adopt a 
limited constant mass-loss rate. 

Low metallicity has the remnants of the previous behaviour in the upper right corner however there is an 
interesting change at the lowest metallicity with the slightly hotter progenitors. This is because the opacity 
is low enough to prevent the formation of red giants. 



4.8.3 Radius, Figure 4.15 



The radius grids reflect the results of surface temperature. It should be noted that with convective overshoot- 
ing the WR models are generally more compact and in both cases WR stars become less compact towards 
lower metallicities. The low metallicity grid also shows the interesting increased compactness at the lowest 
metallicity for some stars. 



4.8.4 Final Mass, Figure 4.16 



The final mass diagram has very similar structure to the luminosity diagram and again reflects how mass 
loss is more severe in the high-metallicity region and drops off as metallicity decreases. The low metallicity 
graph shows less structure and is essentially similar to that with no mass loss. 
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4.8.5 Plateau Length, Figure 4.17 

The plateau length at high metallicity shows the limited region for IIP SNe. The plateau length is on average 
higher for the no overshooting case than when it is removed. This provides a statistical test which probably 
indicates that most stars do have some form of convective overshooting in their cores. 

Low metallicity grids tell a very different story with extremely long plateaus possible. However these are 
at too high a red-shift to observe today. However since their duration is longer there is more of a chance 
of spotting one. At the lowest metallicity here the plateau lengths drop off quickly due to the increasing 
compactness of the stars. 

4.8.6 Hydrogen Envelope Mass, Figure 4.18 

The high metallicity grid mimics the plateau length diagram closely. This also shows there are maxima at 
different metallicities of how massive the envelopes are and therefore how much envelope can be ejected. 
Although on top of this we must also consider the outer mantle of the core that is also ejected. The low 
metallicity plot shows that at lower metallicities there is more hydrogen and so smaller cores. 

4.8.7 Ejected Mass, Figure 4.19 

The mass ejected by the SNe is also interesting because it demonstrates that a WR star ejects very little 
mass because it is denser and so more tightly bound than a red giant progenitor. Also it is possible to see, 
along with figure 4.18, that for red giants a large fraction of the core mantle is ejected. In the low-metallicity 
diagrams the compact nature of the stars gives a more tightly bound envelope so less material is ejected. 

4.8.8 SNe Type, Figure 4.20 

The SNe type diagram demonstrates the different nature of the grids with and without overshooting well with 
a clear shift in the minimum mass for a type I SN. The line for the changeover from IIP to IIL SNe is quite 
constant. However the line separating IIL and type I SNe moves more and there is a cut-off in type I SN at 
just below Z = 0.001. The separation between lb and Ic is quite difficult to observe. At the lowest initial 
masses for type I SNe there is likely a band of lb SNe whose width increases with decreasing metallicity. 
However at LMC metallicity and below there are no more Ic SNe. The low metallicity grids are structurally 
less captivating. 
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Figure 4.13: Contours of pre-SN luminosity. Contours are in log(L/L ). Left no overshooting, right over- 
shooting. 
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Figure 4.14: Contours of pre-SN spectral type. Left no overshooting, right overshooting. The effective 
temperatures taken for each spectral type are as follows, r e g-(0) > 27000K, r e g(B) > 11000K, r e g-(A) > 
7600K, r eff (F) > 6000K, r eff (G) > 4900K, r eff (K) > 3600K and T eff (M) < 3600K. 
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Figure 4.15: Contours of pre-SN radis. Contours are in log(/?/7? ). Left no overshooting, right overshooting. 
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ure 4.16: Contours of pre-SN mass. Contours are in (M/M ). Left no overshooting, right overshooting. 
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Figure 4.17: Contours of plateau length. Contours are in days. Left no overshooting, right overshooting. 
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Figure 4.19: Contours of mass ejected in SN. Contours are in (M/Mq). Left no overshooting, right over- 
shooting. 
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Figure 4.20: Contours of SNe type. Green - IIP, orange - ILL, red - lib, blue - Ic, purple - Ibc and brown - lb. 
Left no overshooting, right overshooting. 
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Figure 4.21: Contours of remnant of helium core mass. Numerical contour values are helium core mass in 
Mq. Green - direct black holes, yellow - fall back black hole, red and other colours - neutron stars. Left no 
overshooting, right overshooting. 
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Figure 4.22: Contours of SN remnant mass. Contours are in (M/M ). Left no overshooting, right over- 
shooting. 
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Using Figure 4.21 and the details of section 4.4, we can identify the mass regime where SN Ibc become 
unseen. Those SNe that form black holes directly have no display, unless a jet driven SNe forms a black 
hole and produces an observable display. Unseen SNe affect the ratio of type II/Ibc SNe observed. At below 
solar metallicities all massive progenitor type Ibc SNe are unseen. We find very few bright type Ibc SNe in 
our sample because only at the highest metallicity do WR stars have helium cores of less than 5M Q . 



4.8.9 SNe Remnant, figures 4.21 & 4.22 



In shape figure 4.21 is similar to those of Heger et al. (2003) although we do not go to zero metallicity. The 
change between the no-overshooting and overshooting cases is minor with only the transitions moving down 
in mass while the shape remains the same. Noticeably at high metallicity and mass with overshooting it is 
possible to see that neutron stars are formed rather than black holes. This is similar to the findings of Heger 
et al. (2003). The transitions from neutron star to black hole are at similar positions in Heger et al. (2003) 
although our value with overshooting is lower by 5M e . 

Comparing to the energy method for determining the remnant in figure 4.22 we find the same transition point 
between neutron stars and black holes just below 25M Q without overshooting and 21M with overshooting. 
The structure is also quite similar although it is difficult to decide when direct black hole formation may 
occur rather than fall back on to a black hole. This would require a denser remnant and our 1OM remnant 
line seems to follow the He core mass remnant region for direct black holes fairly well. While looking at 
these diagrams it is sensible to ask where different types of gamma-ray bursts (GRBs) may occur as discussed 
by Heger et al. (2003). It is thought that GRBs require a progenitor with a small radius (a few R & ) so the 
jets that form from the material accreted on to the central black hole can punch their way through the stellar 
envelope. If the stellar envelope is too large the jets dissipate and lead to a normal SN. Therefore any of 
our models that loose their hydrogen envelopes to become WR stars and also form a black hole may also 
be the progenitors of GRBs. It will be interesting to determine how many stars that fit this requirement give 
rise to an unseen SN compared to those that give a bright display. This could lead to an observational check 
if it were known how many GRBs have an associated SN. The remnant masses here should be considered 
lower limits to the possible remnant mass. The final mass increases during the course of carbon burning so 
if the model is truncated before then a misestimation might occur, for instance the feature at high masses at 
Z = 10~ 6 . 
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4.9 Discussion 

We have now presented the results from our many models and presented our preferred results. But how do 
these preferred results and all the other tables compare with observations? We have already had one test, the 
comparison of theoretical plateau lengths against observations. However we now discuss more observational 
tests and implications for the remnant mass population. 

4.9.1 Convolution with an IMF - SNe Ratios 

Observations tell us that in SOa-Sb galaxies, the ratio of Type II supernovae over Type Ibc supernovae is 
3.82±2.71. In Sbc-Sd galaxies, it is between 6.14±3.96 (Cappellaro et al., 1997). At solar metallicity we 
find the ratio without overshooting to be 8.2, and with overshooting to be 9.0 for a Kroupa, Tout and Gilmore 
(1993) IMF. We get slightly lower values with the NJV rates of 6.9 and 7.8. These are higher than the value 
of 5 calculated by Heger et al. (2003) but we have a lower minimum mass for a star to go SN. Generally the 
theory agrees with observation but this basic picture must be altered at lower metallicities because most type 
Ibc SNe become faint or have no display thereby raising this ratio. However binary systems are likely to 
dominate Ibc supernova production at low metallicity and this complicates prediction of the SNe type ratio. 

4.9.2 Observations of SN Progenitors 

Smartt et al (2003) compiled details on the sample of well studied SN progenitors (shown in Table 4.2). We 
can calculate various details to check against these observations. Most match well with our overshooting 
models, with the JNHV models giving the best agreement. The one problem is that the various types of type 
II supernovae do not match our IIP predictions. This may be due these stars being binaries, for example 
recent observations (Maund et al., 2004) have confirmed the progenitor to 1993J was part of a binary. There 
are also strong arguments that 1987A was part of a binary (Podsiadlowski, 1992). Except for 1993J the 
ejected masses also match up well to those calculated from observations. As this list grows in length and 
detail more tests against our maps will become possible. 

4.9.3 WR Ratios 

In figure 4.23 we show the predicted WR type ratios from models with observations. Assuming constant star 
formation and using the IMF of Kroupa, Tout and Gilmore (1993) we determine how long our models spend 
in each of the WR types and then determine the relative populations. The theoretical lines are approximate 
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Figure 4.23: Comparison of predicted WR type ratios with observations. Lines are calculated from the grid 
above, 'x' are from calculations on the models in chapter 3 and '+' are from observations in Meynet et al 
(1994). 
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owing to the low resolution in initial mass of our grids. However the JNHV OS results do agree with values 
from chapter 3 the minor differences are due to the lower resolution in mass of our grids in this chapter. 
These models do agree well with the observations. The general trend is for the WC/WR and WC/WN ratios 
to agree better than the WR/O ratio. This is expected because binaries produce WR stars from a lower 
initial mass due to increasing the opportunities for mass-loss therefore increasing the number of WR stars 
relative to O stars. But the ratios between WR types would be constant since this depends mass-loss once 
a star becomes a WR star and being in a binary will only have a slight affect on this. All the mass loss 
prescriptions provide reasonable agreement however the JNHV with overshooting is closest overall if we 
also consider the models from chapter 3. 

4.9.4 Remnant Mass Population 

Figure 4.22 shows that, at the highest metallicities, the remnants are all of very low mass, due to a maximum 
remnant mass cut off, with neutron stars dominant. Then as metallicity decreases more massive remnants 
become possible. Very massive (M > 1OOM ) black holes can form in SN however this can only happen at 
low metallicity (Z < 10~ 3 ). The population of such black holes is very small because they are only formed in 
the most massive stars. At higher metallicities massive black holes are not formed because mass loss limits 
the size of the SN progenitor and thus the size of the remnant. Any black holes observed to be very massive 
in high metallicity environments must therefore have grown via accretion from a binary companion or must 
have formed from an earlier population of stars. 

4.10 Conclusion 

In this chapter we have presented results from numerous stellar evolution calculations with various mass-loss 
rates. We have highlighted how changing the mass loss affects the final outcome of these stars. We have 
also discussed in detail our favoured mass-loss prescriptions using the rates of de Jager, Nieuwenhuijzen & 
van der Hucht (1988) with the rates of Vink, de Koter & Lamers (2001) for pre-WR evolution. These we 
compare with models that have essentially the rates of Nugis & Lamers (2000) with the constant WO rate 
from Dray (2003). Except for the Vink rates which have their own scaling, we scale all the mass-loss rates 
by (Z/Z ) 5 . We have shown that the rates agree with three sets of observations. Although to constrain our 
work more we require more data in these observations to reduce the limited number statistics affecting the 
results. 
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Table 4.2: Comparing values from models to observations of progenitors. Adapted from the table in Smartt 



et al (2003) and details from Smartt el al. (2004). 
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Chapter 



5 



Binary 



Progenitors 



"First things first, but not necessarily in that order.' 
The 4th Doctor Who, Meglos. 



5.1 Introduction 

Even though astronomers would like to believe that stellar evolution is done and dusted nothing could be 
further from the truth. The previous chapter has illustrated one problem, mass loss, the other great unknowns 
are rotation, magnetic fields and duplicity. 

Rotation and magnetic fields are complex issues and affect the evolution of single stars although Maeder & 
Meynet (2003) find they might cancel each other out. A binary companion probably has the greatest effect 
on a star's evolution. Its presence provides the possibility of mass loss, mass gain and other binary specific 
interactions (e.g. irradiation, colliding winds, surface contact, gravitational distortion) and hence a greater 
scope for interesting evolution scenarios. SN2003gd only just fits with our theoretical models as none of our 
single-star models have a low enough luminosity for an exact match. Therefore we undertook this study to 
try and find lower luminosity type IIP SN progenitors from binary evolution. 

The basic physics of binaries is deceptively simple until we see the strange and bizarre situations that require 
complex physics to model. The classic example of a problem binary is Algol. In this system it appears the 
less massive star is more evolved than the more massive star. The explanation is that mass transfer occurred 
to swap the mass ratio. However when attempts are made to model this system not all the mass transferred 
is accepted by the secondary star and a fraction is lost from the system (Sarna, 1993). 
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In this chapter we first describe our treatment of binaries in the stellar evolution code. This also includes a 
discussion of what we haven't included in the code that produced our main results and the differences we can 
expect if we were to. Then we break up our binary models into similar mass ranges and discuss them in detail, 
highlighting changes in the evolution outcomes over these ranges. This study is similar to others performed 
by Podsiadlowski, Joss & Hsu (1992), de Donder & Vanbeveren (2003) and Izzard, Ramires-Ruiz & Tout 
(2004). The first of these is most similar to our own although we use a more up-to-date stellar evolution 
code. The others use a binary population synthesis which does not use full stellar evolution calculations but 
instead uses analytic expressions that are fit to the results of detailed models. This has the advantage that 
they are far more rapid but not all the details of the stars' evolution are known. 

To gain an overall picture of the evolution of SNe progenitors we need to have some idea of how duplicity 
affects them. We discuss the SN progenitors only possible from binaries and not single stars. Because the 
models also provide details of the evolution of massive binaries we discuss the possibility of X-ray binaries, 
from mass-transfer events with a compact remnant. 



5.2 Modelling Binaries 



To truly model a binary a code would have to follow both stars simultaneously. There is a problem however 
for systems with stars that do not have mass ratios close to 1 in that the process becomes costly in CPU time. 
The more massive star evolves with shorter time steps and this wastes CPU time on the evolution of the 
secondary with a longer evolution timescale. In our code we evolve one star at a time. The entire evolution 
of the primary, say 10,000 models, takes up less than 1000 models when we evolve the secondary. 

The best system would be to evolve one star until something interesting happens then switch back and forth 
evolving the stars to best save computer time. The system used for our models is a unique method for dealing 
with the secondary star in the system. 



5.2.1 Basic Details 



The first piece of the binary physics jigsaw is Kepler's laws. These detail the orbits of the planets around 
the Sun but they also apply to stars in orbit around each other. The two stars revolve around each other held 
together by gravity but cannot merge because of angular momentum. Modelling the stars as point masses in 
a circular orbit we set the centripetal force equal to the gravitational force between them, 



, GM { M 2 

\la(ti = 5 — =F graV) (5.1) 
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where /I = M l + ^ the reduced mass, M\ is the mass of the primary and more massive star, M2 the mass of 
the secondary star, a the orbital separation and ft) is the orbital period. From rearranging the equations and 
using the fact that ft) = ^ we arrive at Keplers law, 



where P is the period. During the evolution of the binary as a first assumption we can take the angular 
momentum to be conserved. If this is true then 



With some simple analysis (see appendix ??) we can find that mass transfer affects the orbit. For example 
if the more massive star is transferring mass to the less massive star then the orbital separation shrinks, 
however if the less massive star transfers mass the separation expands. This only occurs if mass transfer is 
conservative; if mass is lost from the system such as in a wind then the escaping material also carries away 
angular momentum. This leads to the widening of the orbit. For the most massive stars that have severe 
winds the orbit can widen significantly and prevent interactions. 



5.2.2 Dealing with the Secondary 

We could treat the secondary as a point source with no variables other than mass. This is acceptable for 
compact objects but is not satisfactory for a star of similar mass to the primary. Our solution was to use the 
single stellar evolution equations of Hurley, Pols & Tout (2000). These are mostly used in binary population 
synthesis calculations such as Hurley, Tout & Pols (2002). They are much faster than stellar evolution 
models, taking a fraction of a second to calculate the evolution of one star rather than tens of minutes. 

Using equations enables us to estimate the details of the secondary to check whether it fills its own Roche 
Lobe and calculate the winds from the star and the amount of angular momentum lost. Once the primary's 
evolution is complete we recalculate the secondary's evolution with the full stellar evolution code, using the 
mass loss and gain from the analytic calculations during the primary's evolution. When we reach the point in 
the secondary evolution when the primary undergoes a SN we reset the binary calculations using a compact 
remmant with the remnant mass from the primary's calculation. 




(5.2) 



(5.3) 
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5.2.3 Roche lobe overflow 

In a binary system the equipotential surfaces are not spherical but are distorted by the gravity of the com- 
panion star and by the rotation of the system. Solving for the equipotential surfaces in a corotating frame 
the shape of the equipotential surfaces can be found. If the stars are small compared to the size of the or- 
bital separation these surfaces are spherical but as the star grows in comparison to the separation the shape 
is distorted to a greater degree becoming more ellipsoidal. The main mass transfer event in binaries occur 
when the radius of the star is of the order of the orbital separation. A star's surface will eventually reach 
the LI Lagrangian point where the gravity of both stars cancels exactly. If a star expands beyond this point 
then material begins to flow towards the other star. The method to include this in the code is to define a 
Roche lobe radius, suc h that the sphere of this radius has the same volume as the material held within 
the Roche Lobe defined by the equipotential surface passing through the LI point. Therefore when the star 
has an outer radius greater than the Roche lobe radius it will transfer material onto the other star. We use the 
Roche lobe radius given by Eggleton (1983), 

= (54) 

a 0.6< ? i /3 +rn(l+< ? ; /3 )' 

where q\ = It is accurate to within 2% for the range < q\ < °°. When R[ > 7?li Roche lobe overflow 
(RLOF) occurs and we have mass transfer in our system. We use a second simple equation to calculate the 
rate at which mass is lost from the primary, 



Mm = FtMOKflj/tfu^Moyr- 1 , (5.5) 



where 



F(M\) = 3 x 10" 6 [mm(Mi,5.0)] 2 (5.6) 

and is chosen by experiment to ensure mass transfer is steady (Hurley, Tout & Pols, 2002). We use this for 
all RLOF calculations unless common envelope evolution occurs which requires a different scheme. 

Once we have taken the mass away from the primary we transfer it to the secondary but not all the mass 
is necessarily accreted. In the code the accreting matter has the same temperature and entropy as the star's 
surface. However this is a simplification since the accreted matter will cause the star to expand due to the 
increased mass of the star increasing energy production. The star will relax to a smaller radius but only on a 
thermal timescale, therefore we assume that the star's maximum accretion rate is determined by its current 
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mass and its thermal timescale. We define a maximum accretion rate for a star such that M2 < M2/TKH where 
Tkh is the Kelvin-Helmholtz timescale. If the accretion rate is greater than this then the additional mass is 
lost from the system so the star only increases in mass at a rate that is determined by its thermal timescale. 
If the secondary is of lower mass it has a longer thermal timescale than its more massive primary. Efficient 
transfer is only seen if the two stars are of nearly equal masses so the thermal timescales are similar and the 
mass ratio changes soon after mass transfer begins so the orbit widens. If the masses are dissimilar we find 
that most of the mass is lost from the system and the separation of the orbit increases but if the increase is 
not rapid enough a common envelope phase quickly occurs. 



5.2.4 Common Envelope Evolution 

If RLOF occurs, but doesn't arrest the expansion of the mass losing star, growth continues until the radius 
is greater than the binary separation and the secondary is engulfed in the envelope of the primary. This 
is common envelope evolution (CEE). When this occurs it is thought that the envelope is lost by some 
unknown dynamical process but the energy for the envelope ejection is transferred from the orbital energy 
thus shrinking the orbital separation in the process. There is a chance that the two stars may coalesce 
before the envelope is ejected. Evidence suggests that the progenitor of SN1987A underwent such a merger 
(Podsiadlowski, 1992). The more common outcome is the two stars ending up in a tighter orbit, commonly 
one helium star, and one main sequence star. If the secondary then evolves a second CEE phase can occur 
with the outcomes of a merger or a very compact binary. 

We require a simple prescription for this complex event. Commonly the energy required to eject the envelope 
is calculated and this is subtracted from the orbital energy to provide the new separation of the orbit. If there 
is not enough energy in the orbit to eject the envelope the stars merge. There is some uncertainty in the 
exact ratio of these two energies because there are other sources of energy that can be put into ejecting the 
envelope such as the energy from nuclear reactions and reionisation energy of the hydrogen in the envelope. 
Therefore an arbitrary constant is included to account for the other sources of energy. There are some other 
prescriptions based on angular momentum, for example Nelemans & Tout (2004). 

In our models when CEE occurs we initiate a CEE wind to remove matter as quickly as possible, mass-loss 
rates of 10 _3 M Q yr _1 are possible but such rapid mass loss can cause the evolution to be unstable. Therefore 
we use a rate of lO _4 M yr _1 . The higher mass-loss rate restricts any thermal response from the star and the 
envelope is lost in a short period of time, of the order of lOOOyears, the lower rate only changes the result 
slightly and allows more stable evolution of the star. 

To remove the mass we assume the energy comes from the orbit. For our model we assume that the secondary 
accretes no mass unless it is a compact object which is an acceptable approximation since our calculations 
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show most common envelopes occur in binaries with mismatched masses where the thermal timescales are 
quite different. 

To determine how the mass loss effects the separation we take the small amount of energy, SZibinding to 
remove the mass, 8M\ from the surface as 



GMM 

^binding — „ • P-') 

R\ 

We simplify the orbital energy and assume most of the orbital energy is stored in the core of the primary and 
the secondary. Therefore, by the Virial theorem, the orbital energy is 



_ GM lc M 2 

£orbit — Z ■ P-°) 

2a 



Differentiating this we find, 



GMi C M 2 

d£ OT bit = j dc?. (5.9) 

a 1 



Equating equation 5.7 and 5.9 and rearranging we find, 



a 2 Mi dMi 



da = (5.10) 

Ri M 2 Mic 



In this formalism AM\ is negative so the orbit shrinks. We stop CEE once we have reached the helium core 
mass before the CEE phase began, at this point we have removed the hydrogen envelope. This is slightly 
different to the standard approach to just remove the envelope and calculate the final separation using the 
assumption that binding = «CE(^orb,/ — £orb,()> obtaining the final separation from the final orbital energy. 
G!ce is a free parameter to fix the efficiency of the energy transfer. We could not use this process for a 
full stellar evolution code as we cannot just remove the envelope! In comparison our method is equivalent 
to assuming Oce = 1- Although as we do not remove the envelope in one go the effective result is that 
the effective Oce slightly less than 1 and our final orbit is larger than from the standard approach. Binary 
population synthesis codes also use another parameter Ace to quantify the effect of structure of the envelope 
on the binding energy. We do not have to consider this since with a detailed model we can calculate the 
binding energy accurately. But in our model we calculate the binding energy of the material begin removed 
which is perhaps closer to reality. 
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5.2.5 Compact Objects 

These are white dwarfs, neutron stars and black holes. There are changes for when we evolve with these 
objects. First the maximum accretion rate is derived from the Eddington accretion limit. For black holes we 
allow accretion ten times this limit to reflect the absence of a surface. 

The other main difference is that we allow compact objects to accrete matter during CEE. The phase is so 
short and the Eddington rate so constraining we find only little change in the mass of the compact objects. 

5.2.6 Other things: tides, wind accretion, SN, etc... 

We have also included wind accretion in the formalism of Hurley, Tout & Pols (2002). However we find 
it irrelevant in SN progenitor systems because a slow wind for a long time is required to see an effect and 
only AGB stars fit these requirements. WR stars may have the strongest winds but they are too fast for a 
secondary to accrete much material. Winds also carry away angular momentum from the system. We assume 
any mass lost takes away angular momentum from the orbit. 

We deal with SN very simply. We evaluate the ejected mass and if more than half the total mass of the 
binary is ejected the system is disrupted and the secondary evolves as a single star. Otherwise the secondary 
evolves with the remnant of the primary and we leave the separation at the separation of the last model for 
the primary. This gives us the worst case scenario for subsequent evolution. Based on observations of free 
neutron stars with high velocities it is thought that in SN neutron stars receive a kick . We have not included 
this in our simple model but it will increase the eccentricity of the post-SN system and disrupt more of these 
binaries. 

We have experimented with tides however we find they are numerically unstable because the code is evolving 
on a thermal timescale rather than a dynamical timescale near RLOF They do increase the chance of common 
envelope evolution but there is a competing effect if the mass transfer during RLOF is inefficient and large 
amounts of mass are lost from the system. Their inclusion is likely to increase the number of lb SN from 
stars that have lost their hydrogen envelopes. 

We have not included thermohaline mixing where heavier material is accreted on top of light material and 
then mixes into the envelope even if there is no convection at the surface. Most of our mass-transfer is ineffi- 
cient and occurs while the primary ascends the red giant branch and so still has a large fraction of hydrogen 
in the envelope. However, as we shall discuss, some stars below 15M undergo a second RLOF event when 
the star becomes a helium giant. This is likely to affect our results slightly because thermohaline mixing is 
thought to be quite efficient occurring on a thermal timescale (Kippenhahn, Ruschenplatt & Thomas, 1980). 
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Finally all these models are evolved at solar metallicity with convective overshooting and 199 meshpoints. 
The winds for the stars use the JNH rates for pre-WR evolution and the NL rates for WR evolution. 



5.2.7 Mass Ratio and Separation 

We are limited in our resolution in this study by computing time. On the SunGrid machines that we are able 
to use a single star takes one and a half hours to evolve and if we wish to have an answer in a week then if 
we use 20 computers simultaneously we have a maximum number of required models of about 2000 1 . With 
a binary population synthesis code this would take about 4 minutes. In mass ratio we choose our grid to be 
q = = 0.1, 0.3, 0.5, 0.7 and 0.9. In separation we study 10 15 < a/R Q < 10 4 . The lower limit is arbitrary 
and we assume any stars below this merge during the main sequence and evolve as single stars. The upper 
limit comes about from experiment to find the separation when there is no interaction. The initial separations 
we use are \og(a/R l3 ) = 1.5, 1.75, 2, 2.2, 2.4, 2.6, 2.8, 3, 3.25 and 4. This gives us fairly good resolution 
over the entire range of binaries. For a 1OM primary and with ^ = 0.5 this is a range of periods from 5 days 
to 82 years. Our primary masses are My =5, 6, 7, 8, 9, 10, 11, 12, 13, 15, 20, 25, 30, 40, 60, 80, 100, 120, 
150 and 200. We concentrate the resolution at the low end of our range to pick up the minimum mass SNe. 
The number of models in this grid is 2 x 5 x 10 x 20 = 2000 hitting our above target. In future more models 
could be run but for the moment we do not want to take too much time since the code is still a prototype. 



5.3 Evolution of Massive Binaries 

We will now describe the evolution of the binaries in our study. We will pay particular interest to the types of 
interaction and how these might affect the SN progenitors. The evolution of binaries is extremely non-linear, 
similar behaviour may occur over large ranges of masses and separation while a tiny change will cause a 
quite different outcome. 

The mass transfer interactions are our main interest in for the evolution of binaries, the primary losing mass 
and the secondary gaining mass. This will radically affect the outcome of evolution. It may change the SN 
type, prevent the SN or produce a SN where none would have occurred. 



This means 1000 binaries since there are two stars per binary. 
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5.3.1 5<M 1 /M Q <9 

In this range we only see a primary SN if Mi > 8M©. Any star below this mass evolves towards the AGB 
phase and if there is no interaction will lose its envelope during thermal pulses. The interactions that occur 
are roughly equal occunences of CEE and RLOF for the primary while secondaries experience more RLOF 
than CEE. This is because during primary evolution the masses are more mismatched so mass transfer is 
unstable while for the secondaries they have more similar masses to the remnants so mass transfer is stable 
in general. 

It is common for there to be two interactions during the evolution of the primaries. Once they lose their 
hydrogen envelopes by CEE or RLOF they become helium dwarfs then later after core helium burning they 
become helium giants expanding to the point where they can again interact. The second interaction is stable 
RLOF as the masses become roughly equal after the first interaction. After losing mass these stars form 
a COWD. The secondary may also experience this behaviour with the remnant from the primary. This 
behaviour only occurs if the mass is above about 3M Q otherwise the star never becomes a helium giant after 
the first mass transfer event. 

SNe that do occur depend on the mass of our stars. At 7M some SNe do occur but they are mostly lb 
and IIL when interactions have stripped the envelope and prevented 2nd dredge-up and the AGB phase. We 
therefore have helium stars more massive than Mch that may go SN. We have set the minimum lb mass at 
1.6M Q . From 8M Q and above we see many more SNe because the stars would normally go SN if they were 
single stars. Most SNe are type lb if an interaction occurs and IIP if there is no interaction. There is a small 
band of IIL SNe separating the two regions. If we only see lb and IIP and no IIL or lib then we assume that 
the stars in between the separation where the change over occurs give rise to IIL SNe. The closest binaries 
prevent SNe occurring by removing the hydrogen envelope at first dredge-up leaving helium stars too low in 
mass to undergo core-collapse. 

The secondaries are far more intriguing. Many of the secondaries in this range are not massive enough to 
undergo SNe if they were single stars. Despite this some of them do. When the mass ratio is close to one we 
find that mass transfer is very efficient because the thermal timescales of the two stars are similar. This can 
increase the mass of the secondary to a point where it becomes massive enough to go SN. We find these stars 
when q > 0.7 and at separations of a few lOO/? . The more massive the primary the greater the separation 
must be. We find less than 10% of systems exhibit this behaviour. For this to occur the total mass of the 
binary must be a few solar masses greater than the minimum mass for a single star SN. 

Such binaries are similar to those discussed by Davies, Ritter & King (2002) but these stars do not undergo 
a second interaction when the secondary evolves. The reason for this comes from angular momentum con- 
servation. Because the mass transfer is so efficient the angular momentum of the system is almost fully 



147 



5.3 EVOLUTION OF MASSIVE BINARIES 



conserved. In these systems, soon after the mass transfer occurs, enough mass has been accreted on to the 
secondary for it to be the more massive star. As the less massive primary now continues to transfer mass 
the orbit widens. Now the orbit is wider it prevents interaction between the secondary and the remnant of 
the primary. The secondary is now around 8M and evolves towards a type IIP SN. The main interest in 
these stars are that they appear to have a lower luminosity than single stars however upon closer inspection 
our models have failed just as second dredge-up is occurring. They are likely to continue and form extreme 
S-AGB stars. The helium core masses are lower than for single stars of the same mass because they were 
initially low-mass stars. However the CO and ONe cores are above Mch as are those of the extreme S-AGB 
stars mentioned in chapter 3. The details are listed in table 5.1 and all fit the observed values of SN2003gd 
at this time. However if they do experience 2nd dredge-up this will move them far outside of this region. 

5.3.2 9<M l /M Q <\5 

In this range we still find the two binary interactions during the evolution of the primary although by 13M 
they become rare. Primaries still interact by RLOF and CEE in roughly equal probability while for the 
secondaries there are still more RLOF events than CEE events. We also now find none of the secondary SNe 
discussed above because the secondaries of similar mass are massive enough that they would go SN without 
aid from their primary. From 1OM upwards the primaries always undergo some form of SN. If there has 
been interaction these are lb SN if not IIP. The band of SN progenitors where only a small fraction of mass 
has been removed for a IIL SN is again small. However this is likely to be the most common source of IIL 
SNe because the IMF favours low-mass stars. 

The lb SNe in this region are very different to the lb and Ic SNe seen from single stars. They are normally 
helium giants rather than dwarfs due to their lower mass. Their surface mass fractions of helium are over 0.9 
so this would easily be observed in the SN. It is possible therefore that there are two types of lb SN, giants and 
dwarfs that are less and more massive respectively. If we use the prescription of Heger et al. (2003) discussed 
in chapter 4 to determine how bright these SNe should be many of these lb SN are highly luminous, were 
with our preferred rates we do not see any bright Ibc SNe from single stars at solar metallicity. This can 
be understood from the fact that the smaller stars have more tightly bound envelopes than the larger giants. 
Therefore dwarfs eject less mass than a similar sized giant with the same explosion energy. More ejected 
matter will provide a more visible display. 

5.3.3 15<Mi/M Q <30 

We now start to see more CEE events than RLOF events in this range. We also see the end of the primary 
undergoing two interactions with the secondary because the helium stars formed from the first mass-loss 
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event are too massive to become giants. In contrast the secondary can still exhibit this behaviour. 

The total fraction of interactions now increases with even the very wide binaries, R > 10 3 R 0, now interacting 
because the massive giants become extremely large. Most interactions still lead to a lb or Ic SN. Very few 
binaries are only slightly stripped of their envelopes to leave hydrogen remaining and give rise to a IIL SN. 

Stars stripped of their hydrogen envelope in this region also give rise to WR stars and therefore increase the 
WR/O ratio, although only the most massive are likely to progress far along the WR evolution track to WC 
stars. For that the WR star must be massive enough to have a strong wind of its own; this will therefore 
mainly increase the fraction of WN stars. 



5.3.4 3O<Mi/M <lOO 

Stars in this range would normally become WR stars and mostly give rise to Ic SN if they were single stars. 
Most interactions take the form of CEE during the RGB stage of evolution which is when these stars would 
normally lose their envelopes in very strong winds. This makes evolution in this region highly uncertain. 
The combination of strong wind mass loss as well as CEE mass loss rapidly accelerates the formation of the 
WR star but does not change the final outcome. 

If the interactions occur just after the main sequence, evolution is altered to a greater degree because the 
helium core has only just formed. Removal of the envelope at this point gives rise to a slightly less massive 
WR star although the evolution still leads to the same end point. 

Evolution starts to change at 8OM when the mass loss becomes extreme and the time on the RGB is ex- 
tremely brief. This means the binaries must be closer to interact because the separation is increased by the 
loss of angular momentum in the stellar wind and the star does not grow to a large radius on the RGB. 



5.3.5 100 < Mi /M < 200 

The evolution of these stars is less spectacular. Mass loss is so great that, unless the stars have a separation 
below 1O(W?0, there is no interaction at all, even during the evolution of the secondaries! At the time of the 
primary SN the separation is normally greater than 1OOO7? . Therefore we know at solar metallicity that all 
the interesting binary SN progenitors must occur at masses less than 1OOM . That is still quite a large range 
to study. 
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Table 5.1: The final details of the IIP low luminosity and possible S-AGB progenitors. 



M. lM 
lvl\l 1VI 


M~ lM 


a i/ K Q 


\nc,(J It \ 


\c\ci(T .... IK\ 


lOg(K/K ) 


iw final/-" /j '0 


J 


A ^ 




A 




9 71 


S 1 

O. 1 


5 


4.5 


158.49 


4.46 


3.54 


2.67 


7.7 


6 


4.2 


158.49 


4.60 


3.53 


2.75 


8.2 


6 


4.2 


251.19 


4.52 


3.54 


2.71 


7.9 


6 


5.4 


251.19 


4.60 


3.54 


2.75 


8.7 


7 


4.9 


398.11 


4.46 


3.54 


2.68 


7.7 


7 


6.3 


398.11 


4.57 


3.54 


2.73 


8.4 


8 


7.2 


251.19 


4.51 


3.54 


2.69 


7.9 



5.4 Oddities 

We now discuss some of the oddities and outcomes of the above evolution that are important for SN progen- 
itors. 

5.4.1 Low Luminosity Progenitors 

As discussed in section 5.3.1 we did find some progenitors that did have final low luminosities but it is 
possible that they may actually undergo second dredge-up before core collapse. However this is not certain 
and will require more study of these systems. All these stars fit well with the observations of SN2003gd, 
the low luminosity red giant progenitor. One feature that is likely to affect these results of these models is 
the composition of the material we accrete onto the secondary. Currently we accrete material with the same 
composition as the star itself but in reality this would be slightly different and may even prevent second 
dredge-up. 

The only progenitors that were of significantly lower luminosity were those of lb SN where having removed 
the hydrogen envelopes. The removal of hydrogen burning as an energy source decreases the luminosity. 
None of these models could explain SN2003gd since it was a type IIP SN but they do have implications for 
observations of low mass lb progenitors. 

5.4.2 Type I progenitors 

Of main interest to type I progenitors is the fact that we find another class of distinct progenitors that have 
no hydrogen. The three types are as follows, 
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• Surface < 0.3, R/Rq < 10 and M > 6M & , from single stars M^m > 31M . 

• Surface > 0.3, R/Rq < 10 and M > 6M , from single stars 28 < M initial < 31M ( 

• Surface > 0.9, R/Rq > 10 and M < 6M , from single stars M in iti a i < 25M . 



There is probably a continuum of objects between the first two types but the third type are fairly homogeneous 
although the more massive they are the smaller their radius. Our dilemma is now how to assign these to the 
observed types of lb and Ic? This is difficult and anything we say is a guess. Many more detailed observations 
of type Ibc SNe are required to have a large enough sample to begin to distinguish the SNe from different 
progenitor types. Further the three types above do not include the more massive (M > 2OM ) helium-rich 
progenitors at lower metallicity. However they may be too massive to have a bright display, ejecting little 
mass, and will be difficult to observe. 

The guess we have made is that the main determinant is the fraction of helium in the surface since it is 
a good indication of the amount of helium in the progenitor. Therefore we stick with the lb, Ic and Ibc 
definitions we made in the previous chapters. SN1993J was a type lib SN and provides some evidence for 
this. The progenitor was a subgiant with a large fraction of helium in the envelope. In late times it had the 
appearance of a lb SN so take away the hydrogen and we move towards a typical lb SN where we need a 
large star and high surface helium fraction. However there is probably a continuum of objects and more 
detailed observations may indicate the observable differences in the SNe of these objects. 



5.4.3 Type II progenitors with stripped envelopes 

We find very few stars that retain hydrogen after interaction with a companion. Most retain only a tiny 
fraction of hydrogen. This could be due to having very strong limits on common envelope interaction. 
However this is necessary because the amount of hydrogen a star needs to become a giant is very small. 
Podsiadlowski, Joss & Hsu (1992) found that IIL SNe were most common when the mass ratio was close to 
1. The resolution of our study is quite poor and this has probably limited out ability to find these objects. 
Our total percentage of IIL/IIb SN however does reasonably match the results of Podsiadlowski, Joss & Hsu 
(1992) despite its approximate nature. In our calculations below we take IIL SN to be all those that lose the 
majority of their envelopes but retain a small fraction, such that Mh < 2M , and become type IIL SNe. We 
make no differentiation between IIL and lib with our data. 
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5.4.4 ULXs 

ULXs are ultra luminous X-ray objects observed in other galaxies. They appear to be compact objects 
accreting at super Eddington rates onto a stellar mass black hole. The initial suggestion was that they are 
intermediate mass black holes (IMBH) with masses 100 — 1000M Q . The problem with this is that there are 
not enough IMBHs in star forming regions where ULXs seem to be concentrated. King (2004) suggests that 
ULXs are normal X-ray binaries but are accreting at super-Eddington rates owing to thermal timescale mass 
transfer. By observing the Cartwheel galaxy they place limits on the system requirements for an ULX. These 
are that > 15M© and that, because the emission is thought to be observed along a jet axis, there must be 
many more of them, around 3000 in the Cartwheel galaxy. The perfect source would be massive stars in a 
CEE phase involving a black hole. 

CEE is known to occur with white dwarfs because there are many double helium dwarf systems known 
Nelemans & Tout (2004). It is hypothesised that neutron stars can also undergo CEE with main-sequence 
or giant stars and if the two coalesce a giant is formed with a central neutron star, Thorne-Zyfkow objects 
(Thorne & Z ytkow, 1975). However little is said about the interaction of a massive star and a black hole if a 
CEE occurs. When RLOF occurs on to a black hole an accretion disc which is hot enough to produce X-rays 
is formed . The more material the black hole accretes the greater the X-ray luminosity. These X-rays affect 
the structure of the donor star. Phillips & Podsiadlowski (2002) calculate the effect the accreting luminosity 
has on the shape of the donor. They show the structure is shrunk and the star is pushed away from the black 
hole. Material moves further out towards the second Lagranian point. Therefore it is feasible that black holes 
can avoid CEE because rapid mass transfer on to a black hole produces X-rays that prevent or limit the rate 
of mass transfer. This is probably the true limit to accretion on to black holes rather than the Eddington limit. 

Such systems are limited to those where the primaries form black holes initially. The stars must then be close 
enough after the first SN to interact. Such stars are found in the range 20 < Mi /M Q < 80. Less massive stars 
do not form black holes, more massive stars form systems that are too wide. At the low-mass end of this 
range few stars have a large enough secondary to fit the requirements of King (2004) while at the higher-mass 
end most of the stars interact but some become WR stars and thus the situation is further complicated. To 
fully study this region of possible systems a more accurate treatment of CEE and accretion on to black holes 
is required. 

5.4.5 GRBs 

GRBs are thought to be closely associated with SNe since SN1998bw and GRB980425 were observed to be 
so associated. There are a few others where it is thought the two events are also related. Binaries are the 
most likely progenitor for two reasons. First in the collapsar model, where a black hole is formed and an 
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accretion disc forms to fuel the black hole in its production of the jets that give rise to the observed display, 
the star must be rotating rapidly enough to form the accretion disc otherwise there is not enough angular 
momentum to form a disk and material accretes onto the central black hole spherically. Binary systems that 
have undergone a CEE phase fit the bill since the orbit is shifted to smaller radius and more rapid rotation if 
we assume the stars become tidally locked (Izzard, Ramires-Ruiz & Tout, 2004). 

The second reason is that if the progenitors are WR systems these are extremely dusty and dirty environments 
with a lot of circumstellar material from the WR wind. A standard WR wind is not observed in the afterglows. 
If we assume the GRB jets are produced along the rotation axis we can clear this material away but if not the 
gravity of the secondary star and the rotation may be necessary to clear a region through which the jet can 
shine. 

Importantly if we make the assumption that a black hole must form directly for a GRB to occur then we 
would not see any such systems at solar metallicity. At LMC metallicity we would see an initial mass of 
35M e or greater, while at SMC metallicity the mass limit drops further. Therefore it is likely that we need a 
low-metallicity environment for GRBs. This agrees with observations with the lowest redshift GRBs being 
GRB980425 with a redshift of z = 0.0085 andGRB031203 with a redshift of z = 0.1055. Then also the CEE 
evolution would be more common at low metallicity because the winds from the stars are lower so there 
would be more massive giants that can engulf their companions and undergo CEE to shrink their orbit and 
spin them up. 



5.5 Binary SN Population Results 



Rather than trying to fit individual systems with our binary models we look at the SN ratios that we predict. 
For each primary mass we work out how many of the systems give rise to SNe. The remainder give rise to 
white dwarfs or are lost in a merger event. To calculate the relative rates we use a few different mass ratio 
and separation distributions. We use three main mass distributions; the first is flat in q, whilst the others 
are based on the distribution of Hogeveen (1992) where the probability is flat when q < 0.3 and above this 
range the probability is proportional to q~ 2 . This therefore favours low-q systems and is used by de Donder 
& Vanbeveren (2003) and we denote it as Hog. It does however go against the observational evidence that 
most massive star binaries have components near equal mass. Although it will be difficult to detect low-mass 
secondaries to massive stars because by homology relations L <=c M 3 so a secondary with a tenth the mass it 
is three orders of magnitude fainter. To see the effect of biasing similar masses in binaries we use an Antihog 
distribution which is flat when q > 0.7 but proportional to (1 — q)~ 2 above this range. 

For separation we use a flat distribution in In a and compare this to that derived from Hogeveen (1992) where 
the probability of a separation is proportional to a ' . Again we use a Antihog distribution to favour large 
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separation systems. Combining all these systems we are able to draw some conclusions about how close we 
can get to the observed SN ratios. 

For comparison we use the observational details from Cappellaro et al. (1997) and Cappellaro (2004) that 
the current ratio of type II SNe to type Ibc SN is 5.00 ± 3.45 and that the rate of IIP to IIL SNe is the range 
roughly equal. However Gaskell (1992) suggests that IIL SNe could be less common as low as only a tenth 
of all type II SNe. We can also place some limits by counting the number of different SN that are have been 
observed. If we do this we find that there are roughly twice as many IIP SN than IIL. There are no firm limits 
on the ratio of type lb to Ic, mainly because the distinction between the two types is poorly defined. 

The results are given in table 5.2. We list the percentage of all binaries that undergo type II SN, type Ibc SN, 
merge or produce white dwarfs. We also list the ratios of some of the SN types. We combine all the possible 
separation and mass distributions to compare the affect of varying these. Tables 5.3 and 5.4 present the same 
information but for the primary and secondary star separately. It is possible to see the different trends from 
the two halves of the binaries. 

When we compare between the different distributions we see that the main effect is the separation distribu- 
tion. This is because with the separation we have more resolution and the range is much wider than for mass 
ratio. The mass ratio distribution however can provide fine tuning. A bias to high or low mass binaries only 
affects the ratios to a small degree. The largest impact is on the overall percentage of SN from binaries, as 
would be expected because if we bias to low-mass secondaries most of these will have too low a mass for a 
SN to occur. 

When comparing separation distributions we see that the flat and Hog distributions are similar in the ratios. 
Also the flat distribution has relatively more SN and the Il/Ibc ratio is close to unity. The Hog distribution 
removes a number of SN and also removes a large fraction of type II SN by biasing close binaries that 
interact. It also features more mergers which eliminates the second SNe from secondaries. 

The Il/Ibc ratio is perhaps the most interesting. Only the Antihog separation distribution results are within 
the error bars of the observed systems. This gives only one conclusion, single stars are very important, de 
Donder & Vanbeveren (2003) states that most massive stars are likely to be in binaries and thus single-star 
studies are only of academic interest. However if our models are to be believed this means that in fact we 
must include a relatively large fraction of single stars to produce what is observed. Thus massive single star 
evolution is of much interest. 

The Ib/Ic ratio in all systems is greater than one. This indicates that there are more lbs than Ics. This can 
be understood in that lbs are produced by the more populous low-mass stars while the Ics come from higher 
mass stars. Observations of this ratio are poorly constrained. Also massive type I progenitors might not have 
any observable display. The details of type I SNe are also confused by many different progenitors that fit the 
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requirements for these SN. 

The IIP/IIL ratio is the most worrying. For none of our systems can we get lower than a ratio of 2. 12. There 
are two solutions to this problem. Either there are fewer IIP SNe in our models than we think; some are Iln 
SN and some are Ilpec, so it is possible to reduce this ratio if we say around a quarter to a half of what we 
think are type IIPs are something else so our ratio is the ratio of all other type II SN to IIL SN. The value 
for this ratio Podsiadlowski, Joss & Hsu (1992) provide is 4.75. All are within the to \. One problem 
with observations are how complete are they? If there are low luminosity IIP SNe that form black holes they 
could be missed and thus decrease the fraction of IIL SNe by increasing the total type II population. The 
other less favourable scenario is to produce more type IILs from low mass S-AGB stars that go SN. 

To decide which schemes are the best is difficult. Those which match the observations best are those with a 
flat separation distribution in In a with a flat or Antihog mass ratio distribution. 



5.6 Conclusion 



We have produced some results from simple binary evolution calculations of massive stars. From the results 
we have made predictions of the percentage outcomes and SN ratios. These fail to match observation indicat- 
ing that we do need to include a non-negligible fraction of single stars to match observations. Also the ratio 
of IIP/IIL ratios is of particular interest to us. We also find only a few low luminosity binary SN progenitors 
that still have their hydrogen envelopes. Although there is uncertainty in the final evolution of these objects 
since they are possibly examples of extreme Super- AGB stars that may undergo second dredge-up before 
core collapse and therefore be very luminous prior to SN. 
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Table 5.2: Results for evolution of primary and secondary. 



1 Mi 


a/R Q 


SNII SNIbc Merge White 
1% 1% 1% Dwarfs/% 


Ib/Ic Mbc IIP/IIL 


Flat 
Hog 
Antihog 


Flat 
Flat 
Flat 


14.01 14.14 2.72 69.13 
11.76 12.17 3.8 72.27 
17.56 15.93 2.38 64.13 


3.62 0.99 3.43 
3.41 0.97 4.16 
3.55 1.1 3.09 


Flat 
Hog 
Antihog 


Hog 
Hog 
Hog 


6.37 18.52 7.67 67.43 
5.98 15.98 10.74 67.31 
8.73 20.75 6.77 63.74 


4.82 0.34 2.69 
4.47 0.37 4.95 
4.7 0.42 2.12 


Flat 
Hog 
Antihog 


Antihog 
Antihog 
Antihog 


20.79 8.31 0.45 70.46 
17.84 7.31 0.62 74.23 
24.45 9.22 0.38 65.95 


1.78 2.5 5.89 
1.71 2.44 6.45 
1.71 2.65 5.57 


Table 5.3: Results for evolution of primary. 


Mi 

i = m 


a/R Q 


SNII SNIbc Merge White 
/% /% /% Dwarfs/% 


Ib/Ic Mbc IIP/IIL 


Flat 
Hog 
Antihog 


Flat 
Flat 
Flat 


8.75 10.9 30.35 
9.22 10.3 30.48 
8.75 11.24 30.01 


3.37 0.8 3.87 
3.21 0.89 4.45 
3.34 0.78 3.59 


Flat 
Hog 
Antihog 


Hog 
Hog 
Hog 


3.7 14.6 31.7 
4.96 13.81 31.23 
3.59 14.9 31.51 


4.13 0.25 4.38 
4.05 0.36 6.41 
3.81 0.24 3.87 


Flat 
Hog 
Antihog 


Antihog 
Antihog 
Antihog 


14.2 6.39 29.41 
14.28 6.15 29.57 
14.25 6.53 29.23 


1.73 2.22 6.05 
1.65 2.32 6.51 
1.75 2.18 5.79 


Table 5.4: Results for evolution of secondary. 




a/R Q 


SNII SNIbc Merge White 
/% /% /% Dwarfs/% 


Ib/Ic Mbc IIP/IIL 


Flat 
Hog 
Antihog 


Flat 
Flat 
Flat 


5.27 3.24 2.72 38.77 
2.54 1.87 3.8 41.79 
8.81 4.69 2.38 34.12 


4.71 1.62 2.87 
4.91 1.36 3.33 
4.14 1.88 2.7 


Flat 
Hog 
Antihog 


Hog 
Hog 
Hog 


2.67 3.92 7.67 35.73 
1.01 2.17 10.74 36.08 
5.14 5.85 6.77 32.23 


10.75 0.68 1.58 
10.56 0.47 2.03 
9.79 0.88 1.51 


Flat 
Hog 
Antihog 


Antihog 
Antihog 
Antihog 


6.59 1.93 0.45 41.04 
3.56 1.16 0.62 44.67 
10.21 2.7 0.38 36.72 


1.93 3.42 5.58 
2.1 3.08 6.23 
1.63 3.79 5.28 
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Chapter 6 



Mixing the Progenitors 



"Typical human, you can always count on them to mess things up." 
The 7th Doctor Who, Remembrance of the Daleks. 



6.1 Introduction 

In this penultimate chapter we bring together the results from previous chapters to produce an overall picture 
of our work. First we add details on the lowest-mass progenitors, using results for S-AGB progenitors from 
chapter 3, to the results from chapter 4. On top of this data we plot details of observed progenitors to date 
and discuss the implications of our models for observations for the progenitors of SNe. Next we use all our 
single star data to discuss the implications for observation of progenitors and the mass distribution of the 
remnant black holes. 

Finally we combine the results from chapter 5 with the detailed single star models from chapter 3 to gain 
some idea on the relative populations of progenitors from single and binary stars, trying to fix the relative 
population to observed SN ratios. 



6.2 The full range of single star progenitors 

We have used the results for S-AGB SN progenitors from chapter 3 to extend our preferred grid from chapter 
4. We have not included mass loss in these models and have created models at every integer mass up to 
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5 6 7 8 9 10 15 20 25 30 35 40 50 60 70 80 1 00 1 20 1 50 200 5 6 7 8 910 15 20 25 30 35 40 50 60 70 80 100 120 150 200 

Initial Stellar Mass / M s Initial Stellar Mass / M e 



Figure 6.1: Grids for comparison with observations of SN progenitors. For all grids the thick solid white 
line is the minimum mass for SNe. The dashed line separates type II SNe with Mh > 2M from those with 
Mh < 2Mq. The dashed-dotted line separates type II and type I SNe. Top-left the contours are for pre-SN 
luminosity in log(L/L ). Top-right contours are for stellar type determined by surface temperature. Bottom- 
left the contours are for pre-SN radius in log(/?//? Q ). Bottom-right contours are for pre-SN mass of the 
progenitor in M G . The filled circles represent IIP SN, the circle with filled centre a IIL, the upright triangle 
Iln, the inverted triangle lib, the star symbol Ilpec and the square a Ic. The IIP SNe are from top to bottom 
1999gi, 1999em, 2001du and 2003gd. The IIL SN is 1980K. The lib SN is 1993J. The Ilpec is 1987A. The 
Ic SN is 2002ap. 
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Plateau length Sl\c Type 




5 6 7 8 9 10 15 20 25 30 35 40 50 60 70 80 100 120 150 200 5 6 7 8 9 10 15 20 25 30 35 40 50 60 70 80 100 120 150 200 

Initio! Stellar Moss / M e Initio! Stellar Moss / M e 



Figure 6.2: Grids for comparison with observations of SN progenitors. For all grids the thick solid white 
line is the minimum mass for SNe. The dashed line separates type II SNe with Mh > 2M & from those with 
Mh < 2M Q . The dashed-dotted line separates type II and type I SNe. Top-left the contours are the length of 
the plateau phase in days. The top-right panel shows the boundaries between SN types. The thick black lines 
in this grid represent the change from lb to Ibc to Ic SNe. Bottom-left the helium core mass. Bottom-right 
the final remnant mass, contours are in M . The filled circles represent IIP SN, the circle with filled centre a 
IIL, the upright triangle Iln, the inverted triangle lib, the star symbol Ilpec and the square a Ic. The IIP SNe 
are from top to bottom 1999gi, 1999em, 2001du and 2003gd. The IIL SN is 1980K. The lib SN is 1993J. 
The Ilpec is 1987A. The Ic SN is 2002ap. 
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1OM . We assume the maximum mass S-AGB star lies O.5M above the highest mass S-AGB model we 
have produced, 8.5M Q at Z = Z . The minimum mass S-AGB SN progenitor is then O.5M below our 
highest mass model, 7.5M Q at Z . The results are shown in figures 6.1 and 6.2. On the figures we have 
plotted the location of SN progenitors observed to date. The primary problem is that, in the type IIP region, 
we have SNe that are other type lis. There is evidence that two of these, 1987A and 1993J were in binary 
systems (Maund et al., 2004; Podsiadlowski, 1992). Therefore the remaining two SNe, 1980K and 1997bs, 
could also be binaries. 

SN1980K, as already discussed in chapter 3, indicates that IIL SNe do not necessarily come from massive 
stars that have been stripped of a large fraction of their hydrogen envelope but they could also come from 
S-AGB stars that go SN. A binary model with RLF or CEE would also produce a IIL SN. The problem 
would be fitting a binary to the change in the mass-loss rate over the short period of 10,000 years prior to the 
SN. However it may be possible to observe the companion of this system in the future to test this possibility. 
The progenitor would not have been a single massive star in the range of 15 < Mi n m a \/M < 20 that has lost 
a great amount of envelope before exploding because our models in that region have too great a hydrogen 
envelope remaining for such a SN. Furthermore our IIL single star progenitors might be classified instead as 
Iln SN as these stars have high mass loss leading up to their eventual SN event. This would make a dense 
circumstellar environment that could give the Iln appearance. Therefore the most massive single stars, that 
have only small amounts of hydrogen in their envelope, could be Iln progenitors not type IIL. Cappellaro et 
al. (1997) say that type Iln are about 5 — 10% of all observed type II progenitors. This certainly agrees with 
our models at solar metallicity. The range of models that fit this star are 23 to 26M . 

The type IIP SN progenitors are all in the region we expect IIP progenitors from our models and the plateau 
lengths are in reasonable agreement. However SN2003gd requires special attention as it draws our eye to 
the minimum mass for S-AGB SN progenitors. From the derived luminosity of the progenitor of 2003gd, 
log (L/L Q ) = 4.3 ± 0.3, we have a problem that our models only just go under the upper limit of 4.6. We can 
look at this from two views. The first is that we can place quite stringent limits on the mass of 2003gd since 
it must be around 8.5M otherwise it would have gone through second dredge-up and be a great deal more 
luminous. We could also say that the host galaxy must be more distant than other observations indicate. 

The more correct view we should take is that this progenitor does agree with our models although there is 
now a mystery that can only be solved by more observations of progenitors in this region. Second dredge- 
up before core collapse can be prevented by making convection less efficient at the base of the convective 
envelope. Therefore if 2003gd was around 8M and did not undergo second dredge-up this has important 
ramifications for stellar structure models and the details of convection. 

The only type Ic SN in the field, 2002ap, is in the right place. From the graph we can see that the star's mass 
is also tightly constrained. This is also the only SN that is definitely in the region where we expect a black 
hole was formed; the remnant should have been around 6M . However any position to the right of this is 
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also possible since there is only a limit on it mass from observations. However that would require a binary 
system to explain a mass lower than the theoretical Ic region. 

The super-star SN1987A, the best observed SN since the invention of the telescope, can be seen right on the 
limit of where it may or may not have formed a black hole by fall back on to a neutron star. A remnant has 
yet to be observed at 1987 A although it might have been ejected from the system and so not been discovered 
but there is a chance it could have been a black hole. However if it was the result of a binary the primary 
would have been a lower-mass star to begin with and so would be in the region that would not have formed 
a black hole, although this relies on our very approximate methods to decide on the remnant left by the SN. 

When observers find more SN progenitors we will be able to place them on this diagram and draw more 
conclusions. This will force us to refine our models over time changing the mass-loss rates, our details of 
convection and other physical details in the code. Only time will tell the quality of our current models. 

6.2.1 Remnant masses 

Figure 6.3 shows the logarithm of the fraction of remnants with mass M rem in the ranges listed. At the highest 
metallicities the remnants are all of a few M Q with neutron stars dominant. This means there is a maximum 
mass for remnants at lower metallicity. Then as metallicity decreases the maximum remnant mass increases 
until the production of black holes greater than 1OOM becomes possible. The remnant masses used for 
these diagrams come from the very simple calculation explained in chapter 3. They are at best estimates of 
the distribution. 

Very massive (M > 1OOM ) black holes can form in SNe but this can only happen at low metallicity (Z < 
10 ). The population of such black holes is very small as they are only formed in the most massive stars. 
At higher metallicities massive black holes are not formed because mass loss limits the size of the SN 
progenitor and thus the size of the remnant. Any black holes observed to be very massive in high metallicity 
environments must therefore have grown via accretion from a binary companion or must have formed from 
an earlier population of stars. 

6.2.2 HR Density Plots 

Figure 6.4 shows where we can expect single star progenitors to lie on the HR diagram. Red giants are 
the dominant group. If we decrease metallicity there are fewer WR star progenitors and the luminosities 
increase. The temperatures used here are from hydrostatic models. We know that the temperatures should 
be slightly lower but we can clearly see WR stars form a separate population from the red giants. 
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Figure 6.4: HR diagram progenitor population plots. Displaying the population over the HR diagram for 
single star SN progenitors. Contours are in % of total population. The yellow lines are the evolution tracks 
of 10, 20, 40 and 120 M Q stars at the metallicity of the contour plot. 
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Table 6.1: Type lb SNe progenitor details compared to the details for the a similar star that does not lose its 
hydrogen envelope. 



Initial Mass 

/M Q 


With Hydrogen Envelope 
log(r/K),log(L/L ) 


Without Hydrogen Envelope 
log(r/K),log(L/L Q ) 


11 


3.5,4.8 


4.0, 4.5 


12 


3.5,4.9 


4.3, 4.7 


13 


3.5,4.9 


4.4, 4.8 


15 


3.5,5.1 


4.6, 4.9 


20 


3.5, 5.3 


4.8,5.1 



We must ask the question where would binary progenitors lie on this diagram? Type II progenitors can be 
blue rather than red super giants or have their envelopes removed to become type Ibc progenitors. If the type 
II progenitors become blue super giants through a merger event or mass-gain near the end of their evolution 
this would move them over to the right of the diagram but the luminosity would remain the same because it 
is fixed by the hydrogen burning shell. The situation is quite different if the hydrogen envelope is removed 
completely because the primary source of a star's luminosity is the hydrogen burning shell. This event would 
produce a range of low-luminosity progenitors in between the two populations in the diagram. 

Table 6.1 gives the details of some example lb progenitors with equivalent results for single stars that retain 
their hydrogen envelope. This has important implications when we try to fit the masses of progenitors to type 
lb SNe. Stars that are in binaries that are greater than 20M© are probably massive enough that interactions 
lead to the helium stars becoming WR stars and therefore increase the WR population. It is not possible to use 
the standard stellar evolution models from the hydrogen main sequence alone to identify type I progenitors. 
A sequence of models from the helium main sequence should be constructed and evolved and then these 
used for comparison to observations of the progenitor. Once the helium star initial mass is known we would 
be able to then determine scenarios for the source of the star and its initial mass. 



6.3 Combining Single Stars and Binaries 



Now we take the SN type ratios for single stars from chapter 3 and the ratios for binary stars from chapter 
5 and mix them together. This should lead to a better agreement with what we see in observations. We 
use the same SN type classification scheme for both sets of data. For type II SNe we assume all stars with 
Mh < 2M Q are type IIL SNe and the rest are plain type II. We split the Ic and lb SNe to those with a surface 
mass fraction of helium either greater than or less than 0.5. 

To summarise the observations of Cappellaro et al. (1997) and Cappellaro (2004) we know that the II/Ibc 
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Table 6.2: Mixing different populations of single and binary stars. Not including a S-AGB star contribution 
to IIL SN. 



Single 


Binary 








Fraction 


Fraction 


Vb/lc 


Mbc 


II/IIL 





1 


3.62 


0.99 


3.43 


0.1 


0.9 


3.28 


1.70 


6.82 


0.2 


0.8 


2.94 


2.41 


10.21 


0.3 


0.7 


2.60 


3.12 


13.61 


0.4 


0.6 


2.25 


3.84 


17.00 


0.5 


0.5 


1.91 


4.55 


20.40 


0.6 


0.4 


1.57 


5.26 


23.78 


0.7 


0.3 


1.23 


5.97 


27.18 


0.8 


0.2 


0.89 


6.68 


30.57 


0.9 


0.1 


0.55 


7.39 


33.96 


1 





0.21 


8.10 


37.35 



Table 6.3: Mixing different populations of single and binary stars. Including a S-AGB star contribution to 
IIL SN. 



Single 


Binaiy 








Fraction 


Fraction 


Ib/Ic 


Mbc 


II/IIL 





1 


3.62 


0.99 


3.43 


0.1 


0.9 


3.28 


1.70 


4.65 


0.2 


0.8 


2.94 


2.41 


5.87 


0.3 


0.7 


2.60 


3.12 


7.09 


0.4 


0.6 


2.25 


3.84 


8.31 


0.5 


0.5 


1.91 


4.55 


9.53 


0.6 


0.4 


1.57 


5.26 


10.75 


0.7 


0.3 


1.23 


5.97 


11.97 


0.8 


0.2 


0.89 


6.68 


13.19 


0.9 


0.1 


0.55 


7.39 


14.41 


1 





0.21 


8.10 


15.63 
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ratio must be 5.00 ± 3.45 and the II/IIL ratio must be a few. However depending on galaxy type the II/Ibc 
ratio can take values from 3.82 ±2.71 to 6.14 ±3.96. This makes the situation very difficult to fit with 
models because the value is so uncertain. In fact ratios calculated from any of our mass-loss schemes when 
mixed with binaries agree with one of these values and indicate that binaries are very important for models 
of SNe ratios. They are much more important than mass-loss prescription or inclusion of extra mixing from 
rotation. 

When we look at table 6.2 and use the ratio of Cappellaro (2004) to determine the binary fraction we see that 
it should be around 45% but, due to the errors in the measurement, the errors in our estimation are ±40%. 
The expected binary fraction is expected to be higher than this and to not be constant (Mermilliod, 2001). 
If we assume a fraction of 70% we find a ratio of 3.12 which is on the low side of the observations. The 
binary fraction is likely to vary with initial mass and the distribution of mass ratio and separation. Also we 
can favour the lower ratios if we assume that some Ibc SNe, that we have determined occur, are not observed 
because they are too compact to eject much material and remain unseen. Because of this is it likely that we 
should favour the higher binary fractions. 

The type II/IIL ratio is a problem. From single stars we have only a tiny fraction of possible type ILL 
progenitors and some of these are quite likely to be Iln progenitors. If we include the lowest mass S-AGB, 
stars between 7.5 and 7.6M Q , as in table 6.3, we still find poor agreement. The other option is to lower the 
mass for massive single star IIL progenitors. However this doesn't produce a large change until it is lowered 
to below 18M S and this would produce too many WR stars for the WR ratios. 

There are just too few IIL SNe at the binary fraction that agrees with the II/Ibc ratio. We must get a larger 
population of these models. The main source must be binary stars. In fact the need to produce more IIL SNe 
has important implications for our models. We require a finer grid to find the exact region where IIL SN 
could occur. The other possibility is that CEE could be less efficient than we think and that not all hydrogen 
is removed in such an interaction. It is also possible, that during the CEE, hydrogen is mixed into the helium 
core by an unknown dynamical process that would leave enough hydrogen for a signal in the SNe spectra. 
However we must remember that Cappellaro et al. (1997) are very vague in their estimation of the relative 
IIL population so the error bars are extremely large. Gaskell (1992) suggest the IIL fraction of all type II 
SNe could be a little as one tenth. Our results indicate the later however we will need to wait for more 
observations to tie this ratio down. Also there could be selection effects in the observations that are not 
understood. 

The other solution would be to increase the mass-loss rates to produce more IIL progenitors from single 
stars but this would then reduce the agreement with the WR ratios and furthermore decrease the II/Ibc ratio 
that would in turn decrease the binary fraction needed to reproduce observed values. Therefore this problem 
must be addressed by further observations. 
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6.4 Conclusion 



In this chapter we have combined and analysed our data from the previous sections to draw some conclu- 
sions on the total population of SN progenitors and shown that it does fit observations well and that from 
observations we can provide vital constraints on stellar evolution models. 
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Chapter 7 



Conclusion 



"Answers are easy. It's asking the right questions which is hard." 
The 4th Doctor Who, The Face of Evil. 



7.1 Summary 

This dissertation is long and a final summary to highlight some of the important points is necessary. Also we 
should draw some final overall conclusions on the progenitors of SN progenitors. 

The preliminary chapters 1 and 2 provided an overview to stellar evolution and supernovae. Chapter 2 
described in detail the stellar evolution code used and presented the refinements made to the Cambridge 
code as part of this study. The most important of these was the construction of new and unique opacity tables 
that refine the models of stars during the latest stages of evolution. 

Later chapters covered the evolution of SN progenitors. Chapter 3 contains a detailed look at the lowest 
mass stars to go SNe, S-AGB stars undergoing ONe core collapse by electron capture. The exact nature of 
these stars can only be determined by observations of SN progenitors which can provide tough constraints on 
stellar physics, most critically on convection. This chapter also contained a high resolution study of single 
star progenitors at solar, LMC and SMC metallicities using our favoured mass-loss scheme. By using the 
results we made suggestions on how to relate these models to types of SNe, showing that there is much 
to study in the region of type I SNe with a diverse range of possible progenitors and that the difference 
between type lb and Ic SNe may be more than just the absence or presence of helium. For type II progenitors 
we discussed that IIL progenitors come not only from binary stars but from the S-AGB stars, while Iln 
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progenitors may come from single WN stars that still have hydrogen envelope but have lost a large fraction 
via severe stellar winds. 

In chapter 4 we have described the multiple data sets created to compare the diverse mass-loss schemes to 
fix on one that fits best to observations. We showed that the inclusion of metallicity scaling of mass loss 
for WR stars makes important changes to the population of SNe progenitors and that this scaling should be 
included. From arguments using observations we then settled on our favoured mass-loss prescription that 
uses the rates of de Jager, Nieuwenhuijzen & van der Hucht (1988), Vink, de Koter & Lamers (2001) and 
Nugis & Lamers (2000). 

Next in chapter 5 we moved on to study binary progenitors. After describing our method we discuss the 
results showing that the type I progenitors differ from those from single stars, adding further evidence to the 
idea that there are a few different populations of type I SNe that we do not currently distinguish. 

Finally in chapter 6 we collated some of our information to draw conclusions over the entire mass range 
for SNe and compare this directly to observations. We also then mixed together our results from single 
and binary stars to see that the evidence points toward the fact that the interacting binary fraction of stars is 
around 45% to best fit the observations. However the errors are large at ±40% and higher values are required 
to fit the II/IIL fraction. This shows that understanding of single and binary stars is of equal weight when 
we consider how the population of SN progenitors will be spread over the HR diagram and the distribution 
of black hole masses from supernovae. We importantly showed that black holes more massive than 100 M Q 
are only possible at metallicities below Z < 0.001. We also pointed out that the determination of the type I 
SN progenitors needs extra consideration to allow for the lower luminosity once the hydrogen envelope has 
been removed and that this could lead to overestimation of the initial masses. 



7.2 Conclusion 

Our main conclusions are to define what we can expect the progenitors of SNe to be, all the limits below are 
for solar metallicity with convective overshooting. For type II SNe they are, 

• IIP, these stars are probably single stars with a mass initially below about 25M . They have only lost 
a small fraction of their envelopes and are red giants when the SN occurs. At the low-mass end of this 
range, around 8 to 8.5M Q , they may have undergone second dredge-up shortly before core-collapse to 
become extreme S-AGB stars. There may also be a low luminosity class of IIPs where Min^ai > 20 M 
and a black hole is formed rather than a neutron star and affecting the SN evolution. 

• lib, are mostly likely to be from massive stars in a binary when interactions leave only a tiny fraction 



170 



7.2 CONCLUSION 



of the original hydrogen. It is also probable that this type of progenitor must be in the same mass range 
as progenitors of type lb SN. There is probably a smooth progression from IIL to lib to type lb SNe. 

• IIL, a large fraction of stars that experience only light mass loss via an interaction with a binary 
companion to lose some but not their entire envelope. The population of IIL SNe relative to all type 
II SNe places constraints on binary interactions and the binary mass ratio and separation distributions. 
A small fraction of the progenitors are massive single stars that have lost their envelope by stellar 
winds although the models in this class might be more likely to give rise to Iln. There is also a chance 
that low-mass stars in the tiny range 7.5 to 7.6M give rise to a IIL SNe rather than IIP because they 
experience second dredge-up a few thousand years before core-collapse. This is enough time for a 
wind to remove sufficient envelope to prevent a plateau phase in the light curve. The interesting thing 
is that all these progenitors must produce a similar amount of 56 Ni 0.1M q and all of these models 
do. Observations and further modelling of SNe are required to deal with this group. It is likely there 
is not a unique progenitor and observations already indicate possible subgroups in this class. 

• Iln, are mostly likely to be from hydrogen rich WR stars about 25-26 M Q for single stars. Such stars 
have a dense hydrogen rich circumstellar environment that the SN ejecta collides with to retard its 
expansion. Binary stars could also provide some of these stars from a RLOF or CEE event close to the 
occurrence of the SN. 

• Ilpec, well everything else really. We should introduce a new group of SN that have undergone mergers 
and therefore become blue stars such as SN1987A. Most type Ilpecs probably suffer from lack of 
detailed observations to specify what type they are accurately so many in this bin probably belong to 
one of the other groups. For example there are weak IIP SNe that are thought to form black holes 
rather than neutron stars that would get misplaced into this group. Only more detailed observations 
will provide more clues about this region. 

For type I SNe we should not provide a list as for type II SNe. The observational scheme does not provide 
enough possible groups to place our progenitors in. The answer to which stars give rise to which type I SN 
is difficult because we are trying to discover what these stars don't have pre-SN. The answer is probably not 
just to do with the presence or absence or material but also the size of the object, its structure and the nature 
of the remnant formed. 

However lb progenitors are probably from the evolution of binary stars. They are low-mass helium stars that 
become giants, similar to the progenitors of lib SNe but have been fully stripped of their hydrogen envelopes. 
Type Ic SNe are probably more compact stars because there is evidence that the ejecta velocities are higher 
for Ic SNe than for lb SNe so, if we assume a constant energy from the SN, then less mass is ejected but 
at a higher velocity. This would be correct if the star is more compact because the material would be more 
tightly bound. 
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Therefore our other compact progenitors are probably observationally all classed as type Ic. The more 
massive progenitors that form a black hole directly in their centre could be the progenitors to hypernovae but 
these stars would need to be rapidly rotating for this to occur. 

There are broadly three classes of progenitors that fit into the Ic group. The first are stars around solar 
metallicity initially more massive than 3OM with final masses less than 15 M e . They form black holes 
by fallback or directly. If the black hole forms directly this may lead to a different explosion mechanism 
and therefore different behaviour. The most massive progenitors in this class may be too tightly bound for 
enough matter to be ejected to provide a display. They also have less than a tenth of a solar mass of helium 
remaining so there would be little to affect the spectra of these objects. 

The second group are single stars at lower metallicity are more massive and retain more helium than their 
solar metallicity counterparts. However they are compact. This has an important effect on their SNe. They 
may only give rise to dim SNe or have no display at all because they are so tightly bound that very little 
material is ejected. If there is any display it is uncertain whether it would be possible to observe the helium. 

The final class are GRBs or hypernovae that fall into the type Ic scheme. These are probably in binary 
systems to ensure rapid pre-SN rotation. They must be more massive than 15M Q so a black hole is formed 
directly. The rotation ensures an accretion disc can form to power jets and give rise to an energetic SNe, 
hypernovae or GRBs. 



7.3 Future Work 

"A code is never fully tested until you run it for the last time." 
Sverre Aarseth, 2004. 

The main path forwards from here can be succinctly put; more models and more observations are required. 
While the more observations point relies on other researchers the more models is easy to do. With 70 years 
of computer time it is possible to calculate over half a million model stars. That is roughly the total evaluated 
for this dissertation. 

It is important to model the correct thing. Current plans are to produce a vast and detailed grid of model 
SNe progenitors over the range used in this thesis but with higher resolution in mass and metallicity than 
used here and to reach zero metallicity. However rather than making a single jump it would be better to 
lower the metallicity in steps to determine exactly when the behaviour alters. At the current time and for the 
foreseeable future there will be no observational constraints here. 
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The binary models also need refinements to produce the remnant and HR population diagrams that we have 
for single stars. That is the easy part. Waiting for the code to run a large grid is the problem. However time 
and multiple computers ease the problem. Using the code in parallel with binary population synthesis so we 
know the interesting places to study will also be necessary. 

Comparing binary models to limits for SN progenitors is difficult owing to the large parameter space, espe- 
cially for Ibc progenitors. Therefore it will be quite useful to produce model helium stars in a similar fashion 
to the models in chapters 3 and 4. This will be more practical than just using the normal tracks because 
evolution of low-mass stars, once they have lost their hydrogen envelopes, is significantly different. Once 
we know the initial helium star mass it will be easier to work out where it came from. 

Another step is to take all these models and feed then into a SN simulation code. It doesn't have to be 
a state of the art 3D magnetohydrodynamic masterpiece of coding either. An existing ID code will be 
sufficient. This will replace the need to use analytical approximations and provide further understanding of 
the difference between SNe. 

One important refinement is to take the burning stages of the progenitor models yet further, all the way to the 
formation of an iron core if possible. However this will require making the code stable at time steps below 
10~ 4 years by fixing the mesh and calculating the opacity of central regions on the fly to account for the 
changing composition. 

Also with this data looking at how the SN ratios change over time will be another useful analysis of the 
above data which may provide a better agreement to observations. 

Otherwise let's just hope for a few galactic SNe next year.... 



"Looks like the end of the line, space-mates!" 
Colonel Dan Dare - Eagle 4/7/1952 - Vol.3 No. 13. 
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"Time is an illusion. Lunchtime doubly so" 
Douglas Adams. 



"The true delight is in the finding out rather than in the knowing." 

Isaac Asimov. 

"I don't pretend we have all the answers. But the questions are certainly worth thinking about." 

Arthur C. Clark. 



